Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 302,225 T→A 5.8% N117Y (AAT→TAT)  paoA ← aldehyde dehydrogenase, Fe‑S subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913302,2250TA 5.8% 265.0 / 5.3 86N117Y (AAT→TAT) paoAaldehyde dehydrogenase, Fe‑S subunit
Reads supporting (aligned to +/- strand):  ref base T (39/42);  new base A (3/2);  total (42/44)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.31e-01

ATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAGCTGCTCGGTTTTGCCGTTCA  >  NC_000913/302077‑302374
                                                                                                                                                    |                                                                                                                                                     
aTTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTa                                                                                                                                                      <  2:149831/150‑1 (MQ=255)
aTTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTa                                                                                                                                                      <  2:1259312/150‑1 (MQ=255)
aTTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTa                                                                                                                                                      >  2:1084199/1‑150 (MQ=255)
aTTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATAATCTGGCGAACCCAGGCCTTCAATGGTGGGGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTa                                                                                                                                                      <  1:662526/150‑1 (MQ=255)
  ttGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAg                                                                                                                                                    <  2:411502/150‑1 (MQ=255)
    gCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGcc                                                                                                                                                  <  1:1084199/150‑1 (MQ=255)
      ccGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTg                                                                                                                                                <  2:1151408/150‑1 (MQ=255)
          aGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCTTAATGCCTGcctg                                                                                                                                            <  2:453676/150‑4 (MQ=255)
          aGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGAc                                                                                                                                            <  2:696477/150‑1 (MQ=255)
          aGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGAc                                                                                                                                            <  1:760745/150‑1 (MQ=255)
           ggTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCTTAATGCCTGcctga                                                                                                                                           <  1:681988/150‑5 (MQ=255)
              gCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCAtt                                                                                                                                        >  2:1465023/1‑150 (MQ=255)
              gCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGACGACTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATGTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCAtt                                                                                                                                        >  1:416203/1‑150 (MQ=255)
               cAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTg                                                                                                                                       <  2:913447/150‑1 (MQ=255)
               cAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGAGCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTg                                                                                                                                       <  2:507979/150‑1 (MQ=255)
                aGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGa                                                                                                                                      >  1:192763/1‑150 (MQ=255)
                aGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTCCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGCCCATTgc                                                                                                                                      >  1:1103263/1‑149 (MQ=255)
                 gTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGAATGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGAc                                                                                                                                     >  2:983611/1‑150 (MQ=255)
                  tAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGAcc                                                                                                                                    <  1:787183/150‑1 (MQ=255)
                  tAGCCGCACTGGAAGCCATCATGCGTGATAAAGGCCGCCGGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGAGCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGAcc                                                                                                                                    <  1:876663/150‑1 (MQ=255)
                   aGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCa                                                                                                                                   <  2:743955/150‑1 (MQ=255)
                      cGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCa                                                                                                                                >  1:1325648/1‑150 (MQ=255)
                      cGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCa                                                                                                                                >  1:1337627/1‑150 (MQ=255)
                       gCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCAc                                                                                                                               <  2:423410/150‑1 (MQ=255)
                               aGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAgg                                                                                                                       >  1:968641/1‑150 (MQ=255)
                                   atcatGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTcc                                                                                                                   <  1:1606441/150‑1 (MQ=255)
                                         cTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTg                                                                                                             >  1:722920/1‑150 (MQ=255)
                                            gATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTcc                                                                                                          >  1:1126506/1‑150 (MQ=255)
                                                aaGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGa                                                                                                      >  2:899375/1‑150 (MQ=255)
                                                aaGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGa                                                                                                      <  2:449503/150‑1 (MQ=255)
                                                   gccgccTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCg                                                                                                   <  1:1502459/150‑1 (MQ=255)
                                                   gccgccTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCCTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCg                                                                                                   <  2:134203/150‑1 (MQ=255)
                                                       ccTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAAc                                                                                               <  1:54407/150‑1 (MQ=255)
                                                           cATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCttt                                                                                           >  1:622006/1‑150 (MQ=255)
                                                             tGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCt                                                                                         >  1:391876/1‑150 (MQ=255)
                                                              ggggTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCtt                                                                                        <  2:436252/150‑1 (MQ=255)
                                                                 gTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGt                                                                                     <  2:1020278/150‑1 (MQ=255)
                                                                  tGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGAAATACTAAAAACGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTa                                                                                    >  2:1382340/1‑150 (MQ=255)
                                                                       aTTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGAt                                                                               <  1:1147820/150‑1 (MQ=255)
                                                                         tATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCa                                                                             <  2:528081/150‑1 (MQ=255)
                                                                          aTCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCaa                                                                            >  2:354551/1‑150 (MQ=255)
                                                                            cTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAAcc                                                                                              >  1:1463443/1‑130 (MQ=255)
                                                                            cTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAAcc                                                                                              <  2:1463443/130‑1 (MQ=255)
                                                                             tGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATg                                                                         >  1:907669/1‑150 (MQ=255)
                                                                              ggCGAACCCAGGCCTTCAATGGTGGAGATCTCGGCCCACTGATGCAAAACTGAAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGc                                                                        >  2:397978/1‑150 (MQ=255)
                                                                               gCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGTGACCATTGAGCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCa                                                                       <  1:556625/150‑1 (MQ=255)
                                                                                   aCCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGAtt                                                                   <  2:977489/150‑1 (MQ=255)
                                                                                    cccAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGAttt                                                                  >  1:560935/1‑150 (MQ=255)
                                                                                    cccAGGCCTTCAATGGTGGTGATCTCGGCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGAtttt                                                                 >  1:1176657/1‑150 (MQ=255)
                                                                                      cAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATc                                                                              <  2:603909/136‑1 (MQ=255)
                                                                                      cAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATc                                                                              >  1:603909/1‑136 (MQ=255)
                                                                                           cTTCAATGGTGGAGATCTCGGCCCACTGATGCATGACTCCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCAGTGTCCCTGATCGCAACCTTTCTTGGTACCGATCAAAAGCAGAGTTTCACGCa                                                           >  2:1553783/1‑150 (MQ=255)
                                                                                             tCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCaaa                                                         <  2:942435/150‑1 (MQ=255)
                                                                                                   gtggtgATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTcc                                                   >  1:205040/1‑150 (MQ=255)
                                                                                                   gtggtgATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTcc                                                   <  1:121882/150‑1 (MQ=255)
                                                                                                    tggtgATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCa                                                  >  1:1339992/1‑150 (MQ=255)
                                                                                                        gATCTCGGCCCCATGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACAATTGACCAGCACGGTACAGGCTCAGCACAGTCCGTGATCGCAACCTATCTAGGTACCGATCTAATGCAGATTTTCACGCAAAGTGTCCAGTAg                                              >  2:427015/1‑150 (MQ=255)
                                                                                                           ctcGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGt                                           >  2:1232815/1‑150 (MQ=255)
                                                                                                            tcGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTg                                          >  2:162231/1‑150 (MQ=255)
                                                                                                             cGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTgg                                         >  1:1627512/1‑150 (MQ=255)
                                                                                                                 ccccTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATTGCAACCTTTCTTGGTACCTATCAAATGCAGATTTTCAAGCAAAGTGTCCAGTATCGTGGTTCg                                     <  2:192763/150‑1 (MQ=255)
                                                                                                                     tGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTa                                 >  1:711705/1‑150 (MQ=255)
                                                                                                                        tGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATcc                              <  2:227475/150‑1 (MQ=255)
                                                                                                                         gCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCa                             >  2:758430/1‑150 (MQ=255)
                                                                                                                           ctGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGCTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCAcc                           <  2:550750/149‑1 (MQ=255)
                                                                                                                                 gCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAgc                     <  1:555341/150‑1 (MQ=255)
                                                                                                                                  cAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAgct                    >  1:198359/1‑150 (MQ=255)
                                                                                                                                    aGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAgctgc                  <  2:1325648/150‑1 (MQ=255)
                                                                                                                                     gCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAgctgct                 <  2:711705/150‑1 (MQ=255)
                                                                                                                                         gtcgcaggctTAATGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAGCTGCTCGGt             >  2:539917/4‑150 (MQ=255)
                                                                                                                                         gtcgcaggctTAATGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAGCTGCTCGGt             >  1:398697/4‑150 (MQ=255)
                                                                                                                                         caggcaggCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCgg                                    <  1:917085/127‑1 (MQ=255)
                                                                                                                                         caggcaggCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCgg                                    >  2:917085/1‑127 (MQ=255)
                                                                                                                                          tcgcaggctTAATGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCTACCTCAAGCTGCTCGGtt            >  1:779498/3‑150 (MQ=255)
                                                                                                                                             caggcaTTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAGCTGCTCGGTTTTg         >  2:182317/1‑150 (MQ=255)
                                                                                                                                              aggcaTTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAGCTGCTCGGTTTTGc        <  2:622006/150‑1 (MQ=255)
                                                                                                                                              aggcaTTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAGCTGCTCGGTTTTGc        <  1:1142748/150‑1 (MQ=255)
                                                                                                                                              aggcaTTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAGCTGCTCGGTTTTGc        <  2:895993/150‑1 (MQ=255)
                                                                                                                                               ggcaTTAATCATGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAAACATTCTTGGTACCGAACAAATGAAGATTTTCACGCAAAGTGTCAAGTAGCGTGGTTCGGGTATCCACCTCAAGCTGCTAGGTTTTGca       >  2:1071918/1‑149 (MQ=255)
                                                                                                                                               ggcaTTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAGCTGCTCGGTTTTGcc       >  1:255053/1‑150 (MQ=255)
                                                                                                                                                gcaTTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAGCTGCTCGGTTTTGCCg      <  2:513568/150‑1 (MQ=255)
                                                                                                                                                gcaTTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAGCTGCTCGGTTTTGCCg      <  1:899375/150‑1 (MQ=255)
                                                                                                                                                 caTTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTa                                 <  1:1324117/122‑1 (MQ=255)
                                                                                                                                                 caTTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTa                                 >  2:1324117/1‑122 (MQ=255)
                                                                                                                                                 caTTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAGCTGCTCGGTTTTGCCGt     >  1:1198773/1‑150 (MQ=255)
                                                                                                                                                 caTTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAGCTGCTCGGTTTTGCCGt     >  2:690667/1‑150 (MQ=255)
                                                                                                                                                    tAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAGCTGCTCGGTTTTGCCGTTCa  >  2:760103/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAGCTGCTCGGTTTTGCCGTTCA  >  NC_000913/302077‑302374

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: