Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,626,586 T→C G225G (GGT→GGC radA → DNA recombination protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,626,5860TC95.2% 162.6 / ‑4.3 63G225G (GGT→GGCradADNA recombination protein
Reads supporting (aligned to +/- strand):  ref base T (0/3);  new base C (25/34);  total (26/37)
Fisher's exact test for biased strand distribution p-value = 2.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACG  >  NC_000913/4626350‑4626832
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                          
ccGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGGTGGACACGTAACCAAAGATGGCTCGCTGGCGGGccc                                                                                                                                                                                                                                              >  1:359096/1‑251 (MQ=255)
              cccgtccgatatggAAACCAGCATCGTGCAGAGCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGAGGGACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAac                                                                                                                                                                                                                               <  2:200376/239‑1 (MQ=255)
                       gtCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATc                                                                                                                                                                                                                      <  1:111345/250‑1 (MQ=255)
                          ggAAACCAGCATCGCGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTTAAGGGCTGGAacagaacgggag                                                                                                                                                                                                                     >  1:200376/1‑240 (MQ=255)
                          ggAAACCAGCATCGAGCACATCTGCCTGTGTGCAGAAGAAGAGCAACCGAAGCTGAGGGTAATTGACGCGACCCAGTGGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTTAAACTGCGGCTTATCTTACGCGCTTCCCCAAAACGCGCGGTGTGGCGAGTGTCATGGTTGGACTCGTAACCAAAGATGGCTCGCTGGCTGGCCCTAAAGTGCTGGAACACTGTATCgcct                                                                                                                                                                                                                  <  2:425107/251‑4 (MQ=255)
                                  gCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTg                                                                                                                                                                                                          <  1:338391/251‑1 (MQ=255)
                                    aTCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGGGATGCAAATGGCGGATGGACAGTCATAGCCTGGCAGCGTGGCGCAGGGGCGGGAAACGGAGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGAGGGGCGATTGTCATGGGTGGAAACGTAACCAAAGAGGGCTCGCTGGCGGGCCCTAAAGAGCTGGAAAACTGTATCGACAGTTCGGGGCt                                                                                                                                                                                                        >  1:176099/1‑251 (MQ=255)
                                            gATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTa                                                                                                                                                                                                                        <  2:235181/227‑1 (MQ=255)
                                            gATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTa                                                                                                                                                                                                                        >  1:235181/1‑227 (MQ=255)
                                                 gCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGCCTGTTCGGTGCTTTTGGATGGCg                                                                                                                                                                                            >  1:76635/1‑249 (MQ=255)
                                                    agATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGcc                                                                                                                                                                                       <  1:40660/250‑1 (MQ=255)
                                                      aTTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTAACGCATGTTCGGTGCTTTTGGATGGCGATGCCgg                                                                                                                                                                                     >  1:347776/1‑250 (MQ=255)
                                                                 agagCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTc                                                                                                                                                                          <  1:148646/251‑1 (MQ=255)
                                                                  gagCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCg                                                                                                                                                                         <  1:447822/251‑1 (MQ=255)
                                                                    ccaccgtaaGCTGATGGTAATTGACGCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTACCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCa                                                                                                                                                                       <  1:127369/244‑1 (MQ=255)
                                                                        ccGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGGGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGGTTTTGGATGGCGATGCCGAATCCCGGTTTCGCActtt                                                                                                                                                                   >  1:423422/1‑248 (MQ=255)
                                                                          gAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTTGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTGCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCActttac                                                                                                                                                                 <  1:104364/251‑6 (MQ=255)
                                                                            aGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCCAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCCGGCGGGGCCCAAAATGCTGGAAAACTGTATCGCCTTTTTGGTGCTTTTTGGTGGCGATGCCGACTCCCCTTTTCGCActtttcgc                                                                                                                                                               >  2:144475/1‑244 (MQ=255)
                                                                                gATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGcc                                                                                                                                                           <  1:380788/251‑1 (MQ=255)
                                                                                     tAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGGTTGTCATGGTTGGCCACGTAACCAAAGATGGCTCGCGGGCGGGCCCTAAAGTGCTGGAACACTGTATAGACTGTTTGGGGGGTTTGGGTGGGGATGCCGACTCCCCtttttgccctttaccccccccaca                                                                                                                                                       >  2:9160/1‑230 (MQ=255)
                                                                                          gACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAAcc                                                                                                                                                  <  2:305367/250‑1 (MQ=255)
                                                                                                      gTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGt                                                                                                                                     <  1:298075/251‑1 (MQ=255)
                                                                                                      gTGATGCATATGGCGGATGTACAGACATCGCCTGGCAGCGGGGAGCAGGTGCGGGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACCCGTAACCAAAGATGGCCCGCTGGCGGGCCCTTAAGTGCTGGGACACTGGTACGCATGTTCGGTGCCTTTGGAAGGGCGATCAGCCTCCCGTTTTCCCAACTTACGCCGCCATAAAAACCCCttaggggaggg                                                                                                                                     >  2:396544/1‑242 (MQ=255)
                                                                                                         aTGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGaa                                                                                                                                  <  1:35723/251‑1 (MQ=255)
                                                                                                           gCATAGGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATg                                                                                                                                <  2:372209/251‑1 (MQ=255)
                                                                                                             atatGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACCCGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGGATCGACTGTTCGGGGCTATTGAATGGCGTTGCCGACTCCCGGTTTCGCACTTTACGCAGCCATAAAAACCGCTTCAGGGCGGTGAATg                                                                                                                                >  1:457861/1‑249 (MQ=255)
                                                                                                               atGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGGGTGGCGATTGTCATGGTTGGAAACGTAACCAAAGATGGCTCGCTGGGGGGCCCTAAAGAGGGGGAGCACTGTTTCGACTGGTCGGGGCGTTTGGATGGGGAGGCCGGCTCCCGTTTTTGCAATTTACGCCGCCCTAAAAACCCCTTCgggccgggggatgggat                                                                                                                            >  2:292824/1‑236 (MQ=255)
                                                                                                                   cGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGGGCTGGAACCCTGGATCGACTGTTCGGTGCTTTTGGGGGGCGGTGCCGACCCCCGTTTTCGCACTTTACGCAGCCATAAAAACCggctcgggcgggggaatgggctgggg                                                                                                                        >  2:45882/1‑226 (MQ=255)
                                                                                                                          accgTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTAACTGACTCGCTTCTCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGTTCGCTGGCGGGCCCTAAAGTTCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATGAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTcg                                                                                                                 <  2:176099/248‑1 (MQ=255)
                                                                                                                                     cTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGa                                                                                                      >  1:51192/1‑250 (MQ=255)
                                                                                                                                            gggggCGCTGGGGCGTGAAACGGTGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTTGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGc                                                                                              <  2:16368/248‑1 (MQ=255)
                                                                                                                                              tGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATAGTCATGGATGGACACGTAACCAAAGATGGGTCGCTGGCGGGCCCTATAGTGCTGGACCACTGTATCGACTGGTCGGGGCTTTTGGATGGCGATGCCTACTCCCCGTTTCTCACTTTTCGCAGCCATAAAAAACGCTTCGGCGCGGTGGATGAGACGTGCGACTTCGCCTAGGACCGAACAgggctg                                                                                             >  2:293936/1‑248 (MQ=255)
                                                                                                                                                  gcAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTg                                                                                        >  1:137365/1‑251 (MQ=255)
                                                                                                                                                     ggTGCATGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGATTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAg                                                                                     <  1:9771/251‑1 (MQ=255)
                                                                                                                                                         cGTGAAACGGCGGCTTATCGTACGCGCTTCGCCAAATCGCGCGGTGTGTCTATTTTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTGCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCTGCCATAACAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAg                                                                                 <  1:38551/251‑1 (MQ=255)
                                                                                                                                                           tGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGACCACTGTATCGCCTGTTCGGTGCTTTTGGGGGGCGATGCCGACTCCCGTTTTCGCCCTTTACGCAGCCATAAAAACCGCTTCGGCGCGGGGGCTGAGCTGGGCGTCTTCGCGATGCCCGAACAGGGGCTGCGTggaggcaac                                                                                >  2:359403/1‑242 (MQ=255)
                                                                                                                                                           tGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCa                                                                                  <  2:197324/248‑1 (MQ=255)
                                                                                                                                                           tGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCa                                                                                  >  1:197324/1‑248 (MQ=255)
                                                                                                                                                                 ggcggcTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGTAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCtt                                                                         <  1:175006/251‑1 (MQ=255)
                                                                                                                                                                  gtggcttacctGACGCGCCTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTTGACACGTAACCAAAGATGGCGCGCTGGCGTGCCCTAAAGTGCTGGAACACTGGATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTTAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTc                                                                        <  2:293833/241‑1 (MQ=255)
                                                                                                                                                                  gcggcTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTc                                                                        <  1:200037/251‑1 (MQ=255)
                                                                                                                                                                    ggtttagctGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTGGCGGGCGGGCCCTAAAGTGCTGGAACACGGTATCGACTGTGCGGTGCTTTTTGATGGCGATGCCGACTCCCGTTGTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTGCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCgg                                                                      <  2:162350/243‑1 (MQ=255)
                                                                                                                                                                             tacgcgattCGCCAAAACTCGCGGTGGGGCGATTGTCATGGGTGGATACGGAACCAAAGATGGTCCGCTGGTGGGCCGTAAAGTGCTGGACCACTGTATCGACTGTACGTTGTTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGTCCGAACAGGGTCTGCGTGAATACAGCAATATTTCGGCAATTTTc                                                              <  2:394924/243‑1 (MQ=255)
                                                                                                                                                                             gacgcGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAAtt                                                                 <  1:241998/243‑1 (MQ=255)
                                                                                                                                                                                gcgcgtCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGGTGGGGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCGTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTaa                                                          <  1:61136/246‑1 (MQ=255)
                                                                                                                                                                                 cGCTTCGCAAAAACGCGGGGTGTGGCGATTGTCATGGTGGGGCGCGTACCCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACGCCCGTTTTCGCACCTTGCGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAg                                                         <  2:114656/251‑1 (MQ=255)
                                                                                                                                                                                   cTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGCCTGCTCGGTGCCTTTGGATGGCGCTGCCGACTACCGTTTTCGCCATTTACGCACACATAAAAACCACTTCGGCGCGGGTAATGAGCAGGACTGCTTCGCGATAGCCCACCCGGGGCTGCGTGATGTAAGACACACTTAGGCAATTTTCTTTAGTc                                                       >  2:405223/1‑251 (MQ=255)
                                                                                                                                                                                     tCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTc                                                       <  1:144475/249‑1 (MQ=255)
                                                                                                                                                                                        ccAAAACGCGCGGTGTGGCGATTGTAATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCgcgga                                                  <  2:331747/251‑2 (MQ=255)
                                                                                                                                                                                          aaaaCGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGCTGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCCCGATGGCCGAACAGGGGCTCCGTGAAGTCAGCAACCCTTCGGACATTTTTTTAAGTCgcgggg                                                 >  2:389398/1‑248 (MQ=255)
                                                                                                                                                                                               ggggggTGTGGGGATTGTCATGGTTGGACACGGAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAg                                           <  2:9683/247‑1 (MQ=255)
                                                                                                                                                                                               gggcgGTGTGGCGATTGTCAGGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCTCTAAAGTGCTGGAACACTGTATCGACGGTTCGGTGCTTTTGGAGGGCGATGCCGACTCCCGTTTGCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAg                                           <  2:469532/249‑1 (MQ=255)
                                                                                                                                                                                                cgcgGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAgg                                          >  1:413023/1‑250 (MQ=255)
                                                                                                                                                                                                cgcgGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGAGGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGGACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGt                                          <  1:32672/251‑1 (MQ=255)
                                                                                                                                                                                                   ggTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGAc                                       <  1:221702/251‑1 (MQ=255)
                                                                                                                                                                                                            aTTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGGGCTGGAACCCTGGATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCCGCCATAAAAAACGCTTCGGCGCGGTGAATGGGCTGGGCGTCTTCGCGAAGACCGAAAAGGGGCTGCGGGACGTCAGCAACCCTTCGGCAATTTTCTTAAGGCgcggagtgggagggccctccggccg                              >  2:414251/1‑230 (MQ=255)
                                                                                                                                                                                                              tGTCATGGTTGGACACGTAACCAACGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTGTCTTAATTCGCGGAGATGAAGTGACCTCCGGCAGc                             <  1:483537/250‑1 (MQ=255)
                                                                                                                                                                                                              tGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGa                                        <  2:102891/239‑1 (MQ=255)
                                                                                                                                                                                                              tGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGa                                        >  1:102891/1‑239 (MQ=255)
                                                                                                                                                                                                                    ggttggacACGTAACCAAAGATGGCTCGCTGGCGGGCCCTCAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGGGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGGGATGAAGTGCCCTCCGGCAGCTCCGTg                       >  2:473478/8‑250 (MQ=255)
                                                                                                                                                                                                                    ggttggacACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTAGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGCCTTCGCGATGACCGAACCGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCCTAAGTCGCGGAAGTGAAGTGACCTCCGGCAGCTccgtt                       >  2:404505/8‑246 (MQ=255)
                                                                                                                                                                                                                        ggacACGTAACCAAAGATGGCTCGGTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTGAATTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGTtggt                  <  1:130828/248‑1 (MQ=255)
                                                                                                                                                                                                                                        aTGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACAATGTATCGCCTGTTCGGAGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCCTAAAAACCGCTTCGACGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGGTGGAGGTGTGGGAGGGAAcg  >  1:425064/1‑251 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                          
CCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACG  >  NC_000913/4626350‑4626832

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: