Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 153,933 A→C 5.1% N498K (AAT→AAG) ‡ htrE ← putative fimbrial usher protein HtrE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913153,9330AC 5.1% 369.2 / 7.4 114N498K (AAT→AAG) ‡htrEputative fimbrial usher protein HtrE
Reads supporting (aligned to +/- strand):  ref base A (50/58);  new base C (2/4);  total (52/62)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 115 total aligned reads displayed.

GTAGCCAATAGAGTAATTGCTACGATTTTGTCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGT  >  NC_000913/153786‑154081
                                                                                                                                                   |                                                                                                                                                    
gtAGCCAATAGAGTAATTGCTACGATTTTGTCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAAtt                                                                                                                                                    <  2:1053060/150‑1 (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:192349/1‑150 (MQ=255)
      aaTAGAGTAATTGCTACGATTTTGTCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCAt                                                                                                                                              <  1:201236/150‑1 (MQ=255)
       aTAGAGTAATTGCTACGATTTTGTCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATg                                                                                                                                             >  1:933276/1‑150 (MQ=255)
         agagTAATTGCTACGATTTTGTCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTAATTACTCACGCATgaa                                                                                                                                           <  1:1209569/150‑3 (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:1460781/1‑150 (MQ=255)
                   cTACGATTTTGTCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCa                                                                                                                                 >  2:1428173/1‑150 (MQ=255)
                      cGATTTTGTCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGa                                                                                                                              >  1:664310/1‑150 (MQ=255)
                      cGATTTTGTCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGa                                                                                                                              >  1:942097/1‑150 (MQ=255)
                      cGATTTTGTCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGa                                                                                                                              >  1:988046/1‑150 (MQ=255)
                       gATTTTGTCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGAt                                                                                                                             <  1:921689/150‑1 (MQ=255)
                             gTCCGGAAGCCCAGTAATCGGACCAACTTCCTGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGtt                                                                                                                       <  1:1478683/150‑1 (MQ=255)
                             gTCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGtt                                                                                                                       >  1:901637/1‑150 (MQ=255)
                                cGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCgg                                                                                                                    <  1:1487032/150‑1 (MQ=255)
                                 ggAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCggg                                                                                                                   >  1:1216083/1‑150 (MQ=255)
                                      cccAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGtt                                                                                                              <  2:936157/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)
                                           tAATCGGACCAACTTGCGGAAAGATAAAATGAAGCTTATTCTTTTTTCTCAAATTTCAACGTTTGGTTAATACTGACCGTAAGCTTATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACt                                                                                                         <  2:1262901/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)
                                                       cTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTaga                                                                                             >  1:630823/1‑150 (MQ=255)
                                                        ttCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTagag                                                                                            <  1:375783/150‑1 (MQ=255)
                                                         tCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGt                                                                                           >  1:1423686/1‑150 (MQ=255)
                                                          ccGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGtt                                                                                          >  2:49885/1‑150 (MQ=255)
                                                          ccGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGtt                                                                                          <  1:238679/150‑1 (MQ=255)
                                                          ccGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGtt                                                                                          <  2:933276/150‑1 (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)
                                                                     aaaTGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCAtt                                                                               >  2:100747/1‑150 (MQ=255)
                                                                      aaTGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTa                                                                              <  2:664310/150‑1 (MQ=255)
                                                                      aaTGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTa                                                                              <  1:578765/150‑1 (MQ=255)
                                                                            ccGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCa                                                                        >  1:1517450/1‑150 (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:306968/1‑150 (MQ=255)
                                                                                     tttttCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTc                                                               <  2:1098390/150‑1 (MQ=255)
                                                                                     tttttCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTc                                                               <  2:589022/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)
                                                                                                   cAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCg                                                 <  1:49885/150‑1 (MQ=255)
                                                                                                     aCGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGAt                                               <  1:173783/150‑1 (MQ=255)
                                                                                                     aCGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGAt                                               <  1:17534/150‑1 (MQ=255)
                                                                                                     aCGGATGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGAt                                               <  2:1586665/150‑1 (MQ=255)
                                                                                                      cGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATa                                              >  2:1395132/1‑150 (MQ=255)
                                                                                                      cGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATa                                              <  2:630823/150‑1 (MQ=255)
                                                                                                           ggttAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCcg                                         >  1:1626100/1‑150 (MQ=255)
                                                                                                            gttAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTc                                                         <  2:1339767/133‑1 (MQ=255)
                                                                                                            gttAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTc                                                         >  1:1339767/1‑133 (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)
                                                                                                                      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)
                                                                                                                              cTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCTATATTCAGTGAAGTACt                      <  2:1423686/150‑1 (MQ=255)
                                                                                                                              cTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACt                      >  2:531894/1‑150 (MQ=255)
                                                                                                                              cTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACt                      >  2:1053282/1‑150 (MQ=255)
                                                                                                                                gATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTg                    <  2:273926/150‑1 (MQ=255)
                                                                                                                                gATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTg                    <  1:1457417/150‑1 (MQ=255)
                                                                                                                                   aatccATGCGTAATTACTCACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACAAAGGTAATTCTGTGTCGATTAGCGATAGGCGGCGATTTTAAGTGAAGTACTTGttt                 >  1:1310157/5‑150 (MQ=255)
                                                                                                                                    atccATGCGTAATTACTCACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTttc                >  2:552583/4‑150 (MQ=255)
                                                                                                                                        aTGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCg            <  2:29628/150‑1 (MQ=255)
                                                                                                                                        aTGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCg            <  2:1379670/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)
                                                                                                                                               aGTAATTACGCATGGATATCGGTTCAACATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAATGTAATTCTGTGGCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACt     <  2:1545082/150‑1 (MQ=255)
                                                                                                                                                  aaTTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGt  >  1:241085/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GTAGCCAATAGAGTAATTGCTACGATTTTGTCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGT  >  NC_000913/153786‑154081

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: