Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,637,234 C→A 7.7% K41N (AAG→AAT) ‡ der ← GTPase; multicopy suppressor of ftsJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,637,2340CA 7.7% 212.0 / 6.5 77K41N (AAG→AAT) ‡derGTPase; multicopy suppressor of ftsJ
Reads supporting (aligned to +/- strand):  ref base C (39/32);  new base A (4/2);  total (43/34)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.97e-01

GCATCCACCATAAACAGTACGACGTCCGCTTCTTCAATCGCCAGCAGCGACTGTTCCGCCATGCGGGTTTCTACACCGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCAT  >  NC_000913/2637085‑2637383
                                                                                                                                                     |                                                                                                                                                     
gcATCCACCATAAACAGTACGACGGCCGCTTCTTCAATCGCCAGCAGCGACTGTTCCGCCATGCGGGTTTCTACACCGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTAc                                                                                                                                                       <  2:1181107/150‑1 (MQ=255)
 ccacccccaTAAACCGTCCGACGGCCGGTCCTACAATCGCCAGCAGCGCCTGTTCCGCCTTGCGGGCTTCAAACCCGGCTTCTGTGCCATCAATCCCGCGCGAATCAATACAGATAACCTCCCGGCCTTCAATTTCCCCACGACCGTACt                                                                                                                                                      <  2:1181448/144‑1 (MQ=255)
   tCCACCATAAACAGTACGACGTCCGCTTCTTCAATCGCCAGCAGCGACTGTTCCGCCATGCGGGTTTCTACACCGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCCCGACCGTACTTa                                                                                                                                                    >  2:290560/1‑150 (MQ=255)
    ccaccaTAAACAGTACGACGTCCGCTTCTTCAATCGCCAGCAGCGACTGTTCCGCCATGCGGGTTTCTACACCGTCTTCTGTGCCATCAATCCCGCAGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTAc                                                                                                                                                   >  1:37812/1‑150 (MQ=255)
     caccaTAAACAGTACGACGTCCGCTTCTTCAATCGCCAGCAGCGACTGTTCCGCCATGCGGGTTTCTACACCGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACg                                                                                                                                                  <  2:26074/150‑1 (MQ=255)
        caTAAACAGTACGACGTCCGCTTCTTCAATCGCCAGCAGCGACTGTTCCGCCATGCGGGTTTCTACACCGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAAGTTCCGCCCGACCCTACTTACGGTc                                                                                                                                               >  1:1278000/1‑150 (MQ=255)
        caTAAACAGTACGACGTCCGCTTCTTCAATCGCCAGCAGCGACTGTTCCGACATGCGGGTTTATACACCGTCTTATGTGCCATCAATCCCGCCGGTATCAATACAGATAACCTCACGGGATTCAATTTCCGCACGACCGTCCTTACGGTc                                                                                                                                               >  1:1283049/1‑150 (MQ=255)
          tAAACAGTACGACGTCCGCTTCTTCAATCGCCAGCAGCGACTGTTCCGCCATGCGGGTTTCTACACCGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCAc                                                                                                                                             <  2:656660/150‑1 (MQ=255)
           aaaCAGTACGACGTCCGCTTCTTCAATCGCCAGCAGCGACTGTTCCGCCATGCGGGTTTCTACACCGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACg                                                                                                                                            <  1:1159282/150‑1 (MQ=255)
            aaCAGTACGACGTCCGCTTCTTCAATCGCCAGCAGCGACTGTTCCGCCATGCGGGTTTCTACACCGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTAAGTACGGTcgtgc                                                                                                                                           <  1:475277/150‑5 (MQ=255)
             ccAGTACGACGTCCGCTTCTTCAATCGCCAGCAGCGACTGTTCCGCCATGCGGGTTTCTACACCGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTAAGTACGGTcgtgcg                                                                                                                                          <  2:1214354/149‑6 (MQ=255)
              cAGTACGACGTCCGCTTCTTCAATCGCCAGCAGCGACTGTTCCGCCATGCGGGTTTCTACACCGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACCGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGt                                                                                                                                         <  1:494794/150‑1 (MQ=255)
                gTACGACGTCCGCTTCTTCAATCGCCAGCAGCGACTGTTCCGCCATGCGGGTTTCTACACCGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCa                                                                                                                                       <  1:1136656/150‑1 (MQ=255)
                    gacgTCCGCTTCTTCAATCGCCAGCAGCGACTGTTCCGCCATGCGGGTTTCTACACCGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGAcc                                                                                                                                   <  1:24928/150‑1 (MQ=255)
                      cgTCCGCTTCTTCAATCGCCAGCAGCGACTGTTCCGCCATGCGGGTTTCTACACCGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACACg                                                                                                                                 >  2:1258197/1‑150 (MQ=255)
                         ccGCTTCTTCAATCGCCCGCAGCGACGGTTCCGCCATGCGGGTTTCTACACCGTATTCTGTGCCATCAATGCCGCCGGTATCAATACAGATAACCTCACGGCCTTCAATTACCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGa                                                                                                                              <  2:57628/150‑1 (MQ=255)
                           gCTTCTTCAATCGCCAGCAGCGACTGTTCCGCCATGCGGGTTTCTACACCGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGaaa                                                                                                                            >  2:279612/1‑150 (MQ=255)
                             ttcttcAATCGCCAGCAGCGACTGTTCGGCCATGCGGGTTTCTACACCGTCTTCTGTGCCATCAATCCCGCCGGTATCAATAAAGATAAACTCCCGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATc                                                                                                                          <  2:932582/150‑1 (MQ=255)
                              tcttcAATCGCCAGCAGCGACTGTTCCGCCATGCGGGTTTCTACACCGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCCCGGCCTTCAAGTTCCGCACGACCGAACTTACGGTCACGAGTCAGACCCGGGAAATcc                                                                                                                         >  1:817750/1‑150 (MQ=255)
                                 tcAATCGCCAGCAGCGACTGTTCCGCCATGCGGGTTTCTACACCGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCa                                                                                                                      >  2:780128/1‑150 (MQ=255)
                                    aTCGCCAGCAGCGACTGTTCCGCCATGCGGGTTTCTACACCGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAAcc                                                                                                                   <  1:173842/150‑1 (MQ=255)
                                       gCCAGCAGCGACTGTTCCGCCATGCGGGTTTCTACACCGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGCTAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAATTCCGCACCCAgc                                                                                                                >  2:60286/1‑150 (MQ=255)
                                          agcagcGACTGTTCCGCCATGCGGGTTTCTACACCGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGCCCCGGGAAATCCGCAACCAGCGCa                                                                                                             >  1:678259/1‑150 (MQ=255)
                                          agcagcGACTGTTCCGCCATGCGGGTTTCTACACCGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCa                                                                                                             <  2:1216810/150‑1 (MQ=255)
                                                 aCTGTTCCGCCATGCGGGTTTCTACACCGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATACACTCCCGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCggg                                                                                                      >  1:1242205/1‑150 (MQ=255)
                                                   tGTTCCGCCATGCGGGTTTCTACACCGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTg                                                                                                    <  1:712544/150‑1 (MQ=255)
                                                    gTTCCGCCATGCGGGTTTCTACACCGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACGGTACTTACGGTCACGAGTCAGCCCCGGGAAATCCGCAGCCAGCGCATCTCGGGTGc                                                                                                   >  2:399756/1‑150 (MQ=255)
                                                      tCCGCCATGCGGGTTTCTACACCGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGa                                                                                                 >  2:1298942/1‑150 (MQ=255)
                                                      tCCGCCATGCGGGTTTCTACACCGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGa                                                                                                 >  2:921267/1‑150 (MQ=255)
                                                      gccGCCATGCGTGTTTCTACCCCGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATATCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGa                                                                                                 <  2:310459/149‑1 (MQ=255)
                                                                 ggTTTCTACACCGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGtt                                                                                      >  2:711465/1‑150 (MQ=255)
                                                                           ccGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTg                                                                            >  1:23612/1‑150 (MQ=255)
                                                                            cGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAACTAACGTgg                                                                           <  1:1036049/150‑1 (MQ=255)
                                                                                   tgtgCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTcc                                                                    >  1:160219/1‑150 (MQ=255)
                                                                                   tgtgCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTcc                                                                    <  1:279612/150‑1 (MQ=255)
                                                                                   tgtgCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCATTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTcc                                                                    >  1:54163/1‑150 (MQ=255)
                                                                                    gtgCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCt                                                                   >  2:407167/1‑150 (MQ=255)
                                                                                     tgCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTa                                                                  <  1:671790/150‑1 (MQ=255)
                                                                                      gCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTAc                                                                 >  1:742481/1‑150 (MQ=255)
                                                                                        cATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGt                                                               >  1:1178900/1‑150 (MQ=255)
                                                                                         aTCACTCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGCGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGtt                                                              >  2:671164/1‑150 (MQ=255)
                                                                                                 cgccgGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGcc                                                      >  1:156683/1‑150 (MQ=255)
                                                                                                  gccgGTATCAATACAGATAAACTCACGGCCTTCACTTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGccc                                                     <  2:89047/150‑1 (MQ=255)
                                                                                                   ccgGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCCGCCCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTATATAACGTCGATGATCCTACGTTAGGGCGCCCg                                                    >  1:942162/1‑150 (MQ=255)
                                                                                                     ggTATCAATACAGATAAACTCACGGCCTTCAATTTCCGAACGACCGTACTAACTGTCACGAGTCAGACCCGGCAAATCCGCACCCAGCTCATCTCAGGTGCGCGTGCCACGGATACATAACGTGGAAGATCCTACGGTAGGGCCCCCgcc                                                  >  1:213354/1‑148 (MQ=255)
                                                                                                       tATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTGAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACaa                                                >  2:237951/1‑150 (MQ=255)
                                                                                                        aTCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAg                                               <  1:921267/150‑1 (MQ=255)
                                                                                                        aTCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAg                                               <  1:1298942/150‑1 (MQ=255)
                                                                                                           aaTACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAgcgc                                            >  1:433974/1‑150 (MQ=255)
                                                                                                            aTACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGcgcg                                           <  2:390693/150‑1 (MQ=255)
                                                                                                            aTACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGcgcg                                           >  1:1213958/1‑150 (MQ=255)
                                                                                                             tACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGa                                          <  1:780128/150‑1 (MQ=255)
                                                                                                               cAGATAAACTCACGGCCTTCAATTTCCGCACGTCCGTACTTCCGGTCACGAGTCAGACCAGGGAAATCCGCAAACAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGAcc                                        <  2:1050199/150‑1 (MQ=255)
                                                                                                                aGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACca                                       <  1:298354/150‑1 (MQ=255)
                                                                                                                           ccgCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACAAGCGCATCTCGTGTGCGAGTTAGACGTTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATg                            <  1:91989/148‑1 (MQ=255)
                                                                                                                            ggCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCTGTTGAGAGTTAGACGGTTAAATACCGTGGATTTTCCTACGTTTGGGCGCCCGACAAGCGCGACCCCAGGTACCCTGt                           >  2:207247/1‑150 (MQ=255)
                                                                                                                            ggCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGt                           >  2:368620/1‑150 (MQ=255)
                                                                                                                                ttCAATTTCCGCACGAACGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGCCGGTTAAATAAAGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTaa                       >  1:1311530/1‑150 (MQ=255)
                                                                                                                                 tCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTaaa                      >  2:846675/1‑150 (MQ=255)
                                                                                                                                 tCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAAGCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTaaa                      <  1:60286/150‑1 (MQ=255)
                                                                                                                                   aaTTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACATGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGc                    <  2:909073/150‑1 (MQ=255)
                                                                                                                                    aTTTCCGCACCACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGcc                   <  2:959756/150‑1 (MQ=255)
                                                                                                                                       tCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTaaa                      <  1:1267600/144‑1 (MQ=255)
                                                                                                                                       tCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTaaa                      >  2:1267600/1‑144 (MQ=255)
                                                                                                                                        ccGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCAt               >  2:1000271/1‑150 (MQ=255)
                                                                                                                                         ctcgtgACCGTAAGTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCAtt              >  1:1146981/6‑150 (MQ=255)
                                                                                                                                         cGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCTGGTGCGAGTTAGACGGTTAACTAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCAtt              >  1:243920/1‑150 (MQ=255)
                                                                                                                                          gCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCAttt             <  2:1047132/150‑1 (MQ=255)
                                                                                                                                          gCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCAttt             >  2:1300026/1‑150 (MQ=255)
                                                                                                                                           cgtgACCGTAAGTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGCCGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCAtttt            >  1:916523/4‑150 (MQ=255)
                                                                                                                                           cgtgACCGTAAGTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCAtttt            >  2:208095/4‑150 (MQ=255)
                                                                                                                                           cgtgACCGTAAGTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGGACCATGTTTAACGCCTCAtttt            >  2:205450/4‑150 (MQ=255)
                                                                                                                                           cACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCAtttt            >  2:1050320/1‑150 (MQ=255)
                                                                                                                                             cGACCGTACTTCCGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGTCGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCGCATTTTta          <  1:255702/150‑1 (MQ=255)
                                                                                                                                             cGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTta          >  2:1015669/1‑150 (MQ=255)
                                                                                                                                                ccGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATaa       >  2:365740/1‑150 (MQ=255)
                                                                                                                                                  gTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTTGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAat     <  1:411702/150‑1 (MQ=255)
                                                                                                                                                   tACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAatc    <  2:485200/150‑1 (MQ=255)
                                                                                                                                                     cTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAatcat  <  2:432503/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCATCCACCATAAACAGTACGACGTCCGCTTCTTCAATCGCCAGCAGCGACTGTTCCGCCATGCGGGTTTCTACACCGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCAT  >  NC_000913/2637085‑2637383

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