Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,505,004 G→C 5.2% R339R (CGC→CGG yhfS ← putative aminotransferase YhfS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,505,0040GC 5.2% 326.5 / 5.2 95R339R (CGC→CGGyhfSputative aminotransferase YhfS
Reads supporting (aligned to +/- strand):  ref base G (55/35);  new base C (3/2);  total (58/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

CCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAGCACTCTGGCGGC  >  NC_000913/3504857‑3505153
                                                                                                                                                   |                                                                                                                                                     
ccccTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATgtg                                                                                                                                                     <  2:308260/150‑3 (MQ=255)
ccccTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCg                                                                                                                                                     >  1:1325682/1‑150 (MQ=255)
  ccTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGaa                                                                                                                                                   <  1:708048/150‑1 (MQ=255)
  ccTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGaa                                                                                                                                                   >  1:332101/1‑150 (MQ=255)
  ccTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGaa                                                                                                                                                   >  2:1054973/1‑150 (MQ=255)
  ccTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGaa                                                                                                                                                   >  2:457695/1‑150 (MQ=255)
  ccTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGaa                                                                                                                                                   >  1:291612/1‑150 (MQ=255)
   cTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAAt                                                                                                                                                  <  1:664924/150‑1 (MQ=255)
   cTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAAt                                                                                                                                                  >  2:405203/1‑150 (MQ=255)
      aCCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGc                                                                                                                                               <  1:204624/150‑1 (MQ=255)
       cccTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGca                                                                                                                                              >  1:1160551/1‑150 (MQ=255)
       cccTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGca                                                                                                                                              >  1:1157234/1‑150 (MQ=255)
       cccTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGca                                                                                                                                              >  1:18237/1‑150 (MQ=255)
       cccTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGCAATTTAGCATCACGGTGCTATCCAGCATTAAATACTGGCAATACACTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTCCGATACGGGTTAATGCGAATCGca                                                                                                                                              >  2:1427808/1‑150 (MQ=255)
         cTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGcaca                                                                                                                                            >  2:1210350/1‑150 (MQ=255)
          tAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACaa                                                                                                                                           >  2:722891/1‑150 (MQ=255)
           aaCCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGCATTAACCCGAATCGcagcgg                                                                                                                                          <  2:888293/150‑5 (MQ=255)
            aCCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATg                                                                                                                                         <  1:1361217/150‑1 (MQ=255)
            aCCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGCATTAACCCGAATCGcag                                                                                                                                             >  2:155629/1‑145 (MQ=255)
            aCCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGCATTAACCCGAATCGcag                                                                                                                                             <  1:155629/146‑2 (MQ=255)
             cccTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGt                                                                                                                                        >  2:1167835/1‑150 (MQ=255)
              ccTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGtt                                                                                                                                       >  2:1240366/1‑150 (MQ=255)
                tctcCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCt                                                                                                                                     <  2:452825/150‑1 (MQ=255)
                 ctcCCCAAAGGGGCGAGGGGACCGTTCAGTACTGAATTTAGCTTGACGGTGCTATCGAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTg                                                                                                                                    <  1:754761/150‑1 (MQ=255)
                  tcCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGa                                                                                                                                   >  1:346085/1‑150 (MQ=255)
                    cccAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGAtt                                                                                                                                 >  2:1263183/1‑150 (MQ=255)
                    cccAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGAtt                                                                                                                                 >  2:896583/1‑150 (MQ=255)
                       aaaGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATtgtg                                                                                                                              <  2:1144789/150‑1 (MQ=255)
                       aaaGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATtgtg                                                                                                                              <  1:1004544/150‑1 (MQ=255)
                          ggggCGAGGGTTCCGTTCAGTACTGAATTTAGCTTGACGGTGCTATCCAGGATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCAGCTTCTCTTCACGGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGAt                                                                                                                           <  2:1119337/150‑1 (MQ=255)
                            ggCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTc                                                                                                                         >  2:812036/1‑150 (MQ=255)
                                 ggggACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTg                                                                                                                    >  2:1407820/1‑150 (MQ=255)
                                       cGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGaaa                                                                                                              <  1:660699/150‑1 (MQ=255)
                                        gTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAAc                                                                                                             <  2:1002600/150‑1 (MQ=255)
                                         ttCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACg                                                                                                            >  1:1390037/1‑150 (MQ=255)
                                           cAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGtt                                                                                                          >  1:1244628/1‑150 (MQ=255)
                                             gTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTcc                                                                                                        <  1:127538/150‑1 (MQ=255)
                                                  gAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGaaa                                                                                                   >  1:636078/1‑150 (MQ=255)
                                                   aaTTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAg                                                                                                  >  1:1399169/1‑150 (MQ=255)
                                                    aTTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCTATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGaaatg                                                                                                 <  2:1483478/150‑3 (MQ=255)
                                                     tttAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGc                                                                                                <  2:260835/150‑1 (MQ=255)
                                                        aGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCTCCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCGTGGCGAAACGTTCCGGAAAGGCGAt                                                                                             <  1:1427808/150‑1 (MQ=255)
                                                        aGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGAt                                                                                             <  2:346085/150‑1 (MQ=255)
                                                         gCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGTTTAATGCTAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATa                                                                                            <  2:726495/150‑1 (MQ=255)
                                                         gCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATa                                                                                            >  2:603160/1‑150 (MQ=255)
                                                         gCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATa                                                                                            <  1:1240366/150‑1 (MQ=255)
                                                          cTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGCTTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATaa                                                                                           >  2:42443/1‑150 (MQ=255)
                                                            tGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATaaaa                                                                                         >  2:698888/1‑150 (MQ=255)
                                                                 gTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAgcg                                                                                    >  1:656741/1‑150 (MQ=255)
                                                                   gCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAgcggc                                                                                  >  2:390143/1‑150 (MQ=255)
                                                                    cTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAgcggcg                                                                                 <  2:141600/150‑1 (MQ=255)
                                                                      aTCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCggg                                                                               >  1:668972/1‑150 (MQ=255)
                                                                       tCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGa                                                                              >  2:22504/1‑150 (MQ=255)
                                                                        ccAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGAt                                                                             >  2:596059/1‑150 (MQ=255)
                                                                        ccAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGAt                                                                             >  2:1136396/1‑150 (MQ=255)
                                                                                aaaTACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATAtt                                                                     >  2:988577/1‑150 (MQ=255)
                                                                                  aTACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTc                                                                   >  1:1054616/1‑150 (MQ=255)
                                                                                  aTACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTc                                                                   <  2:744736/150‑1 (MQ=255)
                                                                                       ggCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGActc                                                              >  1:1293832/1‑150 (MQ=255)
                                                                                       ggCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGActc                                                              >  2:77233/1‑150 (MQ=255)
                                                                                        gCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGTGGGATTTCATATTTCGActct                                                             <  1:603160/150‑1 (MQ=255)
                                                                                         cAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGGGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTg                                                            <  1:1196390/150‑1 (MQ=255)
                                                                                          aaTACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGc                                                           <  1:22504/150‑1 (MQ=255)
                                                                                             aCTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGcacc                                                        <  1:1167835/150‑1 (MQ=255)
                                                                                               tctcACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGcaccca                                                      >  1:1189509/1‑150 (MQ=255)
                                                                                                ctcACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGcacccac                                                     <  2:291612/150‑1 (MQ=255)
                                                                                                  cACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATTTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATTTTTCGACTCTGCACCCACCg                                                   >  2:928701/1‑150 (MQ=255)
                                                                                                    cGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCggg                                                 <  1:1265821/150‑1 (MQ=255)
                                                                                                     gCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGcac                                                         >  2:20998/1‑141 (MQ=255)
                                                                                                     gCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGcac                                                         <  1:20998/141‑1 (MQ=255)
                                                                                                      cAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTa                                               >  2:435084/1‑150 (MQ=255)
                                                                                                      cAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTa                                               >  1:339074/1‑150 (MQ=255)
                                                                                                       aaaaTCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTaa                                              >  2:262886/1‑150 (MQ=255)
                                                                                                       aaaaTCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTaa                                              <  2:656741/150‑1 (MQ=255)
                                                                                                        aaaTCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTaag                                             <  1:457695/150‑1 (MQ=255)
                                                                                                        aaaTCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTaag                                             <  1:1054973/150‑1 (MQ=255)
                                                                                                          aTCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTaaggc                                           <  1:469021/150‑1 (MQ=255)
                                                                                                               cagcaCCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGAGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTaaggcaaggc                                      <  1:405203/150‑1 (MQ=255)
                                                                                                                 gcgcggTCTGTTCACCGCTGCGATTCGGGTTAATGCGAATCGCGCATTGTTGTGATTGTGGATTCGCGTGGCGAAACGTTGCGGAAAGGCGATAAAAGAGCGGCGGGGTTTCATATTTCGGCTCTGCACCGGCCGGGTAAGGCGCGGCAc                                    >  2:416742/6‑150 (MQ=255)
                                                                                                                 gcaCCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCAc                                    >  1:592587/1‑150 (MQ=255)
                                                                                                                    ccGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACcgc                                 >  2:426218/1‑150 (MQ=255)
                                                                                                                        ctctTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCtt                             >  2:271135/1‑150 (MQ=255)
                                                                                                                             cACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTgcg                        <  1:271135/150‑1 (MQ=255)
                                                                                                                              cccGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGCCTCTGCACCCACCGGGTACCGCAAGGCACCGCGCTTTTgcgc                       >  1:1469820/2‑150 (MQ=255)
                                                                                                                                ttgtGCGATTCGCATTAACCCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCct                     >  1:567596/4‑150 (MQ=255)
                                                                                                                                cGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTAGACTCTGCACCAACAGGGTAAGGCAAGGCAACGCGATTTTgcgaca                     >  2:720079/1‑147 (MQ=255)
                                                                                                                                 tgtGCGATTCGCATTAACCCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCctc                    >  2:493226/3‑150 (MQ=255)
                                                                                                                                 gCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCctc                    >  1:701954/1‑150 (MQ=255)
                                                                                                                                      gATTCGGGATAATGCGAATCGCACAAAGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTCTACCCACCGGGTAAGGCAAGGCACCGCGCTATTGCGCCTCTTCCa               >  2:1316781/1‑150 (MQ=255)
                                                                                                                                          cGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCATGGCACCGCGCTTTTGCGCCTCTTCCAGCAc           <  2:1399169/150‑1 (MQ=255)
                                                                                                                                           gggTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAGCAct          >  2:100981/1‑150 (MQ=255)
                                                                                                                                            ggTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAGCActc         >  1:1080114/1‑150 (MQ=255)
                                                                                                                                                aaTGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAGCACTCTggc     >  1:553201/1‑150 (MQ=255)
                                                                                                                                                 aTGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAGCACTCTggcg    >  2:712698/1‑150 (MQ=255)
                                                                                                                                                 aTGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAGCACTCTggcg    <  2:1157234/150‑1 (MQ=255)
                                                                                                                                                  tGCGAATCGCACAATGTTCTGATTGTGGATTCGACTGGAGAAACGTTCCGTAAAGGCGATAAAACAGCGGCGGAATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCATTTGCGCCTCTTCCACCCCTCAcgcgg   >  2:235499/1‑150 (MQ=255)
                                                                                                                                                   gCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAGCACTCTggcggc  <  2:1293832/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAGCACTCTGGCGGC  >  NC_000913/3504857‑3505153

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: