Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,318,080 A→T 6.5% intergenic (‑74/+133) 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,0800AT 6.5% 409.0 / 16.0 124intergenic (‑74/+133)rimP/metYribosome maturation factor for 30S subunits/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base A (70/44);  new base T (4/4);  total (74/50)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.66e-01

Only 100 of 125 total aligned reads displayed.

AACCAGTTCAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAGTGGTTGCGGGGGCCGG  >  NC_000913/3317935‑3318228
                                                                                                                                                 |                                                                                                                                                    
aaCCAGTTCAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATc                                                                                                                                                   >  3:857635/1‑149 (MQ=255)
aaCCAGTTCAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATc                                                                                                                                                   >  3:681640/1‑149 (MQ=255)
  ccAGTTCAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCgg                                                                                                                                                 >  1:250312/1‑149 (MQ=255)
  ccAGTTCAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCgg                                                                                                                                                 >  1:604754/1‑149 (MQ=255)
  ccAGTACAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCgg                                                                                                                                                 <  4:470460/149‑1 (MQ=255)
   cAGTTCAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCggg                                                                                                                                                <  4:701073/149‑1 (MQ=255)
   cAGTTCAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCggg                                                                                                                                                <  2:34636/149‑1 (MQ=255)
      ttCAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATTTATCggggtt                                                                                                                                             <  3:606695/149‑3 (MQ=255)
        cAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATTTATCGGGGtttt                                                                                                                                           <  2:469195/149‑1 (MQ=255)
        cAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCttt                                                                                                                                           >  2:455973/1‑149 (MQ=255)
          aaaCCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTtat                                                                                                                                         >  2:785048/1‑149 (MQ=255)
          aaaCCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTATGTAGTCAGATAACAAAAAACCCCGATAGATCGGGGCTTtat                                                                                                                                         <  2:914077/149‑1 (MQ=255)
           aaCCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATTTATCGGGGttttttg                                                                                                                                        <  2:756000/149‑4 (MQ=255)
           aaCCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATTTATCGGGGttttttg                                                                                                                                        >  1:756000/1‑146 (MQ=255)
           aaCCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTtatat                                                                                                                                       >  1:199374/1‑149 (MQ=255)
                aGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCttt                                                                                                                                           <  3:342544/141‑1 (MQ=255)
                aGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCttt                                                                                                                                           >  4:342544/1‑141 (MQ=255)
                 gggCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTg                                                                                                                                  >  4:563621/1‑149 (MQ=255)
                    ccTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAAc                                                                                                                               <  2:604754/149‑1 (MQ=255)
                     cTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCttt                                                                                                                                           <  1:614012/136‑1 (MQ=255)
                     cTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCttt                                                                                                                                           >  2:614012/1‑136 (MQ=255)
                     cTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAAcc                                                                                                                              >  4:477769/1‑149 (MQ=255)
                     cTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGCGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAAcc                                                                                                                              <  4:45640/149‑1 (MQ=255)
                      tCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTtatat                                                                                                                                       <  4:228599/139‑1 (MQ=255)
                      tCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTtatat                                                                                                                                       >  3:228599/1‑139 (MQ=255)
                               gcAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGc                                                                                                                    >  1:803980/1‑149 (MQ=255)
                                ctgTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCa                                                                                                                   >  1:458422/3‑149 (MQ=255)
                                cAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCa                                                                                                                   >  3:782383/1‑149 (MQ=255)
                                    aaTCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTg                                                                                                               >  3:836171/1‑149 (MQ=255)
                                     atcatcTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGc                                                                                                              >  3:971836/1‑149 (MQ=255)
                                         tctcTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTc                                                                                                          >  3:29275/1‑149 (MQ=255)
                                         tctcTGTTAATTTTTGCTCTAATATGGACAAGCCCACCCCCAAGACATAAAAAAAGGACCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATACAGCGGGGCTTTATATAACTGAACCCGATAACCGCAACTGCGGTc                                                                                                          <  4:389759/149‑1 (MQ=255)
                                             tGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAAc                                                                                                                       <  4:552514/132‑1 (MQ=255)
                                             tGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGa                                                                                                      >  4:270043/1‑149 (MQ=255)
                                             tGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGa                                                                                                      >  2:503737/1‑149 (MQ=255)
                                               ttAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAAGCGCAACTGCGGTCTGGAGc                                                                                                    >  4:717528/1‑149 (MQ=255)
                                               ttAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCt                                                                                                         <  4:590021/144‑1 (MQ=255)
                                               ttAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCt                                                                                                         >  3:590021/1‑144 (MQ=255)
                                                  aTTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACt                                                                                                 >  1:475530/1‑149 (MQ=255)
                                                   tttttGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCgg                                                                                                                                                 >  3:815966/1‑100 (MQ=255)
                                                   tttttGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCAGTCTGGAGCACtt                                                                                                >  3:921238/1‑149 (MQ=255)
                                                     tttGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTc                                                                                              >  3:943433/1‑149 (MQ=255)
                                                       tGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCa                                                                                            <  1:443486/149‑1 (MQ=255)
                                                          tctAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGaa                                                                                         >  2:803917/1‑149 (MQ=255)
                                                            tAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCg                                                                                                             >  1:278632/1‑127 (MQ=255)
                                                            tAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAgg                                                                                       >  1:149985/1‑149 (MQ=255)
                                                            tAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAgg                                                                                       >  3:235884/1‑149 (MQ=255)
                                                            tAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAgg                                                                                       >  2:484744/1‑149 (MQ=255)
                                                              aTGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGAt                                                                                     <  2:226847/149‑1 (MQ=255)
                                                              aTGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGAt                                                                                     >  4:123217/1‑149 (MQ=255)
                                                               tgtgGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGAtt                                                                                    >  1:522480/1‑149 (MQ=255)
                                                                   gACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGAtttttt                                                                                >  3:110036/1‑149 (MQ=255)
                                                                   gACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGAtttttt                                                                                >  2:894423/1‑149 (MQ=255)
                                                                    aCAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTc                                                                               >  1:1060224/1‑149 (MQ=255)
                                                                      aaGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCaa                                                                             >  1:210041/1‑149 (MQ=255)
                                                                       aGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGACCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTCTAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCaaa                                                                            <  3:717528/149‑1 (MQ=255)
                                                                        gCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAAt                                                                           >  2:892649/1‑149 (MQ=255)
                                                                        gCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAAt                                                                           >  3:490909/1‑149 (MQ=255)
                                                                            aCCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCa                                                                       <  3:293118/149‑1 (MQ=255)
                                                                            aCCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCa                                                                       <  1:329319/149‑1 (MQ=255)
                                                                            aCCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCa                                                                       >  4:799245/1‑149 (MQ=255)
                                                                            aCCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCa                                                                       >  3:887807/1‑149 (MQ=255)
                                                                            aCCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCa                                                                       <  2:340931/149‑1 (MQ=255)
                                                                               cccAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCaaa                                                                            <  3:890113/141‑1 (MQ=255)
                                                                               cccAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCaaa                                                                            >  4:890113/1‑141 (MQ=255)
                                                                               cccAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTa                                                                    <  2:250312/149‑1 (MQ=255)
                                                                                ccAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTAc                                                                   >  3:553952/1‑149 (MQ=255)
                                                                                 cAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACGAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCAACTACg                                                                  >  2:317897/1‑149 (MQ=255)
                                                                                   aGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGaa                                                                >  2:1006764/1‑149 (MQ=255)
                                                                                    gACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATACAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAg                                                               <  4:857635/149‑1 (MQ=255)
                                                                                     aCATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAgg                                                              >  3:238959/1‑149 (MQ=255)
                                                                                        tAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCa                                                                       <  2:963886/137‑1 (MQ=255)
                                                                                        tAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCa                                                                       >  1:963886/1‑137 (MQ=255)
                                                                                            aaaaGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGt                                                       >  3:308471/1‑149 (MQ=255)
                                                                                              aaGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGATGGCCGAAGTCt                                                     <  1:558467/149‑1 (MQ=255)
                                                                                                 ggCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTca                                                  <  4:112592/149‑1 (MQ=255)
                                                                                                 ggCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTca                                                  <  2:605945/149‑1 (MQ=255)
                                                                                                 ggCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCGGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTca                                                  <  3:270043/149‑1 (MQ=255)
                                                                                                           ccAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATAttt                                        >  4:858340/1‑149 (MQ=255)
                                                                                                           ccAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATAttt                                        >  1:875774/1‑149 (MQ=255)
                                                                                                           ccAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATAttt                                        >  1:859022/1‑149 (MQ=255)
                                                                                                             aGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGa                                      <  1:441122/149‑1 (MQ=255)
                                                                                                               ttATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGaaa                                    >  2:1034692/1‑149 (MQ=255)
                                                                                                                  ttCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAg                                 <  4:553952/149‑1 (MQ=255)
                                                                                                                  ttCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAg                                 <  3:821847/149‑1 (MQ=255)
                                                                                                                  ttCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGAAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAg                                 >  4:755027/1‑149 (MQ=255)
                                                                                                                   tCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAgg                                >  3:461466/1‑149 (MQ=255)
                                                                                                                    cTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGa                               <  2:220957/149‑1 (MQ=255)
                                                                                                                    cTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGa                               <  2:664526/149‑1 (MQ=255)
                                                                                                                         gTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTa                          >  3:551176/1‑149 (MQ=255)
                                                                                                                          tCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTaa                         <  4:836171/149‑1 (MQ=255)
                                                                                                                               ttatataaagCCCCGATTTATCGGGGTTTTATATAATTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGaa                    >  2:568666/7‑149 (MQ=255)
                                                                                                                                aaCAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGaaa                   <  3:755027/149‑1 (MQ=255)
                                                                                                                                   aaaaaaCCCCGGTAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAGTg                >  3:595148/1‑149 (MQ=255)
                                                                                                                                   aaaaaaCCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAGTg                >  4:242665/1‑149 (MQ=255)
                                                                                                                                   aaaaaaCCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAGTg                >  2:245128/1‑149 (MQ=255)
                                                                                                                                   aaaaaaCCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAGTg                >  2:435033/1‑149 (MQ=255)
                                                                                                                                       agccccGATTTATCGGGGTTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAGTGGTTg            >  1:976857/3‑149 (MQ=255)
                                                                                                                                               tAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAGTGGTTGCGGGGGcc    >  4:1040876/1‑149 (MQ=255)
                                                                                                                                                 aaTCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTAGAAAAAGGACTCTAAGGGAAAGTGGTTGCGGGGGCCgg  <  1:317897/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
AACCAGTTCAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGGGCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAGTGGTTGCGGGGGCCGG  >  NC_000913/3317935‑3318228

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