Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 1,850,958 T→A M1K (ATG→AAG) † relA → ATP:GTP 3'‑pyrophosphotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,850,9580TA99.2% 861.9 / NA 246M1K (ATG→AAG) †relAATP:GTP 3'‑pyrophosphotransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (120/124);  minor base C (0/2);  total (120/126)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 246 total aligned reads displayed.

CCTCAGACGGCGCATGTCGAGGCCATGGCGTTATTCGAAGCGGGCTAGCTGACTGGCCCCTTGGTGCGGTTGCCACGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGATGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAA  >  NC_002947/1850813‑1851102
                                                                                                                                                 |                                                                                                                                                
 cTCAGACGGCGCATGTCGAGGCCATGGCGTTATTCGAAGCGGGCTAGCTGACTGGCCCCTTGGTGCGGTTGCCACGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTAc                                                                                                                                             >  2:2844826/1‑150 (MQ=255)
 cTCAGACGGCGCATGTCGAGGCCATGGCGTTATTCGAAGCGGGCTAGCTGACTGGCCCCTTGGTGCGGTTGCCACGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTAc                                                                                                                                             >  1:5723228/1‑150 (MQ=255)
 cTCAGACGGCGCATGTCGAGGCCATGGCGTTATTCGAAGCGGGCTAGCTGACTGGCCCCTTGGTGCGGTTGCCACGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTAc                                                                                                                                             >  2:3169516/1‑150 (MQ=255)
 cTCAGACGGCGCATGTCGAGGCCATGGCGTTATTCGAAGCGGGCTAGCTGACTGGCCCCTTGGTGCGGTTGCCACGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTAc                                                                                                                                             >  2:5451819/1‑150 (MQ=255)
  tCAGACGGCGCATGTCGAGGCCATGGCGTTATTCGAAGCGGGCTAGCTGACTGGCCCCTTGGTGCGGTTGCCACGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACa                                                                                                                                            <  1:1373933/150‑1 (MQ=255)
  tCAGACGGCGCATGTCGAGGCCATGGCGTTATTCGAAGCGGGCTAGCTGACTGGCCCCTTGGTGCGGTTGCCACGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACa                                                                                                                                            >  2:3619122/1‑150 (MQ=255)
  tCAGACGGCGCATGTCGAGGCCATGGCGTTATTCGAAGCGGGCTAGCTGACTGGCCCCTTGGTGCGGTTGCCACGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACa                                                                                                                                            >  2:3787027/1‑150 (MQ=255)
  tCAGACGGCGCATGTCGAGGCCATGGCGTTATTCGAAGCGGGCTAGCTGACTGGCCCCTTGGTGCGGTTGCCACGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACa                                                                                                                                            >  1:325720/1‑150 (MQ=255)
   cAGACGGCGCATGTCGAGGCCATGGCGTTATTCGAAGCGGGCTAGCTGACTGGCCCCTTGGTGCGGTTGCCACGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAg                                                                                                                                           <  1:5291842/150‑1 (MQ=255)
   cAGACGGCGCATGTCGAGGCCATGGCGTTATTCGAAGCGGGCTAGCTGACTGGCCCCTTGGTGCGGTTGCCACGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAg                                                                                                                                           <  1:1282950/150‑1 (MQ=255)
   cAGACGGCGCATGTCGAGGCCATGGCGTTATTCGAAGCGGGCTAGCTGACTGGCCCCTTGGTGCGGTTGCCACGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAg                                                                                                                                           <  1:3631177/150‑1 (MQ=255)
     gACGGCGCATGTCGAGGCCATGGCGTTATTCGAAGCGGGCTAGCTGACTGGCCCCTTGGTGCGGTTGCCACGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGt                                                                                                                                         >  2:1103928/1‑150 (MQ=255)
         gcgcATGTCGAGGCCATGGCGTTATTCGAAGCGGGCTAGCTGACTGGCCCCTTGGTGCGGTTGCCACGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTgaga                                                                                                                                     <  1:2312010/150‑1 (MQ=255)
                tCGAGGCCATGGCGTTATTCGAAGCGGGCTAGCTGACTGGCCCCTTGGTGCGGTTGCCACGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCAcc                                                                                                                              >  2:1803585/1‑150 (MQ=255)
                tCGAGGCCATGGCGTTATTCGAAGCGGGCTAGCTGACTGGCCCCTTGGTGCGGTTGCCACGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCAcc                                                                                                                              <  1:5681096/150‑1 (MQ=255)
                     gCCATGGCGTTATTCGAAGCGGGCTAGCTGACTGGCCCCTTGGTGCGGTTGCCACGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCg                                                                                                                         <  2:3680680/150‑1 (MQ=255)
                     gCCATGGCGTTATTCGAAGCGGGCTAGCTGACTGGCCCCTTGGTGCGGTTGCCACGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCg                                                                                                                         <  1:4126230/150‑1 (MQ=255)
                       cATGGCGTTATTCGAAGCGGGCTAGCTGACTGGCCCCTTGGTGCGGTTGCCACGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGt                                                                                                                       <  2:5945883/150‑1 (MQ=255)
                        aTGGCGTTATTCGAAGCGGGCTAGCTGACTGGCCCCTTGGTGCGGTTGCCACGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTc                                                                                                                      <  1:4801427/150‑1 (MQ=255)
                          ggCGTTATTCGAAGCGGGCTAGCTGACTGGCCCCTTGGTGCGGTTGCCACGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCaa                                                                                                                    <  2:1353317/150‑1 (MQ=255)
                              ttATTCGAAGCGGGCTAGCTGACTGGCCCCTTGGTGCGGTTGCCACGGACGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACt                                                                                                                >  1:4601576/1‑150 (MQ=255)
                                aTTCGAAGCGGGCTAGCTGACTGGCCCCTTGGTGCGGTTGCCACGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGa                                                                                                              >  1:731101/1‑150 (MQ=255)
                                 ttCGAAGCGGGCTAGCTGACTGGCCCCTTGGTGCGGTTGCCACGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGAc                                                                                                             >  2:3244885/1‑150 (MQ=255)
                                 ttCGAAGCGGGCTAGCTGACTGGCCCCTTGGTGCGGTTGCCACGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGAc                                                                                                             <  2:4955762/150‑1 (MQ=255)
                                 ttCGAAGCGGGCTAGCTGACTGGCCCCTTGGTGCGGTTGCCACGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGAc                                                                                                             <  1:2345334/150‑1 (MQ=255)
                                     acgCGGGCTAGCTGACTGGCCCCTTGGTGCGGTTGCCACGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCa                                                                                                         <  2:1891031/148‑1 (MQ=255)
                                     aaGCGGGCTAGCTGACTGGCCCCTTGGTGCGGTTGCCACGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCa                                                                                                         <  1:5332037/150‑1 (MQ=255)
                                        cGGGCTAGCTGACTGGCCCCTTGGTGCGGTTGCCACGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTa                                                                                                      <  1:4038962/150‑1 (MQ=255)
                                          ggCTAGCTGACTGGCCCCTTGGTGCGGTTGCCACGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATc                                                                                                    >  1:343293/1‑150 (MQ=255)
                                            cTAGCTGACTGGCCCCTTGGTGCGGTTGCCACGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCaa                                                                                                  <  1:2536515/150‑1 (MQ=255)
                                             tAGCTGACTGGCCCCTTGGTGCGGTTGCCACGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACACAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAAt                                                                                                 <  1:3529902/150‑1 (MQ=255)
                                             tAGCTGACTGGCCCCTTGGTGCGGTTGCCACGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAAt                                                                                                 >  2:521485/1‑150 (MQ=255)
                                              aGCTGACTGGCCCCTTGGTGCGGTTGCCACGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGGGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAAtc                                                                                                >  2:4278186/1‑150 (MQ=255)
                                               gCTGACTGGCCCCTTGGTGCGGTTGCCACGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAAtct                                                                                               >  1:375384/1‑150 (MQ=255)
                                                 tGACTGGCCCCTTGGTGCGGTTGCCACGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCg                                                                                             >  1:2052548/1‑150 (MQ=255)
                                                  gACTGGCCCCTTGGTGCGGTTGCCACGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGa                                                                                            <  2:5771309/150‑1 (MQ=255)
                                                   aCTGGCCCCTTGGTGCGGTTGCCACGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGaa                                                                                           <  1:5591122/150‑1 (MQ=255)
                                                   aCTGGCCCCTTGGTGCGGTTGCCACGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGaa                                                                                           <  2:5081527/150‑1 (MQ=255)
                                                   aCTGGCCCCTTGGTGCGGTTGCCACGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGaa                                                                                           <  2:2825358/150‑1 (MQ=255)
                                                     tGGCCCCTTGGTGCGGTTGCCACGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGc                                                                                         >  2:1466184/1‑150 (MQ=255)
                                                          ccTTGGTGCGGTTGCCACGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGtt                                                                                   >  2:4982281/1‑150 (MQ=255)
                                                            ttGGTGCGGTTGCCACGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTg                                                                                  <  1:5246398/150‑1 (MQ=255)
                                                                tGCGGTTGCCACGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGtccaacc                                                                              >  1:2628093/1‑144 (MQ=255)
                                                                   ggTTGCCACGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCAtg                                                                           >  2:4182157/1‑150 (MQ=255)
                                                                     gtGCCACGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCAtgtg                                                                         <  2:2279381/149‑1 (MQ=255)
                                                                         cACGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTaa                                                                     >  1:5117923/1‑150 (MQ=255)
                                                                          aCGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGAATCAATCTCGAAGCATGGTTGGACCAtgt                                                                          >  2:4538285/1‑144 (MQ=255)
                                                                           cGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGc                                                                   <  2:2930997/150‑1 (MQ=255)
                                                                           cGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGc                                                                   <  1:3558700/150‑1 (MQ=255)
                                                                                 ggcggcCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGAt                                                             >  2:883364/1‑150 (MQ=255)
                                                                                     gcCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCgg                                                         <  1:4630536/150‑1 (MQ=255)
                                                                                          aCGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTg                                                    <  2:1742083/150‑1 (MQ=255)
                                                                                           cGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTgg                                                   >  2:766821/1‑150 (MQ=255)
                                                                                              tGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATc                                                <  1:4839338/150‑1 (MQ=255)
                                                                                                 tCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATcgcg                                             >  2:1048346/1‑150 (MQ=255)
                                                                                                     cAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAgc                                         >  2:1151707/1‑150 (MQ=255)
                                                                                                        cacaCTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCt                                      <  2:4126746/150‑1 (MQ=255)
                                                                                                         acacTCGTATGGTGTGCTGCATGGACAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTg                                     >  2:4618881/1‑150 (MQ=255)
                                                                                                           acTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGaa                                   >  1:1935028/1‑150 (MQ=255)
                                                                                                           acTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGaa                                   <  2:3871653/150‑1 (MQ=255)
                                                                                                           acTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGaa                                   <  1:3826184/150‑1 (MQ=255)
                                                                                                             tCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAaag                                 >  2:5652226/1‑150 (MQ=255)
                                                                                                               gTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAaagaa                               <  2:5282182/150‑1 (MQ=255)
                                                                                                               gTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAaagaa                               <  1:2719765/150‑1 (MQ=255)
                                                                                                               gTATGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAaagaa                               <  2:1891839/150‑1 (MQ=255)
                                                                                                                 aTGGTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGc                                                                   <  1:5858278/112‑1 (MQ=255)
                                                                                                                   ggTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCt                           >  1:1446323/1‑150 (MQ=255)
                                                                                                                   ggTGTGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCt                           <  1:1052945/150‑1 (MQ=255)
                                                                                                                    gtgtGCTGCATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTg                          >  2:3126525/1‑150 (MQ=255)
                                                                                                                       tgctgcATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGa                       <  1:1208607/150‑1 (MQ=255)
                                                                                                                       tgctgcATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGa                       <  1:4474425/150‑1 (MQ=255)
                                                                                                                           gcATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGttt                   >  2:4101323/1‑150 (MQ=255)
                                                                                                                            cATGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTg                  <  2:1404244/150‑1 (MQ=255)
                                                                                                                             aTGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGc                 <  1:3506686/150‑1 (MQ=255)
                                                                                                                              tGGAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGct                >  2:2509300/1‑150 (MQ=255)
                                                                                                                               ggAAATGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACTCTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGctc               <  1:1803573/150‑1 (MQ=255)
                                                                                                                               ggAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGctc               >  1:4264856/1‑150 (MQ=255)
                                                                                                                               ggAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGctc               >  2:4008350/1‑150 (MQ=255)
                                                                                                                               ggAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGctc               <  2:2411926/150‑1 (MQ=255)
                                                                                                                                gAAAGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGctct              >  1:1205829/1‑150 (MQ=255)
                                                                                                                                 aaaGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACt                                                                                                                <  2:50501/51‑1 (MQ=255)
                                                                                                                                   aGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGa           <  2:351456/150‑1 (MQ=255)
                                                                                                                                   aGGTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGa           <  1:2844826/150‑1 (MQ=255)
                                                                                                                                     gTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAgg         >  2:2631171/1‑150 (MQ=255)
                                                                                                                                     gTAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAgg         >  1:5146179/1‑150 (MQ=255)
                                                                                                                                      tAAAACAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGt        >  2:591910/1‑150 (MQ=255)
                                                                                                                                      tAAAACAAAGAAGGTACAGGTGAGAGGGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGt        <  1:2975763/150‑1 (MQ=255)
                                                                                                                                       aaaaCAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTa       <  2:4220033/150‑1 (MQ=255)
                                                                                                                                       aaaaCAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTa       <  1:4889679/150‑1 (MQ=255)
                                                                                                                                       aaaaCAAAGAAGCTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTa       <  1:4307080/150‑1 (MQ=255)
                                                                                                                                        aaaCAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTag      <  2:1718856/150‑1 (MQ=255)
                                                                                                                                          aCAACGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTagag    <  1:3911860/150‑1 (MQ=255)
                                                                                                                                          aCAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGa           >  2:525407/1‑143 (MQ=255)
                                                                                                                                          aCAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGa           <  1:525407/143‑1 (MQ=255)
                                                                                                                                           cAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAgaga   <  1:5012673/150‑1 (MQ=255)
                                                                                                                                           cAAAGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGGATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAgaga   <  1:5451819/150‑1 (MQ=255)
                                                                                                                                            aaaGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTCGAGaa  >  2:3273679/1‑150 (MQ=255)
                                                                                                                                            aaaGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGaa  >  2:5075184/1‑150 (MQ=255)
                                                                                                                                            aaaGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGaa  <  1:3210733/150‑1 (MQ=255)
                                                                                                                                            aaaGAAGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGaa  >  1:4960071/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CCTCAGACGGCGCATGTCGAGGCCATGGCGTTATTCGAAGCGGGCTAGCTGACTGGCCCCTTGGTGCGGTTGCCACGGAAGGGCGGCCGCACGGTGATCGCCAGCACACTCGTATGGTGTGCTGCATGGAAAGGTAAAACAAAGATGGTACAGGTGAGAGTGCACCAGCCGGTCAACACTGACGGCAGTATCAATCTCGAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAA  >  NC_002947/1850813‑1851102

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