Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 153,930 G→T 5.2% Y499* (TAC→TAA htrE ← putative fimbrial usher protein HtrE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913153,9300GT 5.2% 395.0 / 7.0 115Y499* (TAC→TAAhtrEputative fimbrial usher protein HtrE
Reads supporting (aligned to +/- strand):  ref base G (51/58);  new base T (2/4);  total (53/62)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.58e-01

Only 100 of 119 total aligned reads displayed.

TTACTGTAGCCAATAGAGTAATTGCTACGATTTTGTCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACT  >  NC_000913/153781‑154078
                                                                                                                                                     |                                                                                                                                                    
  actgtAGCGAATAGAGTAATTGCTACGATTTTGTCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAgca                                                                                                                                                    >  2:1385943/1‑148 (MQ=255)
     gtAGCCAATAGAGTAATTGCTACGATTTTGTCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAAtt                                                                                                                                                 <  2:1053060/150‑1 (MQ=255)
     gtAGCCAATAGAGTAATTGCTACGATTTTGTCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAAtt                                                                                                                                                 >  2:238679/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:927605/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)
                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)
                        cTACGATTTTGTCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCa                                                                                                                              >  2:1428173/1‑150 (MQ=255)
                           cGATTTTGTCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGa                                                                                                                           >  1:664310/1‑150 (MQ=255)
                           cGATTTTGTCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGa                                                                                                                           >  1:942097/1‑150 (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)
                                      ggAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCggg                                                                                                                <  2:1100059/150‑1 (MQ=255)
                                        aaGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGAt                                                                                                              <  1:1219396/150‑1 (MQ=255)
                                             cAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTc                                                                                                         >  1:843401/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: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)
                                                             ttCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTagag                                                                                         <  1:375783/150‑1 (MQ=255)
                                                              tCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGt                                                                                        >  1:1423686/1‑150 (MQ=255)
                                                               ccGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGtt                                                                                       <  1:238679/150‑1 (MQ=255)
                                                               ccGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGtt                                                                                       <  2:933276/150‑1 (MQ=255)
                                                               ccGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGtt                                                                                       <  2:927605/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)
                                                                          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)
                                                                                     aaTCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGt                                                                 <  2:390286/150‑1 (MQ=255)
                                                                                          tttttCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTc                                                            <  2:589022/150‑1 (MQ=255)
                                                                                          tttttCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTc                                                            >  2:306968/1‑150 (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)
                                                                                                        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:630823/150‑1 (MQ=255)
                                                                                                           cGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATa                                           >  2:1395132/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)
                                                                                                                       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)
                                                                                                                                   cTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCTATATTCAGTGAAGTACt                   <  2:1423686/150‑1 (MQ=255)
                                                                                                                                   cTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACt                   >  2:369549/1‑150 (MQ=255)
                                                                                                                                   cTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACt                   >  2:1053282/1‑150 (MQ=255)
                                                                                                                                   cTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACt                   >  2:531894/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)
                                                                                                                                         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                       <  1:721655/132‑1 (MQ=255)
                                                                                                                                                 gTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATa     >  2:956663/1‑150 (MQ=255)
                                                                                                                                                   gAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAAc   >  1:1157337/1‑150 (MQ=255)
                                                                                                                                                   gAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAAc   >  1:922083/1‑150 (MQ=255)
                                                                                                                                                    cgTAATTACGCTTGGATTTCGGTTCAACATCTTCTTCGGGATGTTTCTCTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAAtttct  <  2:1550804/149‑5 (MQ=255)
                                                                                                                                                    aGTAATTACGCATGGATATCGGTTCAACATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAATGTAATTCTGTGGCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACt  <  2:1545082/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TTACTGTAGCCAATAGAGTAATTGCTACGATTTTGTCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACT  >  NC_000913/153781‑154078

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: