Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 153,928 G→A 5.3% S500L (TCA→TTA)  htrE ← putative fimbrial usher protein HtrE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913153,9280GA 5.3% 370.7 / 5.7 115S500L (TCA→TTA) htrEputative fimbrial usher protein HtrE
Reads supporting (aligned to +/- strand):  ref base G (50/58);  new base A (2/4);  total (52/63)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

Only 100 of 118 total aligned reads displayed.

TGTTACTGTAGCCAATAGAGTAATTGCTACGATTTTGTCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAAC  >  NC_000913/153779‑154077
                                                                                                                                                     |                                                                                                                                                     
  ttactgtAGCCAATAGAGTAATTGCTACGATTTTGTCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAg                                                                                                                                                     <  2:223504/150‑1 (MQ=255)
  ttactgtAGCCAATAGAGTAATTGCTACGATTTTGTCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAg                                                                                                                                                     <  2:365420/150‑1 (MQ=255)
    actgtAGCGAATAGAGTAATTGCTACGATTTTGTCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAgca                                                                                                                                                   >  2:1385943/1‑148 (MQ=255)
       gtAGCCAATAGAGTAATTGCTACGATTTTGTCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAAtt                                                                                                                                                >  2:238679/1‑150 (MQ=255)
       gtAGCCAATAGAGTAATTGCTACGATTTTGTCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAAtt                                                                                                                                                >  2:578765/1‑150 (MQ=255)
         aGCCAATAGAGTAATTGCTACGATTTTGTCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTAc                                                                                                                                              >  2:1457417/1‑150 (MQ=255)
            cAATAGAGTAATTGCTACGATTTTGTCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCa                                                                                                                                           >  1:1172056/1‑150 (MQ=255)
             aaTAGAGTAATTGCTACGATTTTGTCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCAt                                                                                                                                          <  1:201236/150‑1 (MQ=255)
             aaTAGAGTAATTGCTACGATTTTGTCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCAt                                                                                                                                          >  1:192349/1‑150 (MQ=255)
              aTAGAGTAATTGCTACGATTTTGTCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATg                                                                                                                                         >  1:933276/1‑150 (MQ=255)
               tAGAGTAATTGCTACGATTTTGTCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACGTGATTTTTCATGCGTAATTACTCACGCATga                                                                                                                                        <  1:202466/150‑2 (MQ=255)
                agagTAATTGCTACGATTTTGTCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTAATTACTCACGCATgaa                                                                                                                                       <  1:1209569/150‑3 (MQ=255)
                 cagTAATTGCTACGATTTTGTCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATCAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTAATTACTCACGCATgaaa                                                                                                                                      <  2:1317404/149‑4 (MQ=255)
                  agTAATTGCTACGATTTTGTCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGAtt                                                                                                                                     >  1:1098390/1‑150 (MQ=255)
                     aaTTGCTACGATTTTGTCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTCCATGCGTAATTACTCACGCATgaaaaatc                                                                                                                                  <  2:976910/150‑8 (MQ=255)
                       ttGCTACGATTTTGTCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGt                                                                                                                                <  1:455192/150‑1 (MQ=255)
                        tGCTACGATTTTGTCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGtt                                                                                                                               <  2:1175361/150‑1 (MQ=255)
                         gCTACGATTTTGTCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTc                                                                                                                              >  2:593240/1‑150 (MQ=255)
                         gCTACGATTTTGTCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTc                                                                                                                              >  2:1460781/1‑150 (MQ=255)
                          cTACGATTTTGTCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCa                                                                                                                             >  2:1428173/1‑150 (MQ=255)
                             cGATTTTGTCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGa                                                                                                                          >  1:988046/1‑150 (MQ=255)
                             cGATTTTGTCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGa                                                                                                                          >  1:942097/1‑150 (MQ=255)
                              gATTTTGTCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGAt                                                                                                                         <  1:921689/150‑1 (MQ=255)
                                ttttGTCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCt                                                                                                                       <  1:1552014/150‑1 (MQ=255)
                                   tGTCCGGAAGCCCAATAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTa                                                                                                                                                   <  1:93817/119‑1 (MQ=255)
                                   tGTCCGGAAGCCCAATAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTa                                                                                                                                                   >  2:93817/1‑119 (MQ=255)
                                    gTCCGGAAGCCCAGTAATCGGACCAACTTCCTGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGtt                                                                                                                   <  1:1478683/150‑1 (MQ=255)
                                    gTCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGtt                                                                                                                   >  1:901637/1‑150 (MQ=255)
                                       cGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCgg                                                                                                                <  1:1487032/150‑1 (MQ=255)
                                        ggAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCggg                                                                                                               <  2:1100059/150‑1 (MQ=255)
                                          aaGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGAt                                                                                                             <  1:1219396/150‑1 (MQ=255)
                                               cAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTc                                                                                                        >  1:843401/1‑150 (MQ=255)
                                                aGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCa                                                                                                       >  2:111342/1‑150 (MQ=255)
                                                 gTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCAc                                                                                                      <  1:23896/150‑1 (MQ=255)
                                                  tAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACt                                                                                                     <  2:762372/150‑1 (MQ=255)
                                                   aaTCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACtt                                                                                                    >  1:235623/1‑150 (MQ=255)
                                                   aaTCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACtt                                                                                                    >  1:390286/1‑150 (MQ=255)
                                                    aTCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACt                                                                                                     >  2:1334575/1‑148 (MQ=255)
                                                    aTCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACt                                                                                                     <  1:1334575/148‑1 (MQ=255)
                                                        gACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTtcatca                                                                                               <  1:1543319/150‑1 (MQ=255)
                                                         aCCAACTTCCGGAAAGATAAAATTAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGTGATGTTTCACTTCATCaa                                                                                              <  1:1565764/150‑1 (MQ=255)
                                                           cAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAAtt                                                                                            <  1:605793/150‑1 (MQ=255)
                                                            aaCTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTa                                                                                           >  1:1441268/1‑150 (MQ=255)
                                                              cTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTaga                                                                                         >  1:630823/1‑150 (MQ=255)
                                                               ttCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTagag                                                                                        <  1:375783/150‑1 (MQ=255)
                                                                 ccGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGtt                                                                                      <  2:927605/150‑1 (MQ=255)
                                                                 ccGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGtt                                                                                      <  2:933276/150‑1 (MQ=255)
                                                                 ccGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGtt                                                                                      >  2:49885/1‑150 (MQ=255)
                                                                   ggAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCTTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAg                                                                                    <  2:980720/150‑1 (MQ=255)
                                                                    gAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGt                                                                                   <  2:1303703/150‑1 (MQ=255)
                                                                     aaaGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTg                                                                                  <  1:1566447/150‑1 (MQ=255)
                                                                       aGATAATATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTCGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGca                                                                                <  2:634714/150‑1 (MQ=255)
                                                                       aGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGca                                                                                <  1:915840/150‑1 (MQ=255)
                                                                        gATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGcat                                                                               >  2:8522/1‑150 (MQ=255)
                                                                          tAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCt                                                                                                                       >  2:428401/1‑108 (MQ=255)
                                                                          tAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCt                                                                                                                       <  1:428401/108‑1 (MQ=255)
                                                                            aaaTGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCAtt                                                                           >  2:100747/1‑150 (MQ=255)
                                                                             aaTGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTa                                                                          <  1:578765/150‑1 (MQ=255)
                                                                                     gTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAg                                                                  <  1:100747/150‑1 (MQ=255)
                                                                                     gTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAg                                                                  >  2:1291200/1‑150 (MQ=255)
                                                                                      tAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTACAGTTAGTGCATCATTAAGACCAAgg                                                                 <  1:1355424/150‑1 (MQ=255)
                                                                                       aaTCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGt                                                                <  2:390286/150‑1 (MQ=255)
                                                                                            tttttCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTc                                                           <  2:589022/150‑1 (MQ=255)
                                                                                            tttttCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTc                                                           <  2:1098390/150‑1 (MQ=255)
                                                                                                 ctcAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTgtgt                                                      >  1:1536014/1‑150 (MQ=255)
                                                                                                  tcAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTc                                                     >  2:148753/1‑150 (MQ=255)
                                                                                                   cAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCg                                                    <  2:942097/150‑1 (MQ=255)
                                                                                                            aCGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGAt                                           <  1:17534/150‑1 (MQ=255)
                                                                                                            aCGGATGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGAt                                           <  2:1586665/150‑1 (MQ=255)
                                                                                                             cGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATa                                          <  2:630823/150‑1 (MQ=255)
                                                                                                             cGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATa                                          >  2:1395132/1‑150 (MQ=255)
                                                                                                                  ggttAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCcg                                     >  1:1626100/1‑150 (MQ=255)
                                                                                                                   gttAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTc                                                     >  1:1339767/1‑133 (MQ=255)
                                                                                                                   gttAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTc                                                     <  2:1339767/133‑1 (MQ=255)
                                                                                                                   gttAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCcgc                                    <  2:901637/150‑1 (MQ=255)
                                                                                                                    ttAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCcgcg                                   >  1:292566/1‑150 (MQ=255)
                                                                                                                     tAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGa                                  >  1:59106/1‑150 (MQ=255)
                                                                                                                        tACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGatat                               >  2:1394328/1‑150 (MQ=255)
                                                                                                                         aCTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATAtt                              >  1:1251492/1‑150 (MQ=255)
                                                                                                                          cTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTc                             <  1:1166924/150‑1 (MQ=255)
                                                                                                                           tGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCa                            <  2:141114/150‑1 (MQ=255)
                                                                                                                             aCCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGt                          <  1:660232/150‑1 (MQ=255)
                                                                                                                             aCCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGt                          >  2:1320537/1‑150 (MQ=255)
                                                                                                                              ccGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTg                         >  2:153464/1‑150 (MQ=255)
                                                                                                                              ccGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTg                         >  2:11217/1‑150 (MQ=255)
                                                                                                                                 tAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAg                      <  2:292566/150‑1 (MQ=255)
                                                                                                                                    ccTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTAc                   <  1:1428173/150‑1 (MQ=255)
                                                                                                                                     cTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACt                  >  2:1053282/1‑150 (MQ=255)
                                                                                                                                     cTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACt                  >  2:531894/1‑150 (MQ=255)
                                                                                                                                     cTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACt                  >  2:369549/1‑150 (MQ=255)
                                                                                                                                       gATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTg                <  1:1457417/150‑1 (MQ=255)
                                                                                                                                       gATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTg                <  2:273926/150‑1 (MQ=255)
                                                                                                                                          aatccATGCGTAATTACTCACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACAAAGGTAATTCTGTGTCGATTAGCGATAGGCGGCGATTTTAAGTGAAGTACTTGttt             >  1:1310157/5‑150 (MQ=255)
                                                                                                                                               aTGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCg        <  2:29628/150‑1 (MQ=255)
                                                                                                                                               aTGCGGCAGTAATTACGCATGGNTTTCGGTTCAAGATGTTGTTCGNGATGTTTCACTTCATCAATTAGAGTTAGTNNNTGATTAAGTCCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCg        <  2:1385375/150‑1 (MQ=255)
                                                                                                                                                   gTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAg                      >  2:721655/1‑132 (MQ=255)
                                                                                                                                                   gTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAg                      <  1:721655/132‑1 (MQ=255)
                                                                                                                                                   gTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATa    >  2:956663/1‑150 (MQ=255)
                                                                                                                                                     gAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAAc  >  1:922083/1‑150 (MQ=255)
                                                                                                                                                     gAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAAc  >  1:1157337/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGTTACTGTAGCCAATAGAGTAATTGCTACGATTTTGTCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAAC  >  NC_000913/153779‑154077

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: