Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,318,079 A→T 7.6% intergenic (‑73/+134) rimP ← / ← metY ribosome maturation factor for 30S subunits/tRNA‑Met

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,318,0790AT 7.6% 682.7 / 34.3 198intergenic (‑73/+134)rimP/metYribosome maturation factor for 30S subunits/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base A (94/89);  new base T (5/10);  total (99/99)
Fisher's exact test for biased strand distribution p-value = 2.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

Only 100 of 199 total aligned reads displayed.

TGCCAACCAGTTCAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAGTGGTTGCGGGG  >  NC_000913/3317931‑3318223
                                                                                                                                                    |                                                                                                                                                
    aaCCAGTTCAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATc                                                                                                                                              >  2:340596/1‑149 (MQ=255)
          ttCAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATTTATCggggtt                                                                                                                                        <  1:1234496/149‑3 (MQ=255)
           tCAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATTTATCGGGGttt                                                                                                                                       <  1:1074913/149‑1 (MQ=255)
           tCAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCtt                                                                                                                                       <  1:348899/149‑1 (MQ=255)
            cAAAACCCAGGGTCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTATTGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATTTATCGGGGtttt                                                                                                                                      <  1:686887/149‑1 (MQ=255)
             aaaaCCCAGGGCCTCAACTGGCGCAGTAATCGTCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAGGGGCCTAATGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTta                                                                                                                                     >  4:1035459/1‑149 (MQ=255)
             aaaaCCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATTTATCGGGGttttt                                                                                                                                     <  4:791891/149‑2 (MQ=255)
             aaaaCCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTta                                                                                                                                     <  1:265319/149‑1 (MQ=255)
             aaaaCCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTta                                                                                                                                     >  2:62045/1‑149 (MQ=255)
              aaaCCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTtat                                                                                                                                    >  3:533618/1‑149 (MQ=255)
               aaCCCAGGGCCTCAACTGGCGCAGTCATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTtata                                                                                                                                   <  1:843928/149‑1 (MQ=255)
               aaCCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTtata                                                                                                                                   <  3:592530/149‑1 (MQ=255)
                     gggCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTg                                                                                                                             <  4:1400761/149‑1 (MQ=255)
                         cTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAAcc                                                                                                                         >  4:644932/1‑149 (MQ=255)
                         cTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAAcc                                                                                                                         >  3:28343/1‑149 (MQ=255)
                           cAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCt                                                                                                                       <  2:1380788/149‑1 (MQ=255)
                                ggCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAAc                                                                                                                  <  1:1421323/149‑1 (MQ=255)
                                   gcAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGc                                                                                                               <  1:80857/149‑1 (MQ=255)
                                   gcAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCCGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGc                                                                                                               >  2:390803/1‑149 (MQ=255)
                                     aGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATaa                                                                                                                                >  1:1551794/1‑130 (MQ=255)
                                     aGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATaa                                                                                                                                <  2:1551794/130‑1 (MQ=255)
                                      gTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAAc                                                                                                            >  2:828072/1‑149 (MQ=255)
                                      gTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAAc                                                                                                            >  3:960301/1‑149 (MQ=255)
                                         atcatcTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGc                                                                                                         >  3:1039722/1‑149 (MQ=255)
                                            atcTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGt                                                                                                      >  3:456716/1‑149 (MQ=255)
                                             tctcTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATTTATCGGGGtttt                                                                                                                                      <  3:893240/116‑1 (MQ=255)
                                             tctcTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTc                                                                                                     <  4:867947/149‑1 (MQ=255)
                                                   ttAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGc                                                                                               >  4:46975/1‑149 (MQ=255)
                                                   ttAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGc                                                                                               >  1:12173/1‑149 (MQ=255)
                                                       tttttGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGa                                                                                                 >  4:461144/1‑143 (MQ=255)
                                                         tttGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTc                                                                                         <  2:26621/149‑1 (MQ=255)
                                                         tttGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTc                                                                                         <  4:887975/149‑1 (MQ=255)
                                                              tctAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGaa                                                                                    >  2:792259/1‑149 (MQ=255)
                                                              tctAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGaa                                                                                    >  1:527637/1‑149 (MQ=255)
                                                              tctAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGaa                                                                                    <  4:1557808/149‑1 (MQ=255)
                                                              tctAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGaa                                                                                    >  1:427915/1‑149 (MQ=255)
                                                               ctAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAg                                                                                   <  2:988101/149‑1 (MQ=255)
                                                                  aTGTGGACGAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGaa                                                                                                                           >  2:891367/1‑106 (MQ=255)
                                                                  aTGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAGATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGAt                                                                                >  1:496842/1‑149 (MQ=255)
                                                                   tgtgGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGAtt                                                                               <  4:1003094/149‑1 (MQ=255)
                                                                    gtgGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGAttt                                                                              <  1:340596/149‑1 (MQ=255)
                                                                     tgGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGAtttt                                                                             <  3:1304825/149‑1 (MQ=255)
                                                                       gACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGAtttttt                                                                           >  2:1165965/1‑149 (MQ=255)
                                                                       gACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGAtttt                                                                             <  2:1674336/147‑1 (MQ=255)
                                                                        aCAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTc                                                                          >  2:1067248/1‑149 (MQ=255)
                                                                          aaGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCaa                                                                        <  3:1604650/149‑1 (MQ=255)
                                                                           aGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCaaa                                                                       <  1:1660582/149‑1 (MQ=255)
                                                                               cACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATccc                                                                   >  3:660112/1‑149 (MQ=255)
                                                                                    ccAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCa                                                                  >  3:15260/1‑145 (MQ=255)
                                                                                    ccAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTAc                                                              >  2:1651135/1‑149 (MQ=255)
                                                                                    ccAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAACACTTTCCAGAAGGATTTTTTCAAATCCCACTAc                                                              >  3:1379641/1‑149 (MQ=255)
                                                                                      aaGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGa                                                            >  4:963935/1‑149 (MQ=255)
                                                                                        gACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAg                                                          <  1:390796/149‑1 (MQ=255)
                                                                                           aTAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGt                                                                                                      <  1:1116778/102‑1 (MQ=255)
                                                                                            tAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCg                                                      >  3:1562065/1‑149 (MQ=255)
                                                                                            tAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCg                                                      >  3:1192900/1‑149 (MQ=255)
                                                                                            tAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCg                                                      <  4:1269794/149‑1 (MQ=255)
                                                                                             aaaaaaaGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGa                                                     >  4:229728/1‑149 (MQ=255)
                                                                                             aaaaaaaGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGa                                                     >  4:1645034/1‑149 (MQ=255)
                                                                                               aaaaaGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAg                                                   >  1:615691/1‑149 (MQ=255)
                                                                                                aaaaGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACg                                                             <  1:725500/138‑1 (MQ=255)
                                                                                                aaaaGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGt                                                  >  2:853574/1‑149 (MQ=255)
                                                                                                aaaaGGGCCTAAAGCCCAGTTATTCTGTAGCCAGATAACAAAAAACCCCGATAAAACGGGGTTTTATATAACTGAACCCTATAACGTCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACACCAAAGGCCGAAGt                                                  >  3:861795/1‑149 (MQ=255)
                                                                                                   aGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCtt                                               <  3:229728/149‑1 (MQ=255)
                                                                                                    gggCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCTGGCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGcc                                                       >  2:158504/1‑140 (MQ=255)
                                                                                                    gggCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGcc                                                       <  1:158504/140‑1 (MQ=255)
                                                                                                      gCCTAAAGCCCATTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTcac                                            >  4:1057800/1‑149 (MQ=255)
                                                                                                       ccTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTcaca                                           <  4:1343820/149‑1 (MQ=255)
                                                                                                        cTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAg                                          <  2:527637/149‑1 (MQ=255)
                                                                                                        cTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAg                                          <  3:719044/149‑1 (MQ=255)
                                                                                                         tAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGa                                                                                                 <  2:730588/93‑1 (MQ=255)
                                                                                                         tAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGt                                                  >  3:1021061/1‑140 (MQ=255)
                                                                                                               ccAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATAttt                                   <  2:496842/149‑1 (MQ=255)
                                                                                                                 aGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGa                                 >  3:479954/1‑149 (MQ=255)
                                                                                                                  gTTATTCTTTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACATTATATTTGaa                                <  1:569480/149‑1 (MQ=255)
                                                                                                                     aTTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAACTCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGaaaaa                             <  1:647113/149‑1 (MQ=255)
                                                                                                                     aTTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGaaaaa                             <  4:1175298/149‑1 (MQ=255)
                                                                                                                     aTTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGaaaaa                             <  2:427915/149‑1 (MQ=255)
                                                                                                                      ttCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAg                            <  4:709332/149‑1 (MQ=255)
                                                                                                                        cTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGa                          <  1:1165965/149‑1 (MQ=255)
                                                                                                                        cTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGa                          <  3:1673841/149‑1 (MQ=255)
                                                                                                                         tGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGAc                         >  2:204817/1‑149 (MQ=255)
                                                                                                                         tGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGAc                         <  4:28343/149‑1 (MQ=255)
                                                                                                                          gTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGAct                        <  2:1352286/149‑1 (MQ=255)
                                                                                                                             gTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGaa                                                    <  1:1006719/118‑1 (MQ=255)
                                                                                                                             gTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGaa                                                    >  2:1006719/1‑118 (MQ=255)
                                                                                                                             gTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTca                                             <  1:310098/125‑1 (MQ=255)
                                                                                                                             gTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTca                                             >  2:310098/1‑125 (MQ=255)
                                                                                                                              tCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTaa                    >  4:587210/1‑149 (MQ=255)
                                                                                                                              tCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTaa                    >  4:897688/1‑149 (MQ=255)
                                                                                                                                aGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTCTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAgg                  >  1:232022/1‑149 (MQ=255)
                                                                                                                                    aaCAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGaaa              <  4:54854/149‑1 (MQ=255)
                                                                                                                                       ataaagccccGATTTATCGGGGTTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAGTg           >  1:1011353/7‑149 (MQ=255)
                                                                                                                                        aaaaaCCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAGTgg          <  1:892827/149‑1 (MQ=255)
                                                                                                                                         aaagCCCCGATTTATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAGTGGt         >  1:1623014/1‑149 (MQ=255)
                                                                                                                                                cGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGaa                                                                                    >  4:685758/1‑67 (MQ=255)
                                                                                                                                                cGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGaa                                                                                    <  3:685758/67‑1 (MQ=255)
                                                                                                                                                cGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAGTGGTTGCgggg  >  1:185053/1‑149 (MQ=255)
                                                                                                                                                cGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAGTGGTTGCgggg  <  1:836682/149‑1 (MQ=255)
                                                                                                                                                 gATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAGTGGTTGCgg    <  2:908099/146‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
TGCCAACCAGTTCAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAGTGGTTGCGGGG  >  NC_000913/3317931‑3318223

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: