Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,376,952 T→A 5.2% C41S (TGT→AGT) ‡ ycjS → putative NADH‑binding oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,376,9520TA 5.2% 278.2 / 4.1 96C41S (TGT→AGT) ‡ycjSputative NADH‑binding oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (50/41);  new base A (3/2);  total (53/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.89e-01

TGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTACGAACATACCCTG  >  NC_000913/1376805‑1377101
                                                                                                                                                   |                                                                                                                                                     
tGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGctgtctgt                                                                                                                                                     <  1:1281513/150‑1 (MQ=255)
  gTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGa                                                                                                                                                   <  1:1368029/150‑1 (MQ=255)
     gCGTAACTGCAAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATACTAGGCGCAGGCCAGGTGGTGGATAAGGTGCATGCTACGAACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAg                                                                                                                                                <  1:253272/150‑1 (MQ=255)
       gTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCCGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTATGACAgcc                                                                                                                                              <  1:840182/150‑1 (MQ=255)
        tACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAgccg                                                                                                                                             <  1:1186473/150‑1 (MQ=255)
         aCCTGCTAAGAGGTCTTTGTGAAATGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAgccgc                                                                                                                                            >  2:930075/1‑150 (MQ=255)
         aCCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGGCACAGACAgccac                                                                                                                                            <  2:1219770/150‑3 (MQ=255)
          ccTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAgccgcc                                                                                                                                           <  1:1052264/150‑1 (MQ=255)
          ccTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGGCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCACAGACAgccacc                                                                                                                                           <  1:1475711/150‑1 (MQ=255)
           cTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCt                                                                                                                                          <  1:792291/150‑1 (MQ=255)
           cTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCt                                                                                                                                          <  1:737409/150‑1 (MQ=255)
            tGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCtt                                                                                                                                         >  1:1106650/1‑150 (MQ=255)
            tGCTAAGAGGTCTTGGTGAAAAGGGCAATGACAAGCTCTCCCCTGCGGGGCGCGATAATAGGCGCAGGCCAGGCGGCGGATAAGGGTCATGCTTCGTACTACTGCACCCCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCttt                                                                                                                                        <  1:872219/150‑1 (MQ=255)
               tAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGTTAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTcc                                                                                                                                      <  2:943294/150‑1 (MQ=255)
               tAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCAATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTcc                                                                                                                                      >  1:312527/1‑150 (MQ=255)
                 agagGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGGTGCGGGTCGCGATAATAGGCGCAGTCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCca                                                                                                                                    <  2:686552/150‑1 (MQ=255)
                  gcggTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCGCCCGCAACGCTCTGGAACTGGTGGCTGGCTGTGCCCGCCGCCTTTCCcag                                                                                                                                   >  1:1145394/3‑150 (MQ=255)
                   agGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCcagg                                                                                                                                  <  1:605345/150‑1 (MQ=255)
                    ggTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCcaggc                                                                                                                                 >  2:322557/1‑150 (MQ=255)
                     gTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCcaggcg                                                                                                                                >  1:1376001/1‑150 (MQ=255)
                     gTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCcaggcg                                                                                                                                >  1:1382651/1‑150 (MQ=255)
                        tttGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCcaggcgcag                                                                                                                             >  1:982084/1‑150 (MQ=255)
                        tttGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCcaggcgcag                                                                                                                             >  2:514816/1‑150 (MQ=255)
                          tgtgAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCcaggcgcaggc                                                                                                                           >  1:104616/1‑150 (MQ=255)
                          tgtgAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCcaggcgcaggc                                                                                                                           <  2:483673/150‑1 (MQ=255)
                          tgtgAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCcaggcgcaggc                                                                                                                           >  1:865179/1‑150 (MQ=255)
                           gtgAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCcaggcgcaggcg                                                                                                                          <  2:1609022/150‑1 (MQ=255)
                           gtgAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCcaggcacaggcg                                                                                                                          >  2:1136033/1‑150 (MQ=255)
                            tgAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCaggcgcaggcgc                                                                                                                         >  1:551877/1‑150 (MQ=255)
                               aaaGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTgg                                                                                                                      <  2:369982/150‑1 (MQ=255)
                                  gTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGGGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAg                                                                                                                   <  2:1299781/150‑1 (MQ=255)
                                    gCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGaa                                                                                                                 >  2:1123748/1‑150 (MQ=255)
                                     cAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGaaa                                                                                                                >  2:878889/1‑150 (MQ=255)
                                     cAATGACAAGCTCTCCGCTGCGGGTCGCGATAAAAGGCGCAGGCCAGCTGGCGCATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTACCCAGGCGCACGCGCTGGCAGaaa                                                                                                                >  1:1118340/1‑150 (MQ=255)
                                       aTGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGaaaaa                                                                                                              >  1:1050423/1‑150 (MQ=255)
                                            aaGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACgg                                                                                                         <  1:659669/150‑1 (MQ=255)
                                              gCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGa                                                                                                       >  1:316927/1‑150 (MQ=255)
                                               ctctCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGaa                                                                                                      >  1:157296/1‑150 (MQ=255)
                                               ctctCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGaa                                                                                                      >  2:321907/1‑150 (MQ=255)
                                               ctctCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGaa                                                                                                      >  2:591251/1‑150 (MQ=255)
                                                   ccGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCa                                                                                                  <  1:1392480/150‑1 (MQ=255)
                                                        gCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCgt                                                                                             <  2:865179/150‑1 (MQ=255)
                                                         cGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCgtg                                                                                            >  2:1141889/1‑150 (MQ=255)
                                                           ggTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGtgtg                                                                                          <  1:347558/150‑1 (MQ=255)
                                                           ggTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGtgtg                                                                                          <  1:1297311/150‑1 (MQ=255)
                                                            gTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTgg                                                                                         >  1:1371169/1‑150 (MQ=255)
                                                              cgcgATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCCTTCCCAGGCGCAGGCGCTGGCAAAAAAATACGGGAATGCATCCGTGTGGGa                                                                                       >  1:1203853/1‑150 (MQ=255)
                                                                cgATAATAGGCGCAGGCCAGGGGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACg                                                                                     <  1:1703419/150‑1 (MQ=255)
                                                                     ataGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGCTCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGGATCCGTGTGGGACGATCCg                                                                                >  1:284825/1‑150 (MQ=255)
                                                                         gcgcAGGCCAGGGGGCGGATAAGGTTCATGCTTCGTACGACTGCACCCTCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGGGGGACGATCCGCAgg                                                                            <  1:250972/150‑1 (MQ=255)
                                                                          cgcAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGc                                                                           <  2:1382651/150‑1 (MQ=255)
                                                                          cgcAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGc                                                                           <  2:1376001/150‑1 (MQ=255)
                                                                           gcAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGcc                                                                          <  2:426640/150‑1 (MQ=255)
                                                                             aGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCAt                                                                        >  2:776497/1‑150 (MQ=255)
                                                                              ggCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCAtg                                                                       >  2:80513/1‑150 (MQ=255)
                                                                              ggCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCAtg                                                                       >  1:1297306/1‑150 (MQ=255)
                                                                              ggCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCAtg                                                                       >  2:1503742/1‑150 (MQ=255)
                                                                              ggCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCAtg                                                                       >  1:729716/1‑150 (MQ=255)
                                                                              ggCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCAtg                                                                       >  2:731416/1‑150 (MQ=255)
                                                                              ggCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCAtg                                                                       >  1:1600477/1‑150 (MQ=255)
                                                                               gCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCAtgc                                                                      <  1:1123748/150‑1 (MQ=255)
                                                                                ccAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATg                                                                                                    >  1:1039046/1‑119 (MQ=255)
                                                                                ccAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATg                                                                                                    <  2:1039046/119‑1 (MQ=255)
                                                                                  aGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCAtgctgc                                                                   <  1:1685498/150‑1 (MQ=255)
                                                                                      ggCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTTGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGTCAGAAAAATACGGGAATGCATCCGTTTGGGACGATCCGCAGTCAATGCTGCTGGc                                                               >  2:309956/1‑150 (MQ=255)
                                                                                       gCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCg                                                              >  1:141849/1‑150 (MQ=255)
                                                                                        cGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCgg                                                             <  1:1503742/150‑1 (MQ=255)
                                                                                          gagaaGGTGCCTGCTTCGTACTACTGCCCCCGCAACGAGCTGGAACTGGGGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCGGGCAGAAAAAGACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTg                                                           <  1:640266/147‑1 (MQ=255)
                                                                                              aGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAAc                                                       <  1:1141889/150‑1 (MQ=255)
                                                                                               ggTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTTTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAAcc                                                      <  1:930075/150‑1 (MQ=255)
                                                                                                gTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGGGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCt                                                     <  1:878889/150‑1 (MQ=255)
                                                                                                       cTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTgg                                              >  2:653890/1‑150 (MQ=255)
                                                                                                       cTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTgg                                              >  1:1208049/1‑150 (MQ=255)
                                                                                                           gTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAg                                          <  2:1106650/150‑1 (MQ=255)
                                                                                                            tactacTGCCCCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGc                                         >  2:837168/1‑150 (MQ=255)
                                                                                                               tactGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTc                                      <  2:157296/150‑1 (MQ=255)
                                                                                                                 ctGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTg                                    <  2:920917/150‑1 (MQ=255)
                                                                                                                 ctGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTg                                    >  1:1568710/1‑150 (MQ=255)
                                                                                                                   gCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCt                                  >  1:614438/1‑150 (MQ=255)
                                                                                                                      cccGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCAc                               >  2:7300/1‑150 (MQ=255)
                                                                                                                      cccGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCAc                               >  1:1009619/1‑150 (MQ=255)
                                                                                                                                tCTGGAACTGGTGGCTGTCTGTGACAGCCGCCCTTCCCAGGCGCAGGCGCTGGCAAAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGttt                     <  2:1203853/150‑1 (MQ=255)
                                                                                                                                 cTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGtttt                    >  2:1317247/1‑150 (MQ=255)
                                                                                                                                  tGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTc                   <  1:776497/150‑1 (MQ=255)
                                                                                                                                   ggAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCATAAAAATACGGGAATGCATCCGTGTGGGACGATACGCAGGCCATGCTGATGGCGGTGAAACAGGATGTGGTTAGCGTCTGCTCACCTAACCGTTATCa                  >  2:223343/1‑150 (MQ=255)
                                                                                                                                   ggAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGTGAATGCATCAGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAAACTGATGTGGTTAGCGTCTGCTCAACTAAACGTTTTCa                  >  2:223340/1‑150 (MQ=255)
                                                                                                                                      aCTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCAc                               >  2:336902/1‑134 (MQ=255)
                                                                                                                                      aCTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCAc                               <  1:336902/134‑1 (MQ=255)
                                                                                                                                       aaggcggCTGTCACAGACAGCCACCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTAc              >  2:1679355/6‑150 (MQ=255)
                                                                                                                                        tggtggCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTACg             <  1:321907/150‑1 (MQ=255)
                                                                                                                                        tggtggCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTACg             <  1:591251/150‑1 (MQ=255)
                                                                                                                                        tggtggCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTACg             >  1:1379957/1‑150 (MQ=255)
                                                                                                                                          gtggCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTACGaa           >  2:749615/1‑150 (MQ=255)
                                                                                                                                          gcggCTGTCACAGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCCTCCGTGTGGGACGATCCGCAGGCCATGCTGCGGGCGGTGAACCCTGATGTGGTTCGCGTCGGCTCACCTAACCGTGTTCATGACgcc           >  1:835994/3‑148 (MQ=255)
                                                                                                                                          gcggCTGTCACAGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTACGaa           >  1:1219726/3‑150 (MQ=255)
                                                                                                                                            ggCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTACGAACa         >  2:1045279/1‑150 (MQ=255)
                                                                                                                                             gCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTACGAACAt        >  1:752527/1‑150 (MQ=255)
                                                                                                                                                  ctgtGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTACGAACATACCCt   >  1:1556370/1‑150 (MQ=255)
                                                                                                                                                   tgtgACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTACGAACATACCCTg  >  2:1642688/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTACGAACATACCCTG  >  NC_000913/1376805‑1377101

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