Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,637,235 T→G 7.8% K41T (AAG→ACG) ‡ 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,2350TG 7.8% 207.8 / 8.4 77K41T (AAG→ACG) ‡derGTPase; multicopy suppressor of ftsJ
Reads supporting (aligned to +/- strand):  ref base T (39/32);  new base G (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 = 4.16e-01

CATCCACCATAAACAGTACGACGTCCGCTTCTTCAATCGCCAGCAGCGACTGTTCCGCCATGCGGGTTTCTACACCGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATC  >  NC_000913/2637086‑2637384
                                                                                                                                                     |                                                                                                                                                     
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                                            >  1:1213958/1‑150 (MQ=255)
                                                                                                           aTACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGcgcg                                            <  2:390693/150‑1 (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)
                                                                                                                                                     ttACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAatcatc  <  1:290560/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CATCCACCATAAACAGTACGACGTCCGCTTCTTCAATCGCCAGCAGCGACTGTTCCGCCATGCGGGTTTCTACACCGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATC  >  NC_000913/2637086‑2637384

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: