Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,748,853 A→G 5.9% E27G (GAA→GGA)  yfjD → UPF0053 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,748,8530AG 5.9% 365.1 / 9.5 118E27G (GAA→GGA) yfjDUPF0053 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (54/56);  new base G (4/3);  total (59/59)
Fisher's exact test for biased strand distribution p-value = 7.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.76e-01

Only 100 of 121 total aligned reads displayed.

TGGCCTACTTCGGCAGCCGAATTGTCCAGCAGTTAATCAGCTAAACCCAGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTT  >  NC_000913/2748710‑2749001
                                                                                                                                               |                                                                                                                                                    
tGGCCTACTTCGGCAGCCGAATTGTCCAGCAGTTAATCAGCTAAACCCAGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCg                                                                                                                                                 <  2:126124/149‑1 (MQ=255)
tGGCCTACTTCGGCAGCCGAATTGTCCAGCAGTTAATCAGCTAAACCCAGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCg                                                                                                                                                 <  2:762470/149‑1 (MQ=255)
   ccTACTTCGGCAGCCGAATTGTCCAGCAGTTAATCAGCTAAACCCAGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGTTTCGGACCCGGaa                                                                                                                                              <  1:525812/149‑1 (MQ=255)
   ccTACTTCGGCAGCCGAATTGTCCAGCAGTTAATCAGCTAAACCCAGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACAACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGaa                                                                                                                                              >  2:687456/1‑149 (MQ=255)
    cTACTTCGGCAGCCGAATTGTCCAGCAGTTAATCAGCTAAACCCAGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAat                                                                                                                                             <  2:709741/149‑1 (MQ=255)
        ttCGGCAGCCGAATTGTCCAGCAGTTAATCAGCTAAACCCAGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGTTTCGGACCCGGaaaaata                                                                                                                                         <  1:204075/149‑6 (MQ=255)
         tCGGCAGCCGAATTGTCCAGCAGTTAATCAGCTAAACCCAGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGTTTCGGACCCGGaaaaatag                                                                                                                                        <  1:349703/149‑7 (MQ=255)
          cGGCAGCCGAATTGTCCAGCAGTTAATCAGCTAAACCCAGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGAc                                                                                                                                       <  2:280478/149‑1 (MQ=255)
             cAGCCGAATTGTCCAGCAGTTAATCAGCTAAACCCAGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAAAGATGCCCCt                                                                                                                                    >  1:286305/1‑149 (MQ=255)
              aGCCGAATTGTCCAGCAGTTAATCAGCTAAACCCAGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTc                                                                                                                                   <  1:322797/149‑1 (MQ=255)
               gCCGAATTGTCCAGCAGTTAATCAGCTAAACCCAGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCa                                                                                                                                  >  2:1362976/1‑149 (MQ=255)
                ccGAATTGTCCAGCAGTTAATCAGCTAAACCCAGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCaa                                                                                                                                 >  1:1007386/1‑149 (MQ=255)
                    aTTGTCCAGCAGTTAATCAGCTAAACCCAGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGCCCCTCAACCGc                                                                                                                             >  2:822691/1‑149 (MQ=255)
                    aTTGTCCAGCAGTTAATCAGCTAAACCCAGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGc                                                                                                                             >  2:634320/1‑149 (MQ=255)
                     ttGTCCAGCAGTTAATCAGCTAAACCCAGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCt                                                                                                                            <  1:645177/149‑1 (MQ=255)
                       gTCCAGCAGTTAATCAGCTAAACCCAGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTAt                                                                                                                          >  1:1466668/1‑149 (MQ=255)
                        tCCAGCAGTTAATCAGCTAAACCCAGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATc                                                                                                                         >  1:1421767/1‑149 (MQ=255)
                          cagcagTTAATCAGCTAAACCCAGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGCAACCGGAATGATGACCCTCAACCGCTATCGt                                                                                                                       >  1:1403954/1‑149 (MQ=255)
                             cagTTAATCAGCTATAACCAGAAAAGGTGTTACCCCTGGTACACATTTCTACTACTACGTTGAGCTTTATTCTGCTCATCATGGTGGTAATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTg                                                                                                                    <  2:1042594/149‑1 (MQ=255)
                             cagTTAATCAGCTAAACCCAGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTg                                                                                                                    <  2:1313371/149‑1 (MQ=255)
                              agTTACTCAGCTAAACCCAGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGc                                                                                                                   <  2:1512041/149‑1 (MQ=255)
                                ttAATCAGCTAAACCCAGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGa                                                                                                                 <  2:396510/149‑1 (MQ=255)
                                ttAATCAGCTAAACCCAGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGa                                                                                                                 <  2:306250/149‑1 (MQ=255)
                                   aTCAGCTAAACCCAGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACat                                                                                                              <  1:404608/149‑1 (MQ=255)
                                      ttcTAAACCCAGAAAAGGAGTTTTCCCTGGAACACATTTCTACTACTTCGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCATTTTTTTCCGGGTCGGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATg                                                                                                           <  1:45720/147‑1 (MQ=255)
                                      atcTAAACCCAGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATg                                                                                                           <  2:813272/147‑1 (MQ=255)
                                      aGCTAAACCCAGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATg                                                                                                           <  1:1048549/149‑1 (MQ=255)
                                       gCTAAACCCAGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATgg                                                                                                          >  1:615666/1‑149 (MQ=255)
                                       gCTAAACCCAGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATgg                                                                                                          <  1:1366595/149‑1 (MQ=255)
                                       gCTAAACCCAGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATgg                                                                                                          >  1:1234155/1‑149 (MQ=255)
                                        cTAAACCCAGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGc                                                                                                         <  2:354670/149‑1 (MQ=255)
                                                aGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAgg                                                                                                 <  2:144859/149‑1 (MQ=255)
                                                aGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAgg                                                                                                 <  1:416718/149‑1 (MQ=255)
                                                 gAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAggg                                                                                                <  1:1569521/149‑1 (MQ=255)
                                                    aaGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTaa                                                                                             <  1:1399918/149‑1 (MQ=255)
                                                       gAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCg                                                                                          >  1:1193945/1‑149 (MQ=255)
                                                        aGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCCGCGACATATGGCGAAACAGGGTAATCGc                                                                                         >  1:82045/1‑149 (MQ=255)
                                                          tttCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTc                                                                                       <  1:1525033/149‑1 (MQ=255)
                                                               ccTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCa                                                                                  <  2:932876/149‑1 (MQ=255)
                                                               ccTGGAACACATTTCTACTACTACGTTGATAATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTCCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGCCATATGGCGACACAGGGTAAACGCTCGGCCa                                                                                  >  2:616852/1‑149 (MQ=255)
                                                                cTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCaa                                                                                 >  2:1274625/1‑149 (MQ=255)
                                                                      cacaTTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGt                                                                           <  1:1148885/149‑1 (MQ=255)
                                                                       acaTTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTc                                                                          >  2:414791/1‑149 (MQ=255)
                                                                       acaTTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCCGTACAAACGCATcc                                                                         >  2:133921/1‑148 (MQ=255)
                                                                         aTTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGa                                                                        >  2:1487376/1‑149 (MQ=255)
                                                                             ctactactaCGTTTATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCTGGTCAGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGaaaaa                                                                    >  2:691332/1‑149 (MQ=255)
                                                                               actactacGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAAtt                                                                  <  1:1588652/149‑1 (MQ=255)
                                                                                    tacGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATtgctgc                                                             >  1:1153951/1‑149 (MQ=255)
                                                                                         tGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTTTCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAg                                                        >  2:1031630/1‑149 (MQ=255)
                                                                                         tGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTTTCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAg                                                        >  2:1097915/1‑149 (MQ=255)
                                                                                         tGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAg                                                        <  2:1234155/149‑1 (MQ=255)
                                                                                          gATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGc                                                       <  2:853559/149‑1 (MQ=255)
                                                                                            tCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCa                                                     <  1:591774/149‑1 (MQ=255)
                                                                                               ttattCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAgcc                                                  >  2:634286/1‑147 (MQ=255)
                                                                                               ttattCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGAc                                                  >  1:1569501/1‑149 (MQ=255)
                                                                                               ttattCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGAc                                                  >  1:10442/1‑149 (MQ=255)
                                                                                               ttattCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCCGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGAc                                                  <  2:82045/149‑1 (MQ=255)
                                                                                                  ttCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGc                                               >  2:1274591/1‑149 (MQ=255)
                                                                                                    cTGATCATCATGGTTGTCATTTCAGCCTAATTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCt                                             <  1:1090864/149‑1 (MQ=255)
                                                                                                     tGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCTCTATCGTCTGCGACATATTGAGAAACAGGGTAATCGCTAGGCCAAACGCGTAGAAAAATTGCTGCGTAAGCCAGCCCGACTg                                            >  1:1024626/1‑149 (MQ=255)
                                                                                                     tGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTg                                            >  1:1018837/1‑149 (MQ=255)
                                                                                                     tGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTg                                            >  1:1410864/1‑149 (MQ=255)
                                                                                                     tGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACCACTCGAAAAATTGCTGCGTAAGCCAGACCGCCtt                                            >  1:986036/1‑148 (MQ=255)
                                                                                                       atcatcATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGAt                                          >  1:106696/1‑149 (MQ=255)
                                                                                                           tgatGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGc                                      <  2:1383608/147‑1 (MQ=255)
                                                                                                           tcatGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGTAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGc                                      <  1:687456/149‑1 (MQ=255)
                                                                                                           tcatGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGc                                      <  2:1466668/149‑1 (MQ=255)
                                                                                                              tggtggTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTg                                   >  2:929898/1‑149 (MQ=255)
                                                                                                              tggtggTCATTTCAGCCTATATTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTg                                   <  2:1223723/149‑1 (MQ=255)
                                                                                                                gtggtCATTTCAGCCTATTTTTCCGCGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGgt                                 <  1:26562/149‑1 (MQ=255)
                                                                                                                   gtCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGa                                                   <  2:598755/128‑1 (MQ=255)
                                                                                                                    tCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTa                             >  2:338166/1‑149 (MQ=255)
                                                                                                                      aattCAGCCTATTTTTCCGGGTCCGAAACCGGAATCATCACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAAAAGCTAGGACAAACGCGACGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGATAAt                           >  2:1389056/3‑149 (MQ=255)
                                                                                                                      aTTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAAt                           <  1:1213264/149‑1 (MQ=255)
                                                                                                                      aTTTCAGCCTATTTTTACGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGTGTAATCGCTCGGCCAAACGCGTCGAAATATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAAt                           <  2:262521/149‑1 (MQ=255)
                                                                                                                        ttCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCACCCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCg                         >  1:992355/1‑149 (MQ=255)
                                                                                                                           aGCCTATTTTTCCGGGTCCGAAACCGGATTGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGATACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTACTCGGCa                      >  1:1167300/1‑149 (MQ=255)
                                                                                                                            gCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCACCTGATAAGCCTGGTGTTAATCGGCaa                     >  2:201095/1‑149 (MQ=255)
                                                                                                                            gCCTATTTTTCCGCGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCaa                     <  2:903680/149‑1 (MQ=255)
                                                                                                                             ccTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGTCAAt                    >  2:1509640/1‑149 (MQ=255)
                                                                                                                              cTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATa                   <  1:634320/149‑1 (MQ=255)
                                                                                                                               tATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATaa                  >  1:1456599/1‑149 (MQ=255)
                                                                                                                               tATTTTTCCGGGTCCGAAACCGGAATGAAGACACTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAACCAGACCCCCTGATAAGCCTGGAGTTAATAGGCAATaa                  >  2:698964/1‑149 (MQ=255)
                                                                                                                               catcattCCGGTTTCGGACCCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGCCCGCCTGATAAGCCAGCTGTTAATCGGCAAGaa                  >  1:1191624/6‑149 (MQ=255)
                                                                                                                               catcattCCGGTTTCGGACCCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATaa                  >  2:1107710/6‑149 (MQ=255)
                                                                                                                               catcattCCGGTTTCGGACCCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTGAATCGGCAATaa                  >  1:1185462/6‑149 (MQ=255)
                                                                                                                                aTTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAAc                 <  1:1445084/149‑1 (MQ=255)
                                                                                                                                   tttCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGTTAATCGCTCGGCCAATCGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTg              <  2:872175/149‑1 (MQ=255)
                                                                                                                                   tttCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATacctca              >  1:735152/1‑144 (MQ=255)
                                                                                                                                   tttCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTg              <  1:1483981/149‑1 (MQ=255)
                                                                                                                                    ttCCGGTTTCGGACCCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTgg             >  2:538841/1‑149 (MQ=255)
                                                                                                                                     tCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGt            <  2:1193945/149‑1 (MQ=255)
                                                                                                                                      ccGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTc           <  2:1182561/149‑1 (MQ=255)
                                                                                                                                        gggTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAACCGGCAATAACCTGGTCaa         <  2:686912/149‑1 (MQ=255)
                                                                                                                                          gTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAACCGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAata       >  1:734981/1‑149 (MQ=255)
                                                                                                                                           tCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAatat      >  1:847/1‑149 (MQ=255)
                                                                                                                                             cGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTc    <  2:286305/149‑1 (MQ=255)
                                                                                                                                              gAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCt   >  1:1508851/1‑149 (MQ=255)
                                                                                                                                               aaaCCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGCCCGCCTGATAAGCCTGGGGTTAATCGGCATTAACCTGGTCAATATTCtt  >  1:764293/1‑149 (MQ=255)
                                                                                                                                               aaaCCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCtt  <  2:1007386/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
TGGCCTACTTCGGCAGCCGAATTGTCCAGCAGTTAATCAGCTAAACCCAGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTT  >  NC_000913/2748710‑2749001

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: