Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,779,756 T→A 5.0% S991C (AGC→TGC)  ypjA ← adhesin‑like autotransporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,779,7560TA 5.0% 257.4 / 2.8 97S991C (AGC→TGC) ypjAadhesin‑like autotransporter
Reads supporting (aligned to +/- strand):  ref base T (42/47);  new base A (2/3);  total (45/52)
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.68e-01

CCAGACACTGCCCCGGCTGGCAGCCGTCAGGTCGGCAAGCGTGGTATCGACGCCATTTCCCAGAGTCAATGTCGCGCTATCGGTAATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGGAGATCTTCTGCCCGGGCCAGAGCCTGAAACTCATCTTTTGACCGCCCAAGGGTATATTGCCCGCCAGAGTTAA  >  NC_000913/2779607‑2779902
                                                                                                                                                     |                                                                                                                                                  
ccAGACACTGCCCCGGCTGGCAGCCGTCAGGTCGGCAAGCGTGGTATCGACGCCATTTCCCAGAGTCAATGTCGCGCTATCGGTAATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCt                                                                                                                                                    >  1:276629/1‑150 (MQ=255)
   gACACTGCCCCGGCTGGCAGCCGTCAGGTCGGCAAGCGTGGTATCGACGCCATTTCCCAGAGTCAATGTCGCGCTATCGGTAATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTcc                                                                                                                                                 <  2:1635554/150‑1 (MQ=255)
     cacTGCCCCGGCTGGCAGCCGTCAGGTCGGCAAGCGTGGTATCGACGCCATTTCCCAGAGTCAATGTCGCGCTATCGGTAATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTg                                                                                                                                               >  1:1330250/1‑150 (MQ=255)
     cacTGCCCCGGCTGGCAGCCGTCAGGTCGGCAAGCGTGGTATCGACGCCATTTCCCAGAGTCAATGTCGCGCTATCGGTAATCCGGACCGCCCCGTCGTGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTg                                                                                                                                               >  1:2033002/1‑150 (MQ=255)
          ccccGGCTGGCGGCCGTCAGGTCGGCAAGCGTGGTATCGACGCCATTTCCCAGAGTCAATGTCGCGCTATCGGTAATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTcgcg                                                                                                                                          >  1:1413438/1‑150 (MQ=255)
          ccccGGCTGGCAGCCGTCAGGTCGGCAAGCGTGGTATCGACGCCATTTCCCAGAGTCAATGTCGCGCTATCGGTAATCCGGACCGCCCCATCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTcgcg                                                                                                                                          <  1:2170556/150‑1 (MQ=255)
            ccGGCTGGCAGCCGGCAGGTCGGCAAGCGTGGTATCGACGCCATTTCCCAGAGTCAATGTCGCGCTATCGGTAATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGAAGCCTGTcgttaa                                                                                                                                        <  2:1070026/150‑5 (MQ=255)
             cGGCTGGCAGCCGTCAGGTCGGCAAGCGTGGTATCGACGCCATTTCCCAGAGTCAATGTCGCGCTATCGGTAATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGAAGCCTGTcgttaat                                                                                                                                       <  1:2144336/150‑6 (MQ=255)
              ggCTGGCAGCCGTCAGGTCGGCAAGCGTGGTATCGACGCCATTTCCCAGAGTCAATGTCGCGCTATCGGTAATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGAAGCCTGTcgttaatg                                                                                                                                      <  1:377360/150‑7 (MQ=255)
                cTGGCAGCCGTCAGGTCGGCAAGCGTGGTATCGACGCCATTTCCCAGAGTCAATGTCGCGCTATCGGTACTCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTg                                                                                                                                    <  1:383358/150‑1 (MQ=255)
                cTGGCAGCCGTCAGGTCGGCAAGCGTGGTATCGACGCCATTTCCCAGAGTCAATGTCGCGCTATCGGTAATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTg                                                                                                                                    <  2:1384764/150‑1 (MQ=255)
                cTGGCAGCCGTCAGGTCGGCAAGCGTGGTATCGACGCCATTTCCCAGAGTCAATGTCGCGCTATCGGTAATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTg                                                                                                                                    >  1:689957/1‑150 (MQ=255)
                   gCAGCCGTCAGGTCGGCAAGCTTGGTATCGACGCCATTTCCCAGAGTCAATGTCGCGCTATCGGTTATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCAGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGAcc                                                                                                                                 <  2:1719685/150‑1 (MQ=255)
                   gCAGCCGTCAGGTCGGCAAGCGTGGTATCGACGCCATTTCCCAGAGTCAATGGCGCGCTATCGGTAATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGAcc                                                                                                                                 <  2:1907703/150‑1 (MQ=255)
                   gCAGCCGTCAGGGCGGCAAGCGTGGTAGCGACGCCATTTCCCAGAGTCACTGGCGCGCTATCGGTAATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCGGTCGCGCCACTGAcc                                                                                                                                 <  1:914049/150‑1 (MQ=255)
                    cAGCCGTCAGGTCGGCAAGCGTGGTATCGACGCCATTTCCCAGAGTCAATGTCGCGCTATCGGTAATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGCCCg                                                                                                                                >  2:63557/1‑150 (MQ=255)
                    cAGCCGTCAGGTCGGCAAGCGTGGTATCGACGCCATTTCCCAGAGTCAATGTCGCGCTATCGGTAATCCGGACCGCCCCTGCGAGTTTAACTGGCGTACCATTATCCCGTGGCGGCATTAACGACAGGCGTCCGGTCGCGCCACTGACCg                                                                                                                                >  1:438002/1‑150 (MQ=255)
                       ccGGCAGGTCGGCAAGCGTGGTATCGACGCCATTTCCCAGAGTCAAAGTCGCGCTATCGGTAATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATg                                                                                                                             <  2:1056600/150‑1 (MQ=255)
                        cGTCAGGTCGGCAAGCGTGGTATCGACGCCATTTCCCAGAGTCAATGTCGCGCTATCGGTAATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGCCCGATGc                                                                                                                            >  2:630376/1‑150 (MQ=255)
                           cAGGTCGGCAAGCGTGGTATCGACGCCATTTCCCAGAGTCAATGTCGCGCTATCGGTAATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGTCGTCAATAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATc                                                                                                                         >  2:633551/1‑150 (MQ=255)
                           cAGGTCGGCAAGCGTGGTATCGACGCCATTTCCCAGAGTCAATGTCGCGCTATCGGTAATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCAGTCGCGCCGCTGCCCGATGCAta                                                                                                                         >  1:2198581/1‑149 (MQ=255)
                             ggTCGGCAAGCGTGGTATCGACGACATTTCCCAGAGTCAATGTCGCGCTATCGGTAATCCGGACCGCCCCTTCTAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCg                                                                                                                       >  1:1094507/1‑150 (MQ=255)
                                 ggCAAGCGTGGTATCGACGCCATTTCCCAGAGTCAATGTCGCGCTATCGGTAATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAg                                                                                                                   >  2:1862285/1‑150 (MQ=255)
                                  gCAAGCGTGGTATCGACGCCATTTCCCAGAGTCAATGTCGCGCTATCGGTAATCCGGCCCGCCCCGTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGc                                                                                                                  >  1:917830/1‑150 (MQ=255)
                                   cAAGCGTGGTATCGACGCCATTTCCCAGAGTCAATGTCGCGCTATCGGTAATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTATCGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCg                                                                                                                 <  1:326009/150‑1 (MQ=255)
                                     aGCGTGTATCGACGCCATTTCCCAGAGTCAATGTCGCGCTATCGGTAATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTAc                                                                                                              >  2:1245867/1‑150 (MQ=255)
                                      gCGTGGTATCGACGCCATTTCCCAGAGTCAATGTCGCGCTATCGGTAATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTAc                                                                                                              >  1:1249763/1‑150 (MQ=255)
                                        gTGGTATCGACGCCATTTCCCAGAGTCAATGTCGCGCTATCGGTAATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCt                                                                                                            >  1:1010615/1‑150 (MQ=255)
                                         tGGTATCGACGCCATTTCCCAGAGTCAATGTCGCGCTCTCGGTAATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTg                                                                                                           <  1:1097992/150‑1 (MQ=255)
                                            tATCGACGCCATTTCCCAGAGTCAATGTCGCGCTATCGGTAATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGt                                                                                                        <  2:2106051/150‑1 (MQ=255)
                                             aTCGACGCCATTTCCCAGAGTCAATGTCGCGCTATCGGTAATCCGGACCGCCCCTTCGAGTTTAACTGGTGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTa                                                                                                       <  2:411383/150‑1 (MQ=255)
                                             aTCGACGCCATTTCCCAGAGACAATGTCGCGCTATCGGTAATCAGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGAGTAACTGCGtt                                                                                                       >  2:1140208/1‑149 (MQ=255)
                                              tCGACGCCATTTCCCAGAGTCAATGTCGCGCTATCGGTAATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAg                                                                                                      <  2:1334284/150‑1 (MQ=255)
                                              tCGACGCCATTTCCCAGAGTCAATGTCGCGCTATCGGTAATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAg                                                                                                      <  2:449516/150‑1 (MQ=255)
                                              tCGACGCCATTTCCCAGAGTCAATGTCGCGCTATCGGTAATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAg                                                                                                      <  1:1899602/150‑1 (MQ=255)
                                               cGACGCCATTTCCCAGAGTCAATGTCGCGCTATCGGTAATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCCCTGACCGATGCATCCGCCAGCGAACCTGCGTAGa                                                                                                     >  1:1416511/1‑150 (MQ=255)
                                                    ccATTTCCCAGAGTCAATGTCGCGCTATCGGTAATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGAtt                                                                                                >  1:1510882/1‑150 (MQ=255)
                                                    ccATTTCCCAGAGTCAATGTCGCGCTATCGGTAATCCGGACCGCCCCTCCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCAGTCTCGCCACTGCCCGATGCATCCGCCAGCGTACCTGAGTAGACGAtt                                                                                                >  2:556156/1‑150 (MQ=255)
                                                     cATTTCCCAGAGTCAATGTCGCGCTATCGGTAATCCGGACAGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGAGTACCTGCGTAGACGATTg                                                                                               >  2:1571089/1‑150 (MQ=255)
                                                      aTTTCCCAGAGTCAATGTCGCGCTATCGGTAATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAAGATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGATTgg                                                                                              <  2:992961/150‑2 (MQ=255)
                                                          cccAGAGTCAATGTCGCGCTATCGGTAATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTa                                                                                          <  1:457407/150‑1 (MQ=255)
                                                          cccAGAGTCAATGTCGCGCTATCGGTAATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTa                                                                                          >  2:157111/1‑150 (MQ=255)
                                                          cccAGAGTCAATGTCGCGCTATCGGTAATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTa                                                                                          >  1:2116964/1‑150 (MQ=255)
                                                          cccAGAGTCAATGTCGCGATATCGGTAATCCGGCCCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCAACCGCCAGCGTACCTGCGTAGACGATTGCCGTa                                                                                          >  2:2096151/1‑150 (MQ=255)
                                                           ccAGAGTCAATGTCGCGCTATCGGTAATCCGGCCCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCTTACCTGCGTAGACGATTGCCGTAc                                                                                         >  1:1368588/1‑150 (MQ=255)
                                                            cAGAGTCAATGTCGCGCTATCGGTAATCCGGCCCGCCCCTTCGAGTTTAACTGGCGTAACATTAACCCGTGGCGTCATTAACGACACGCATCCTGTCGCGCCCCTGCCCCATGCATCCGCCAGCCTACCTGCGTAGACGCTAGCCGTAcc                                                                                        >  2:913733/1‑150 (MQ=255)
                                                              gagTCAATGTCGCGCTATCGGTAATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTGATTTACGACCGGCTTCATGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGCTTGCCGTACCGc                                                                                      >  2:1962722/1‑150 (MQ=255)
                                                                gTCAATGTCGCGCTATCGGTAATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCa                                                                                    <  2:589017/150‑1 (MQ=255)
                                                                  cAATGTCGCGCTATCGGTAATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGc                                                                                  <  1:1750992/150‑1 (MQ=255)
                                                                  cAATGTCGCGCTATCGGTAATCCGGACCGCCCCTTCGAGTTTAACTAGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGc                                                                                  >  2:1515293/1‑150 (MQ=255)
                                                                      gTCGCGCTATCGGTAATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGCCCGATGCATCCGCCAGCGTCCCTGCGTAGACGATTGCCGTACCGCCCGCGAcc                                                                              >  1:500346/1‑150 (MQ=255)
                                                                           gcTATCGGTAATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGGCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCAGCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGgag                                                                         <  1:1109669/150‑1 (MQ=255)
                                                                            cgaTCGGTAATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGAGAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGgaga                                                                        <  2:376156/148‑1 (MQ=255)
                                                                            cTATCGGTAATCCGGACCGCCCCTTCGAGTTTCACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGgaga                                                                        <  2:330410/150‑1 (MQ=255)
                                                                             tATCGGTAATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGGAGAt                                                                       <  2:1010615/150‑1 (MQ=255)
                                                                              aTCGGTAATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGGAGAtc                                                                      <  2:1757237/150‑1 (MQ=255)
                                                                              aTCGGTAATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGGAGAtc                                                                      <  1:528756/150‑1 (MQ=255)
                                                                              aTCGGTAATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGGAGAtc                                                                      <  1:520621/150‑1 (MQ=255)
                                                                              aTCGGTAATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGGAGAtc                                                                      <  2:1662798/150‑1 (MQ=255)
                                                                              aTCGGTAATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGGAGAtc                                                                      <  1:1556679/150‑1 (MQ=255)
                                                                                 ggTAATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCCGCGTACCTGCGTAGACGCTTGCCGTCCCGCCCGCGACCTGGAGAactta                                                                   >  2:415986/1‑149 (MQ=255)
                                                                                      tCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGGAGATCTTCTGccc                                                              <  2:1923816/150‑1 (MQ=255)
                                                                                        cGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGGAGATCTTCTGCCCgg                                                            >  2:2034672/1‑150 (MQ=255)
                                                                                          gACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGGAGATCTTCTGCCCGGGc                                                          >  1:51227/1‑150 (MQ=255)
                                                                                           aCCGCCCCTTCGAGTTTAACTGGCGGAACATTATCCCGTGGCGGCATTAACGACAGGCCTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGGAGATCTTCTGCCCGGGcc                                                         <  1:307315/150‑1 (MQ=255)
                                                                                           aCCGCCCCTTCGAGTATAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGGAGATCTTCTGCCCGGGcc                                                         <  2:286225/150‑1 (MQ=255)
                                                                                                cccTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGGAGATCTTCTGCCCGGGCCAGAGc                                                    <  2:1094507/150‑1 (MQ=255)
                                                                                                 ccTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGGCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGGAGATCTTCTGCCCGGGCCAGAGcc                                                   <  1:1986212/150‑1 (MQ=255)
                                                                                                  cTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGGAGATCTTCTGCCCGGGCCAGAGCCt                                                  >  1:1447737/1‑150 (MQ=255)
                                                                                                   ttCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGGAGATCTTCTGCCCGGGCCAGAGCCTg                                                 <  1:157111/150‑1 (MQ=255)
                                                                                                   ttCGAGTTTAACTGGCGTAACATTATCCCGTGGCGGCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGGAGATCTTCTGCCCGGGCCAGAGCCTg                                                 <  2:1330250/150‑1 (MQ=255)
                                                                                                    tCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGGAGATCTTCTGCCCGGGCCAGAGCCTGa                                                <  1:1140208/150‑1 (MQ=255)
                                                                                                      gAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGGAGATCTTCTGCCCGGGCCAGAGCCTGaaa                                              <  1:503397/150‑1 (MQ=255)
                                                                                                       aGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGGAGATCTTCTGCCCGGGCCAGAGCCTGAAAc                                             <  2:135717/150‑1 (MQ=255)
                                                                                                           tAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGGAGATCTTCTGCCCGGGCCAGAGCCTGAAACTCaa                                         >  2:306070/1‑149 (MQ=255)
                                                                                                           tAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGGAGATCTTCTGCCCGGGCCAGAGCCTGAAACTCAt                                         <  1:1420978/150‑1 (MQ=255)
                                                                                                             ccTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCGTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGGAGATCTTCTGCCCGGGCCAGAGCCTGAAACTCATCt                                       <  2:1816843/149‑1 (MQ=255)
                                                                                                              cTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGAACGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGGAGATCTTCTGCCCGGGCCAGAGCCTGAAACTCATCtt                                      >  1:1531587/1‑150 (MQ=255)
                                                                                                               tGGCGTAACATTATCCCGTGGCGGCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCTCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGGAGATCTTCTGCCCGGGCCAGAGCCTGAAACTCATCttt                                     <  2:1839520/150‑1 (MQ=255)
                                                                                                                ggCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGGAGATCTTCTGCCCGGGCCAGAGCCTGAAACTCATCtttt                                    >  2:459096/1‑150 (MQ=255)
                                                                                                                 gagTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGGAGATCTTCTGCCCGGGCCAGAGCCTGAAACTCATCTTTTg                                   <  2:1085326/148‑1 (MQ=255)
                                                                                                                   gTAACATTATCCCGTGGCGTCATTAACGCCCGGCTTCCTGTCGCGCCACTGACCGAAGCATCCGCCAGAGTACCTGCGTAGACGATGGCCGTTCCGCCCGCGACCTCGAGATCTTCTGCCCGGTCCAGAGCCTGAAACTCATCTGTTGAc                                 >  2:1441975/1‑150 (MQ=255)
                                                                                                                   gTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCTCCACTGACCGATGCATCAGCCAGAGTACCTGCGTAGACGAATGCCGTACCGCCAGCCACCTGGAGATCTTCTGACCGGGCCAGAGCCTGAAACGAAACTTTTGAc                                 >  2:1441978/1‑150 (MQ=255)
                                                                                                                     aaCATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCCGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGGAGCTCTTCTGCCCGGGCCAGAGCCTGAAGCTGATCTGTTGACCg                               >  1:314627/1‑150 (MQ=255)
                                                                                                                      aCATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGGAGATCTTCTGCCCGGGCCAGAGCCTGAAACTCATCTTTTGACCGc                              <  1:1324151/150‑1 (MQ=255)
                                                                                                                       cATTATCCCGTGGCGTCATTAACGATAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCGGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGGAGATCTTCTGCCCGGGCCAGAGCCTGAAACTCATCTTTTGACCGcc                             >  1:1478780/1‑150 (MQ=255)
                                                                                                                       cATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGGAGATCTTCTGCCCGGGCCAGCGCCTGAAACTCATCTTTTGACCGcc                             >  1:766922/1‑150 (MQ=255)
                                                                                                                       cATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGGAGATCTTCTGCCCGGGCCAGAGCCTGAAACTCATCTTTTGACCGcc                             <  2:1042160/150‑1 (MQ=255)
                                                                                                                       cATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGGAGATCTTCTGCCCGGGCCAGAGCCTGAAACTCATCTTTTGACCGcc                             >  1:791463/1‑150 (MQ=255)
                                                                                                                        aTTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGGAGATCTTCTGCCCGGGCCAGAGCCTGAAACTCATCTTTTGACCGccc                            <  1:1519426/150‑1 (MQ=255)
                                                                                                                                    cGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGGAGATCTTCTGCCCGGGCCAGAGCCTGAAACTCATCTTTTGACCGCCCAAGGGTATATTg                <  2:2116964/150‑1 (MQ=255)
                                                                                                                                       cATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGGAGATCTTCTGCCCGGGCCAGAGCCTGAAACTCATCTTTTGACCGCCCAAGGGTATATTGccc             >  1:1276881/1‑150 (MQ=255)
                                                                                                                                        tttAACGACAGGGTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGGAGATCTTCTGCCCGGGCCAGAGCCTGAAACTCATCTTTTGACCGCCCAAGGGTATATTGCCCg            <  2:1640721/149‑1 (MQ=255)
                                                                                                                                        aTTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGGAGATCTTCTGCCCGGGCCAGAGCCTGAAACTCATCTTTTGACCGCCCAAGGGTATATTGCCCg            <  2:276629/150‑1 (MQ=255)
                                                                                                                                        aTTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGGAGATCTTCTGCCCGGGCCAGAGCCTGAAACTCATCTTTTGACCGCCCAAGGGTATATTGCCCg            <  2:519178/150‑1 (MQ=255)
                                                                                                                                         ggcgcGACAGGAAGCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGGAGATCTTCTGCCCGGGCCAGAGCCTGAAACTCATCTTTTGACCGCCCAAGGGTATATTGCCCGc           >  2:995620/5‑150 (MQ=255)
                                                                                                                                          tAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGGAGATCTTCTGCCCGGGCCAGAGCCTGAAACTCATCTTTTGACCGCCCAAGGGTATATTGCCCGcc          >  1:1442246/1‑150 (MQ=255)
                                                                                                                                                  ggCTTCCTCTCGCGCCACTGACCGATGCATCCGCCAGCGAACCTCCGTAGACGATTGCCGTACCGCCAGGCACCTCGAGATCTTCTGCCCGCGCCAGAGCCTGAAACTCATCTAATCACCGCCCAAGGGTATATTGCCCGCCAGAGTTaa  >  2:1292112/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CCAGACACTGCCCCGGCTGGCAGCCGTCAGGTCGGCAAGCGTGGTATCGACGCCATTTCCCAGAGTCAATGTCGCGCTATCGGTAATCCGGACCGCCCCTTCGAGTTTAACTGGCGTAACATTATCCCGTGGCGTCATTAACGACAGGCTTCCTGTCGCGCCACTGACCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGGAGATCTTCTGCCCGGGCCAGAGCCTGAAACTCATCTTTTGACCGCCCAAGGGTATATTGCCCGCCAGAGTTAA  >  NC_000913/2779607‑2779902

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