Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,257,396 G→T 6.0% T139T (ACC→ACA)  ychF ← redox‑responsive ATPase YchF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,257,3960GT 6.0% 221.3 / 3.5 67T139T (ACC→ACA) ychFredox‑responsive ATPase YchF
Reads supporting (aligned to +/- strand):  ref base G (27/36);  new base T (2/2);  total (29/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

GCTTTCTCTTCAGCGCTTAAATCCAGCGCGCGCAGCATACCTGCATTTTCCAACTGGGGCAGGCATTTTTCCAGGACCGCCAGCTCAGCTTTCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTCAATATCGTCAGCCGGGTTAACTTTGCCGGAAACGTGAATGATGTTGTCATTTTCAAAGCAGCGAACAACGTGACCGATCGCTTCGGTTTCACGGATGTTGGTCAGGA  >  NC_000913/1257247‑1257544
                                                                                                                                                     |                                                                                                                                                    
gCTTTCTCTTCAGCGCTTAAATCCAGCGCGCGCAGCATACCTGCATTTTCCAACTGGGGCAGGCATTTTTCCAGGACCGCCAGCTCAGCTTTCGCGTCTTTATCGCCCCCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAg                                                                                                                                                      >  2:372393/1‑150 (MQ=255)
     ctctTCAGCGCTTAAATCCAGCGCGCGCAGCATACCTGCATTTTCCAACT‑GGGCAGGCATTTTTCCAGGACCGCCAGCTCAGCTTTCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCg                                                                                                                                                <  1:1174600/150‑1 (MQ=255)
           aGCGCTTAAATCCAGCGCGCGCAGCATACCTGCATTTTCCAACTGGGGCAGGCATTTTTCCAGGACCGCCAGCTCAGCTTTCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTc                                                                                                                                           >  1:658440/1‑150 (MQ=255)
               cTTAAATCCAGCGCGCGCAGCATACCTGCATTTTCCAACTGGGGCAGGCATTTTTCCAGGACCGCCAGCTCAGCTTTCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGcc                                                                                                                                       >  1:1183940/1‑150 (MQ=255)
                    aTCCGGCGACCGAAGCATAACTGCAATTTCAAACTGGGGCAGGCATTTTTCCAGGACCGCAAGCTCAGATTTAGCGTATTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAgcgc                                                                                                                                  >  2:171758/1‑150 (MQ=255)
                    aTCCAGCGCGCGCAGCATACCTGCATTTTCCAACTGGGGCAGGCATTTTTCCAGGACCGCCAGCTCAGCTTTCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAgcgc                                                                                                                                  >  2:989024/1‑150 (MQ=255)
                    aTCCAGCGCGCGCAGCATACCTGCATTTTCCAACTGGGGCAGGCATTTTTCCAGGACCGCCAGCTCAGCTTTCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAgcgc                                                                                                                                  >  2:784939/1‑150 (MQ=255)
                     tCCAGCGCGCGCAGCATACCTGCATTTTCCAACTGGGGCAGGCATTTTTCCAGGACCGCCAGCTCAGCTTTCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGcc                                                                                                                                 <  1:196784/150‑1 (MQ=255)
                       cAGCGCGCGCAGCATACCTGCATTTTCCAACTGGGGCAGGCATTTTTCCAGGACCGCCAGCTCAGCTTTCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAg                                                                                                                               <  2:41052/150‑1 (MQ=255)
                          cgcgcgcAGCATACCTGCATTTTCCAACTGGGGCAGGCATTTTTCCAGGACCGCCAGCTCAGCTTTCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTc                                                                                                                            <  1:558617/150‑1 (MQ=255)
                              cgcAGCATACCTGCATTTTCCAACTGGGGCAGGCATTTTTCCAGGACCGCCAGCTCAGCTTTCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGgtg                                                                                                                        <  1:1412741/150‑1 (MQ=255)
                                 agcaTACCTGCATTTTCCAACTGGGGCAGGCATTTTTCCAGGACCGCCAGCTCAGCTTTCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTg                                                                                                                     <  1:923214/150‑1 (MQ=255)
                                  gcaTACCTGCATTTTCCAACTGGGGCAGGCATTTTTCCAGGACCGCCAGCTCAGCTTTCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGa                                                                                                                    <  2:236640/150‑1 (MQ=255)
                                       ccTGCATTTTCCAACTGGGGCAGGCATTTTTCCAGGACCGCCAGCTCAGCTTTCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCATGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACt                                                                                                               >  2:558449/1‑150 (MQ=255)
                                       ccTGCATTTTCCAACTGGGGCAGGCATTTTTCCAGGACCGCCAGCTCAGCTTTCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCATGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACt                                                                                                               >  2:553829/1‑150 (MQ=255)
                                       ccTGCATTTTCCAACTGGGGCAGGCATTTTTCCAGGACCGCCAGCTCAGCTTTCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTAGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACt                                                                                                               >  2:1073695/1‑150 (MQ=255)
                                        cTGCATTTTCCAACTGGGGCAGGCATTTTTCCAGGACCGCCAGCTCAGCTTTCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGa                                                                                                                    >  2:559724/1‑144 (MQ=255)
                                        cTGCATTTTCCAACTGGGGCAGGCATTTTTCCAGGACCGCCAGCTCAGCTTTCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGa                                                                                                                    <  1:559724/144‑1 (MQ=255)
                                        cTGCATTTTCCAACTGGGGCAGGCATTTTTCCAGGACCGCCAGCTCAGCTTTCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACtt                                                                                                              >  1:1032091/1‑150 (MQ=255)
                                         tGCATTTTCCAACTGGGGCAGGCATTTTTCCAGGACCGCCAGCTCAGCTTTCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATTAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTc                                                                                                             <  2:546591/150‑1 (MQ=255)
                                         tGCATTTTCCAACTGGGGCAGGCATTTTTCCAGGACCGCCAGCTCAGCTTTCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCATGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTc                                                                                                             <  1:553829/150‑1 (MQ=255)
                                         tGCATTTTCCAACTGGGGCAGGCATTTTTCCAGGACCGCCAGCTCAGCTTTCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCATGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTc                                                                                                             <  1:558449/150‑1 (MQ=255)
                                         agCATTTTCCAACTGGGGCATGCATTTTTCAAGGAACGCCAGCTCAGCTTTAGCGTCTTTATCGACACCATTGGCTTTCTTCTGTACGCGATTAATCGCAAGTTCGCAGGTGTCGAGGTCTGCCAGCTCCAGTTCGGTGTTGATAACtta                                                                                                             >  2:643955/2‑149 (MQ=255)
                                           cATTTTCCAACTGGGGCAGGCATTTTTCCAGGACCGCCAGCTCAGCTTTCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAgg                                                                                                                                                     <  1:710775/108‑1 (MQ=255)
                                           cATTTTCCAACTGGGGCAGGCATTTTTCCAGGACCGCCAGCTCAGCTTTCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAgg                                                                                                                                                     >  2:710775/1‑108 (MQ=255)
                                           cATTTTCCAACTGGGGCAGGCATTTTTCCAGGACCGCCAGCTCAGCTTTCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTCaa                                                                                                           >  1:1281903/1‑150 (MQ=255)
                                            aTTTTCCAACTGGGGCAGGCATTTTTCCAGGACCGCCAGCTCAGCTTTCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTCAat                                                                                                          <  2:1183940/150‑1 (MQ=255)
                                               ttCCAACTGGGGCAGGCATTTTTCCAGGACCGCCAGCTCAGCTTTCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTCAATATc                                                                                                       >  2:745650/1‑150 (MQ=255)
                                                tCCAACTGGGGCAGGCATTTTTCCAGGACCGCCAGCTCAGCTTTCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTCAATATCg                                                                                                      <  1:1064363/150‑1 (MQ=255)
                                                tCCAACTGGGGCAGGCATTTTTCCAGGACCGCCAGCTCAGCTTTCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTCAATATCg                                                                                                      <  1:526413/150‑1 (MQ=255)
                                                    aCTGGGGCAGGCATTTTTCCAGGACCGCCAGCTCAGCTTTCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTCAATATCGTCAg                                                                                                  <  1:989024/150‑1 (MQ=255)
                                                       ggggCAGGCATTTTTCCAGGACCGCCAGCTCAGCTTTCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTCAATATCGTCAGCCg                                                                                               <  2:974274/150‑1 (MQ=255)
                                                       ggggCAGGCATTTTTCCAGGACCGCCAGCTCAGCTTTCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTCAATATCGTCAGCCg                                                                                               >  2:547824/1‑150 (MQ=255)
                                                          gcaggcaTTTTTCCAGGACCGCCAGCTCAGCTTTCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTCAATATCGTCAGCCGGGt                                                                                            <  2:877741/150‑1 (MQ=255)
                                                           caggcaTTTTTCCAGGACCGCCAGCTCAGCTTTCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTCAATATCGTCAGCCGGGtt                                                                                           >  2:295107/1‑150 (MQ=255)
                                                             ggcaTTTTTCCAGGACCGCCAGCTCAGCTTTCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTCAATATCGTCAGCCGGGTTaa                                                                                         >  1:1161464/1‑150 (MQ=255)
                                                              gcaTTTTTCCAGGACCGCCAGCTCAGCTTTCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTCAATATCGTCAGCCGGGTTAAc                                                                                        <  1:275496/150‑1 (MQ=255)
                                                               caTTTTTCCAGGACCGCCAGCTCAGCTTTCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTCAATATCGTCAGCCGGGTTAACt                                                                                       >  1:1328005/1‑150 (MQ=255)
                                                                    ttCCAGGACCGCCAGCTCAGCTTTCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTCAATATCGTCAGCCGGGTTAACTTTGcc                                                                                  <  2:229854/150‑1 (MQ=255)
                                                                     tCCAGGACCGCCAGCTCAGCTTTCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTCAATATCGTCAGCCGGGTTAACTTTGCCg                                                                                 <  2:907961/150‑1 (MQ=255)
                                                                       cAGGACCGCCAGCTCAGCTTTCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTCAATATCGTCAGCCGGGTTAACTTTGCCGGa                                                                               >  1:770249/1‑150 (MQ=255)
                                                                       cAGGACCGCCAGCTCAGCTTTCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTCAATATCGTCAGCCGGGTTAACTTTGCCGGa                                                                               >  1:273700/1‑150 (MQ=255)
                                                                              gCCAGCTCAGCTTTCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTCAATATCGTCAGCCGGGTTAACTTTGCCGGAAACGTGa                                                                        >  1:583421/1‑150 (MQ=255)
                                                                                 agctcagctTTCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTCAATATCGTCAGCCGGGTTAACTTTGCCTGAAACGTGAatg                                                                     <  2:1089805/150‑1 (MQ=255)
                                                                                   ctcagctTTCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTCAATATCGTCAGCCGGGTTAACTTTGCCGGAAACGTGAatgat                                                                   >  2:1174362/1‑150 (MQ=255)
                                                                                     cagctTTCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTCAATATCGTCAGCCGGGTTAACTTTGCCGGAAACGTGAATGAtgt                                                                 >  2:842352/1‑150 (MQ=255)
                                                                                       gctTTCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTCAATATCGTCAGCCGGGTTAACTTTGCCGGAAACGTGAATGAtgttg                                                               <  1:547824/150‑1 (MQ=255)
                                                                                           tCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTCAATATCGTCAGCCGGGTTAACTTTGCCGGAAACGTGAATGATGTTGTCAt                                                           >  1:1251962/1‑150 (MQ=255)
                                                                                             gcgTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTCAATATCGTCAGCCggg                                                                                             <  2:1421459/114‑1 (MQ=255)
                                                                                             gcgTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTCAATATCGTCAGCCggg                                                                                             >  1:1421459/1‑114 (MQ=255)
                                                                                             gcgTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTCAATATCGTCAGCCGGGTTAACTTTGCCGGAAACGTGAATGATGTTGTCAttt                                                         >  2:116313/1‑150 (MQ=255)
                                                                                               gTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTCAATATCGTCAGCCGGGTTAACTTTGCCGGAAACGTGAATGATGTTGTCATTTTc                                                       <  2:1354145/150‑1 (MQ=255)
                                                                                                 cTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTCAATATCGTCAGCCGGGTTAACTTTGCCGGAAACGAGAATGATGTTGTCATTTTCaa                                                     >  2:494341/1‑150 (MQ=255)
                                                                                                        gCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTCAATATCGTCAGCCGGGTTAACTTTGCCGGAAACGTGAATGATGTTGTCATTTTCAAAGCAGCg                                              <  1:295107/150‑1 (MQ=255)
                                                                                                                  gCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTCAATATCGTCAGCCGGGTTAACTTTGCCGGAAACGTGAATGATGTTGTCATTTTCAAAGCAGCGAACAACGTGa                                    <  1:745650/150‑1 (MQ=255)
                                                                                                                   cTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTCAATATCGTCAGCCGGGTTAACTTTGCCGGAAACGTGAATGATGTTGTCATTTTCAAAGCAGCGAACAACGTGAc                                   <  1:784939/150‑1 (MQ=255)
                                                                                                                      tcttctGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTCAATATCGTCAGCCGGGTTAACTTTGCCGGAAACGTGAATGATGTTGTCATTTTCTAAGCAGCGAACAACGTGACCGa                                >  2:1091486/1‑150 (MQ=255)
                                                                                                                        ttctGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTCAATATCGTCAGCCGGGTTAACTTTGCCGGAAACGTGAATGATGTTGTCATTTTCAAAGCAGCGAACAACGTGACCGATc                              <  1:643955/150‑1 (MQ=255)
                                                                                                                        ttctGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTCAATATCGTCAGCCGGGTTAACTTTGCCGGAAACGTGAATGATGTTGTCATTTTCAAAGCAGCGAACAACGTGACCGATc                              <  1:208013/150‑1 (MQ=255)
                                                                                                                               cgcgATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTCAATATCGTCAGCCGGGTTAACTTTGCCGGAAACGTGAATGATGTTGTCATTTTCAAAGCAGCg                                              <  2:374157/127‑1 (MQ=255)
                                                                                                                               cgcgATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTCAATATCGTCAGCCGGGTTAACTTTGCCGGAAACGTGAATGATGTTGTCATTTTCAAAGCAGCg                                              >  1:374157/1‑127 (MQ=255)
                                                                                                                                 cgATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTCAATATCGTCAGCCGGGTTAACTTTGCCGGAAACGTGAATGATGTTGTCATTTTCAAAGCAGCGAACAACGTGACCGATCGCTTCGGtt                     <  2:78526/150‑1 (MQ=255)
                                                                                                                                     gAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTCAATATCGTCAGCCGGGTTAACTTTGCCGGAAACGTGAATGATGTTGTCATATTCAAAGTAGCGAACAACGTGACCGATCGCTTCGGTTTCAc                 >  2:79237/1‑150 (MQ=255)
                                                                                                                                            aCGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTCAATATCGTCAGCCGGGTTAACTTTTCCTGAAAATTTAATTATTTTGTCATTTTCAAATAAGCTAAAAACTTGACCGATCGCTTCGGTTTCACGGATGtt          <  2:658440/150‑1 (MQ=255)
                                                                                                                                            aCGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTCAATATCGTCAGCCGGGTTAACTTTGCCGGAAACGTGAATGATGTTGTCATTTTCAAAGCAGCGAACAACGTGACCGATCGCTTCGGTTTCACGGATGtt          <  1:789188/150‑1 (MQ=255)
                                                                                                                                              gTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTCAATATCGTCAGCCGGGTTAACTTTGCCGGAAACGTGAATGATGTTGTCATTTTCAAAGCAGCGAACAACGTGACCGATCGCTTCGGTTTCACGGATGTTgg        <  2:1251962/150‑1 (MQ=255)
                                                                                                                                                tCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTCAATATCGTCAGCCGGGTTAACTTTGCCGGAAACGTGAATGATGTTGTCATTTTCAAAGCAGCGAACAACGTGACCGATCGCTTCGGTTTCACGGATGTTGGTc      <  2:1281903/150‑1 (MQ=255)
                                                                                                                                                 cGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTCAATATCGTCAGCCGGGTTAACTTTGCCGGAAACGTGAATGATGTTGTCATTTTCAAAGCAGCGAACAACGTGACCGATCGCTTCGGTTTCACGGATGTTGGTCa     <  1:116313/150‑1 (MQ=255)
                                                                                                                                                  gCAGGTGTCGAGGGCTGCCAGCGCCAGTTCGGTGTTGATAACTTCAATATCGTCAGCCGGGTTAACTTTGCCGGAAACGTGAATGATGTTGTCATTTTCAAAGCAGCGAACAACGTGACCGATCGCTTCGGTTTCACGGATGTTGGTCAg    <  1:171758/150‑1 (MQ=255)
                                                                                                                                                    aGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTCAATATCGTCAGCCGGGTTAACTTTGCCGGAAACGTGAATGATGTTGTCATTTTCAAAGCAGCGAACAACGTGACCGATCGCTTCGGTTTCACGGATGTTGGTCAGGa  <  1:1174362/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCTTTCTCTTCAGCGCTTAAATCCAGCGCGCGCAGCATACCTGCATTTTCCAACTGGGGCAGGCATTTTTCCAGGACCGCCAGCTCAGCTTTCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCGATGAATCGCACGTTCGCAGGTGTCGAGGTCTGCCAGCGCCAGTTCGGTGTTGATAACTTCAATATCGTCAGCCGGGTTAACTTTGCCGGAAACGTGAATGATGTTGTCATTTTCAAAGCAGCGAACAACGTGACCGATCGCTTCGGTTTCACGGATGTTGGTCAGGA  >  NC_000913/1257247‑1257544

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