Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,748,855 A→C 5.4% T28P (ACC→CCC)  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,8550AC 5.4% 329.1 / 7.1 118T28P (ACC→CCC) yfjDUPF0053 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (55/56);  new base C (4/3);  total (59/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.53e-01

Only 100 of 122 total aligned reads displayed.

TGGCCTACTTCGGCAGCCGAATTGTCCAGCAGTTAATCAGCTAAACCCAGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGC  >  NC_000913/2748710‑2749003
                                                                                                                                                 |                                                                                                                                                    
tGGCCTACTTCGGCAGCCGAATTGTCCAGCAGTTAATCAGCTAAACCCAGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCg                                                                                                                                                   <  2:762470/149‑1 (MQ=255)
tGGCCTACTTCGGCAGCCGAATTGTCCAGCAGTTAATCAGCTAAACCCAGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCg                                                                                                                                                   <  2:126124/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)
                ccGAATTGTCCAGCAGTTAATCAGCTAAACCCAGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCaa                                                                                                                                   >  1:1007386/1‑149 (MQ=255)
                    aTTGTCCAGCAGTTAATCAGCTAAACCCAGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGCCCCTCAACCGc                                                                                                                               >  2:822691/1‑149 (MQ=255)
                    aTTGTCCAGCAGTTAATCAGCTAAACCCAGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGc                                                                                                                               <  2:424020/149‑1 (MQ=255)
                    aTTGTCCAGCAGTTAATCAGCTAAACCCAGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGc                                                                                                                               >  2:634320/1‑149 (MQ=255)
                       gTCCAGCAGTTAATCAGCTAAACCCAGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTAt                                                                                                                            >  1:1466668/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:306250/149‑1 (MQ=255)
                                ttAATCAGCTAAACCCAGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGa                                                                                                                   <  2:396510/149‑1 (MQ=255)
                                ttAATCAGCTAAACCCAGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGa                                                                                                                   <  1:1494876/149‑1 (MQ=255)
                                ttAATCAGCTAAACCCAGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGa                                                                                                                   <  2:1283425/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                                                                                                             >  2:193585/1‑149 (MQ=255)
                                      aGCTAAACCCAGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATg                                                                                                             <  1:1048549/149‑1 (MQ=255)
                                       gCTAAACCCAGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATgg                                                                                                            >  1:615666/1‑149 (MQ=255)
                                       gCTAAACCCAGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATgg                                                                                                            <  1:1366595/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)
                                                         gTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCt                                                                                          >  1:932876/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)
                                                                   gAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAAcg                                                                                >  1:385219/1‑149 (MQ=255)
                                                                    aaCACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAAcgc                                                                               >  2:335213/1‑149 (MQ=255)
                                                                      cacaTTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAAcgcgg                                                                             >  2:679898/1‑148 (MQ=255)
                                                                      cacaTTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGt                                                                             <  1:1148885/149‑1 (MQ=255)
                                                                       acaTTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTc                                                                            >  2:414791/1‑149 (MQ=255)
                                                                       acaTTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCCGTACAAACGCATcc                                                                           >  2:133921/1‑148 (MQ=255)
                                                                             ctactactaCGTTTATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCTGGTCAGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGaaaaa                                                                      >  2:691332/1‑149 (MQ=255)
                                                                                 tactacGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAAc                                                                                 >  1:1434770/1‑134 (MQ=255)
                                                                                 tactacGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAAc                                                                                 <  2:1434770/134‑1 (MQ=255)
                                                                                  actacGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATtgct                                                                 <  1:1362976/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)
                                                                                                    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)
                                                                                                           tcatGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGc                                        >  2:1232872/1‑149 (MQ=255)
                                                                                                              tggtggTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTg                                     >  2:929898/1‑149 (MQ=255)
                                                                                                              tggtggTCATTTCAGCCTATATTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTg                                     <  2:1223723/149‑1 (MQ=255)
                                                                                                                gtggtCATTTCAGCCTATTTTTCCGCGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGgt                                   <  1:26562/149‑1 (MQ=255)
                                                                                                                   gtCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGa                                                     >  1:598755/1‑128 (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)
                                                                                                                        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                   <  2:106696/149‑1 (MQ=255)
                                                                                                                                aTTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAAc                   <  1:1445084/149‑1 (MQ=255)
                                                                                                                                   tttCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGTTAATCGCTCGGCCAATCGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTg                <  2:872175/149‑1 (MQ=255)
                                                                                                                                   tttCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATacctca                >  1:735152/1‑144 (MQ=255)
                                                                                                                                     tCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGt              <  1:193585/149‑1 (MQ=255)
                                                                                                                                     tCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGt              <  2:1193945/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)
                                                                                                                                               aaaCCGGAATTATGACCCTCAACCGCTATCGTCTTGGAAATATGGCGAAACAGGGTAATCTCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCGAGACCGGCTGATATGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCtt    <  1:698964/149‑1 (MQ=255)
                                                                                                                                               aaaCCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGCCCGCCTGATAAGCCTGGGGTTAATCGGCATTAACCTGGTCAATATTCtt    >  1:764293/1‑149 (MQ=255)
                                                                                                                                                 aCCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGc  <  2:385219/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TGGCCTACTTCGGCAGCCGAATTGTCCAGCAGTTAATCAGCTAAACCCAGAAAAGGAGTTTCCCCTGGAACACATTTCTACTACTACGTTGATCATTATTCTGATCATCATGGTGGTCATTTCAGCCTATTTTTCCGGGTCCGAAACCGGAATGATGACCCTCAACCGCTATCGTCTGCGACATATGGCGAAACAGGGTAATCGCTCGGCCAAACGCGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGC  >  NC_000913/2748710‑2749003

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: