Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 263,909 C→T R33R (CGG→CGA yafW ← antitoxin of the YkfI‑YafW toxin‑antitoxin pair

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913263,9090CT97.4% 250.6 / ‑2.0 78R33R (CGG→CGAyafWantitoxin of the YkfI‑YafW toxin‑antitoxin pair
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (35/41);  total (35/43)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.89e-01

TACAGGGTAACGCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCGGTCAGCCAGATAATGCAGTCGGTTGCCTTCCTGCACCAGACGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTG  >  NC_000913/263772‑264038
                                                                                                                                         |                                                                                                                                 
tACAGGGTAACGCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCg                                                                                                                                 <  2:991897/140‑1 (MQ=255)
tACAGGGTAACGCAGTGGGCATGGTGCGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCg                                                                                                                                 >  1:922402/1‑140 (MQ=255)
   aGGGTAACGCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTc                                                                                                                              <  1:751345/140‑1 (MQ=255)
     ggTAACGCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAg                                                                                                                            <  2:954075/140‑1 (MQ=255)
     ggTAACGCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAg                                                                                                                            >  1:792682/1‑140 (MQ=255)
     ggTAACGCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAg                                                                                                                            >  2:171748/1‑140 (MQ=255)
     ggTAACGCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAg                                                                                                                            <  2:94577/140‑1 (MQ=255)
     ggTAACGCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAg                                                                                                                            >  2:171956/1‑140 (MQ=255)
     ggTAACGCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACCGAAGTTGCCGGTGATGCTGGCTCGGTCAg                                                                                                                            >  1:396816/1‑140 (MQ=255)
     ggTAACGCAGTGGGCATGGTGGGGGCTGAGTCCACTGGTCGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAg                                                                                                                            <  2:391764/140‑1 (MQ=255)
       tAACGCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCGGTCAGcc                                                                                                                          <  1:198056/140‑1 (MQ=255)
       tAACGCAGTGGGCATGGTGGGGGCTGAGTTCACCAGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGAAGTCACTGAACTAGCCGGTGATCCTGGCTCGGTCAGcc                                                                                                                          >  1:367935/1‑140 (MQ=37)
         aCGCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAg                                                                                                                        >  1:632628/1‑140 (MQ=255)
          cGCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGa                                                                                                                       >  1:981523/1‑140 (MQ=255)
          cGCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGa                                                                                                                       >  1:204843/1‑140 (MQ=255)
             aGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCAGTCGATGTCACTGCACTGGCCGGTGACGCTGACTCGGTCAGCCAGATaa                                                                                                                    >  2:277295/1‑140 (MQ=255)
               tGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCGTCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATg                                                                                                                  >  2:23982/1‑140 (MQ=255)
               tGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATg                                                                                                                  >  1:65074/1‑140 (MQ=255)
               tGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATg                                                                                                                  >  1:853737/1‑140 (MQ=255)
               tGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATg                                                                                                                  >  2:254870/1‑140 (MQ=255)
               tGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATg                                                                                                                  >  2:973433/1‑140 (MQ=255)
               tGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATg                                                                                                                  >  1:37801/1‑140 (MQ=255)
               tGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATg                                                                                                                  >  2:596940/1‑140 (MQ=255)
                gggCATGGTGGGGGCTGAGTTCTCCTGTGATCAGCATGGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCGGTCAGCCAGATAATGc                                                                                                                 <  1:340244/140‑1 (MQ=255)
                  gCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAg                                                                                                               <  2:155460/140‑1 (MQ=255)
                   cATGGTGGGGGCTGACTTCACCAGTTGTCTGCATCGACTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAGt                                                                                                              <  1:277295/140‑1 (MQ=255)
                           gggCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGACTCGGTCACCCAGTTAATGCAGTCGGTTGcc                                                                                                      >  2:406538/1‑140 (MQ=255)
                           gggCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATCCTGACTCGGTCGGCCAGATAATGCAGTCGGTTGcc                                                                                                      >  2:423977/1‑140 (MQ=255)
                              cTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAGTCGGTTGCCTTc                                                                                                   <  2:922402/140‑1 (MQ=255)
                              cTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAGTCGGTTGCCTTc                                                                                                   <  2:65074/140‑1 (MQ=255)
                                 aGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAGTCGGTTGCCTTCCTg                                                                                                >  1:696893/1‑140 (MQ=255)
                                 aGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAGTCGGTTGCCTTCCTg                                                                                                >  2:132424/1‑140 (MQ=255)
                                 aGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAGTCGGTTGCCTTCCTg                                                                                                >  1:885155/1‑140 (MQ=255)
                                 aGTTCACCGGGGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAGTCGGTTGCCTTCCTg                                                                                                >  2:52009/1‑140 (MQ=255)
                                       ccGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAGTCGGTTGCCTTCCTGCACCAg                                                                                          <  2:13134/140‑1 (MQ=255)
                                          gtggtCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAGTCGGTTGCCTTCCTGCACCAGACg                                                                                       <  1:107027/140‑1 (MQ=255)
                                             gtCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAGTCGGTTGCCTTCCTGCACCAGACGGGc                                                                                    >  1:810256/1‑140 (MQ=255)
                                                 gCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAGTCGGTTGCCTTCCTGCACCAGACGGGCTccc                                                                                <  2:696893/140‑1 (MQ=255)
                                                    tCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAGTCGGTTGCCTTCCTGCACCAGACGGGCTCCCAgg                                                                             <  2:115877/140‑1 (MQ=255)
                                                    tCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAGTCGGTTGCCTTCCTGCACCAGACGGGCTCCCAgg                                                                             <  1:792643/140‑1 (MQ=255)
                                                      gATTCCATCTGGAGGCTAAAGTACGAGAATGTTTCATTCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAGTCGGTTGCCTTCCTGCACCAGACGGGCTCCCAGgcg                                                                           <  1:423977/140‑1 (MQ=255)
                                                      gATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAGTCGGTTGCCTTCCTGCACCAGACGGGCTCCCAGgcg                                                                           <  2:720677/140‑1 (MQ=255)
                                                      gATTCCATCTGGAGGATAAAGTGCGGGAAGGTTTCATACAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAGTCGGTTGCCTTCCTGCACCAGACGGGCTCCCAGgcg                                                                           <  1:406538/140‑1 (MQ=255)
                                                       aTTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAGTCGGTTGCCTTCCTGCACCAGACGGGCTCCCAGgcgc                                                                          >  2:1007463/1‑140 (MQ=255)
                                                            aTCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAGTCGGGTGCCTTCCTGCACCAGACGGGGCCACAGGCGCAgcgt                                                                     >  2:725927/1‑140 (MQ=255)
                                                                  aggaTAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAGTCGGTTGCCTTCCTGCACCAGACGGGCTCCCAgg                                                                             >  2:779162/1‑126 (MQ=255)
                                                                  aggaTAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAGTCGGTTGCCTTCCTGCACCAGACGGGCTCCCAgg                                                                             <  1:779162/126‑1 (MQ=255)
                                                                  aggaTAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAGTCGGTTGCCTTCCTGCACCAGACGGGCTCCCAgg                                                                             <  1:733973/126‑1 (MQ=255)
                                                                  aggaTAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAGTCGGTTGCCTTCCTGCACCAGACGGGCTCCCAgg                                                                             >  2:733973/1‑126 (MQ=255)
                                                                  aggaTAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAGTCGGTTGCCTTCCTGCACAAGAAGGACTCCCAGGCGCAGTGTAAtctc                                                               >  1:454063/1‑140 (MQ=255)
                                                                   ggaTAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAGTCGGTTGCCTTCCTGCACCAGACGGGCTCCCAGGCGCAGTGTAATCTcc                                                              <  2:204843/140‑1 (MQ=255)
                                                                         aGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAGTCGGTTGCCTTCCTGCACCAGACGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGc                                                        <  1:596940/140‑1 (MQ=255)
                                                                          gTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAGTCGGTTGCCTTCCTGCACCAGACGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCa                                                       <  2:810256/140‑1 (MQ=255)
                                                                          gTGCGGGAATGTTTCATACAGCATCCGGCATACGATGTCACTGAACTTTCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAGTCGGTTGCCTTCCTGCACCAGACGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCa                                                       <  2:367935/140‑1 (MQ=255)
                                                                           tGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAGTCGGTTGCCTTCCTGCACCAGACGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAg                                                      >  2:765527/1‑140 (MQ=255)
                                                                             cGGGAATGTTTCATCCCTCTTCGGGCATTCGGTTTTACTTAACTTGCCAGTGATGCTGGCTCGGTCAGCCAGATAATGCAGTCGGTTGCCTTCCTGCACCAGACGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGc                                                    <  2:495713/140‑1 (MQ=255)
                                                                              gggAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAGTCGGTTGCCTTCCTGCACCAGACGGGTTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGcc                                                   >  1:151419/1‑140 (MQ=255)
                                                                              gggAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAGTCGGTTGCCTTCCTGCACCAGACGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGcc                                                   <  1:944766/140‑1 (MQ=255)
                                                                                   tGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAGTCGGTTGCCTTCCTGCACCAGACGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTgg                                              <  2:792682/140‑1 (MQ=255)
                                                                                   tGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAGTCGGTTGCCTTCCTGCACCAGACGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTgg                                              <  1:765527/140‑1 (MQ=255)
                                                                                   tGTTTCATCCAGCTTCCGGCATCCGATGTAACTGAACATGCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAGTCGGTTGCCTTCCTGCACCAGACGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTgg                                              <  1:372108/140‑1 (MQ=255)
                                                                                     tttCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCTGTCAGCCAGATAATGCAGTCGGTTGCCTTCCTGCACCAGACGGGCTCCCAGGCGCAGTGTAATCTCCCCCTGCAGGCCCCTGGTa                                            >  2:552429/1‑140 (MQ=255)
                                                                                     tttCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAGTCGGTTGCCTTCCTGCACCAGACGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTa                                            >  2:242443/1‑140 (MQ=255)
                                                                                        cATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAGTCGGTTGCCTTCCTGCACCAGACGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAggg                                         <  2:885155/140‑1 (MQ=255)
                                                                                        cATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAGTCGGTTGCCTTCCTGCACCAGACGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAggg                                         <  1:254870/140‑1 (MQ=255)
                                                                                           ccAGCTTCCGGCATTCGGTGTCGCTGAACTTTCCGCTGAAGCTGGCTCGGTCAGCCAGATAATGCAGTCGGTTGCCTTCCTGCACCAGACGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTg                                      <  2:396816/140‑1 (MQ=255)
                                                                                           ccAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAGTCTGTTGCCTTCCTGCACCAGACGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTg                                      <  1:1000670/140‑1 (MQ=255)
                                                                                           ccAGCTTCCGGCATTAGATGTCGCTCAACATGCCGGTCATGCTGGCTCGGTCAGCCAGATAATGCAGTCGGTTGCCTTCCTGCACCAGACGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTg                                      <  1:469014/140‑1 (MQ=255)
                                                                                                   cGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAGTCGGTTGCCTTCCTGCACCAGACGGGCTCCCAGGCGCAGTGGAATTTCCGGCTGCAGGGCCCTCGTAGGGTTGCTCATTAc                              >  2:341139/1‑140 (MQ=255)
                                                                                                     gCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAGTCGGTTGCCTTCCTGCACCAGACGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAg                            <  1:171956/140‑1 (MQ=255)
                                                                                                     gCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAGTCGGATGCCTTCCTGCACCAGACGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAg                            <  2:853737/140‑1 (MQ=255)
                                                                                                       aTTCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAGTCGGTTGCCTTCCTGCACCAGACGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGtt                          <  1:973433/140‑1 (MQ=255)
                                                                                                       aTTCGATGTCACTGAACTTGCCCGTTATGCTGGCTCGATCAGCCAGATAATGCAGTCGGTTGCCTTCCTGCACCAGACGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGtt                          <  1:326013/140‑1 (MQ=255)
                                                                                                         tCGATGTCACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAGTCGGTTGCCTTCCTGCACCAGACGGGTTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTtct                        <  2:151419/140‑1 (MQ=255)
                                                                                                                cACTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAGTCGGTTGCCTTCCTGCACCAGACGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTAt                 >  1:725394/1‑140 (MQ=255)
                                                                                                                 aCTGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAGTCGGTTGCCTTCCTGCACCAGACGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTAtt                >  2:36092/1‑140 (MQ=255)
                                                                                                                   tGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAGTCGGTTGCCTTCCTGCACCAGACGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACa                             >  1:230961/1‑125 (MQ=255)
                                                                                                                   tGAACTTGCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAGTCGGTTGCCTTCCTGCACCAGACGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACa                             <  2:230961/125‑1 (MQ=255)
                                                                                                                        ttGCCGGTGATGCTGGCTCGGTCAGCCAGATAACGAAGTCGGTTGCCTGCCTGCACCAGACGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGtt         <  2:454063/140‑1 (MQ=255)
                                                                                                                         tGCCGGTGATGCTGGCTCGGTCAGCCAGATAATGCAGTCGGTTGCCTTCCTGCACCAGACGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGTTCAGTAGAGTTCTGCACTATTGTCAGTTc        >  1:424014/1‑140 (MQ=255)
                                                                                                                               tGATGCTGGCTCGGTCACCCAGTTAATGCAGTCGGCTGCCTTCCTGCATCAGACGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTg  <  2:424014/140‑1 (MQ=38)
                                                                                                                                         |                                                                                                                                 
TACAGGGTAACGCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCGGTCAGCCAGATAATGCAGTCGGTTGCCTTCCTGCACCAGACGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTG  >  NC_000913/263772‑264038

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