Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,533,713 A→C 5.5% F47V (TTT→GTT)  yddH ← flavin reductase‑like protein YddH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,533,7130AC 5.5% 362.7 / 9.4 124F47V (TTT→GTT) yddHflavin reductase‑like protein YddH
Reads supporting (aligned to +/- strand):  ref base A (61/56);  new base C (4/3);  total (65/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 = 9.98e-01

Only 100 of 130 total aligned reads displayed.

TTCATCACGCCCCGACACACTTCCCACCGCCCACGTCCAGTTAGTTGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGGTGTAATTCTATGGGGATGAATCGGCTCATGGTATCGCCT  >  NC_000913/1533564‑1533861
                                                                                                                                                     |                                                                                                                                                    
ttCATCACGCCCCGACACACTTCCCACCGCCCACGTCCAGTTAGTTGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCaaa                                                                                                                                                      >  1:271518/1‑150 (MQ=255)
ttCATCACGCCCCGACACACTTCCCACCGCCCACGTCCAGTTAGTTGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCaaa                                                                                                                                                      >  1:367560/1‑150 (MQ=255)
      aCGCCCCGACACACTTCCCACCGCCCACGTCCAGTTAGTTGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCAc                                                                                                                                                <  2:1525516/150‑1 (MQ=255)
       cGCCCCGACACACTTCCCACCGCCCACGTCCAGTTAGTTGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCAcc                                                                                                                                               >  1:1622766/1‑150 (MQ=255)
        gCCCCGACACACTTCCCACCGCCCACGTCCAGTTAGTTGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCg                                                                                                                                              >  1:2226286/1‑150 (MQ=255)
         ccccGACACACTTCCCACCGCCCACGTCCAGTTAGTTGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGTTTGAACCTCCAcgcg                                                                                                                                             <  2:931180/150‑4 (MQ=255)
         ccccGACACACTTCCCACCGCCCACGTCCAGTTAGTTGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCgg                                                                                                                                             >  1:949328/1‑150 (MQ=255)
         ccccGACACACTTCCCACCGCCCACGTCCAGTTAGTTGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCgg                                                                                                                                             <  1:1638401/150‑1 (MQ=255)
             gACACACTTCCCACCGCCCACGTCCAGTTAGTTGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGTTTGAACCTCCAcgcgtggc                                                                                                                                         <  1:1492282/150‑8 (MQ=255)
             gACACACTTCCCACCGCCCACGTCCAGTTAGTTGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCAtt                                                                                                                                         >  2:452080/1‑150 (MQ=255)
               cacacTTCCCACCGCCCACGTCCAGTTAGTTGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGa                                                                                                                                       >  2:2185925/1‑150 (MQ=255)
                 cacTTCCCACCGCCCACGTCCAGTTAGTTGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGAcc                                                                                                                                     <  1:841209/150‑1 (MQ=255)
                   cTTCCCACCGCCCACGTCCAGTTAGTTGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAg                                                                                                                                   <  2:548304/150‑1 (MQ=255)
                    ttCCCACCGCCCACGTCCAGTTAGTTGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAgg                                                                                                                                  <  1:995783/150‑1 (MQ=255)
                      cccACCGCCCACGTCCAGTTAGTTGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCt                                                                                                                                <  1:1170478/150‑1 (MQ=255)
                                aCGTCCAGTTAGTTGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATaa                                                                                                                      >  1:1216321/1‑150 (MQ=255)
                                  gTCCAGTTAGTTGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATg                                                                                                                    <  2:1779860/150‑1 (MQ=255)
                                    ccAGTTAGTTGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGtt                                                                                                                  >  2:1654740/1‑150 (MQ=255)
                                     cAGTTAGTTGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTa                                                                                                                 <  1:2179314/150‑1 (MQ=255)
                                      agttagttGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTAc                                                                                                                >  2:1955441/1‑150 (MQ=255)
                                      agttagttGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTAc                                                                                                                <  2:271518/150‑1 (MQ=255)
                                       gttagttGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCAGGTCCAAGTCGATTTATATCCCACAATCGCCACGCGTGGAGGTACAAACTCCACCGGCATTGCCCAGGCTGACGCCATAATGTTAcg                                                                                                               >  2:1311391/1‑150 (MQ=255)
                                        ttagttGCTTCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACTCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTAcgc                                                                                                              <  1:1656134/150‑1 (MQ=255)
                                         tagttGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAACCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTAcgcc                                                                                                             <  2:144091/150‑1 (MQ=255)
                                         tagttGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTAcgcc                                                                                                             <  1:735675/150‑1 (MQ=255)
                                              gctgcAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTgg                                                                                                        <  1:180015/150‑1 (MQ=255)
                                               ctgcAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTggg                                                                                                       <  1:2396901/150‑1 (MQ=255)
                                                  cAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGAtt                                                                                                    <  1:874801/150‑1 (MQ=255)
                                                    aCGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGc                                                                                                  <  1:570784/150‑1 (MQ=255)
                                                         ctggATAACTATGCCAAATTTACGTTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACTCGTGTAGTTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATc                                                                                             <  2:886967/148‑1 (MQ=255)
                                                         cGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATc                                                                                             <  1:1172480/150‑1 (MQ=255)
                                                          gggATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCCGGCTGCCGCCATAAGGTTACGCCGCTGGGATTGCTCATCg                                                                                            >  2:2314135/1‑150 (MQ=255)
                                                              tAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACAACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGaaaa                                                                                        >  2:1341794/1‑149 (MQ=255)
                                                              tAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGGCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAAc                                                                                        <  2:1363079/150‑1 (MQ=255)
                                                                aCGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTTACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTg                                                                                      >  1:724437/1‑150 (MQ=255)
                                                                     gCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGAt                                                                                 <  2:2485991/150‑1 (MQ=255)
                                                                       cAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCa                                                                               <  2:1145159/150‑1 (MQ=255)
                                                                       cAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCa                                                                               >  2:134665/1‑150 (MQ=255)
                                                                       cAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCa                                                                               >  1:1742310/1‑150 (MQ=255)
                                                                       cAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCa                                                                               <  1:1756726/150‑1 (MQ=255)
                                                                       cAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCGCATCGAAACTGGTGATCa                                                                               <  2:1036048/150‑1 (MQ=255)
                                                                         aaTTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAg                                                                             <  1:1072402/150‑1 (MQ=255)
                                                                          aTTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGa                                                                            <  2:769780/150‑1 (MQ=255)
                                                                          aTTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATAGGTCATAGAAACTGGTGATCAAGa                                                                            >  1:627379/1‑150 (MQ=255)
                                                                            ttACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGAcc                                                                          >  2:1287041/1‑150 (MQ=255)
                                                                                   tATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCAGCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACGGTTGGGc                                                                   <  2:1123645/150‑1 (MQ=255)
                                                                                      gTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTAATCAAGACCGTTGGGCCAt                                                                >  2:1901618/1‑150 (MQ=255)
                                                                                       ttCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCCCGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATg                                                               >  2:622720/1‑150 (MQ=255)
                                                                                       ttCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATg                                                               <  1:1941912/150‑1 (MQ=255)
                                                                                        tcaatcaaCTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCa                                                                 >  2:1995801/1‑147 (MQ=255)
                                                                                        tcaatcaaCTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCa                                                                 <  1:1995801/147‑1 (MQ=255)
                                                                                        tcaatcaaCTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGa                                                              >  2:2220306/1‑150 (MQ=255)
                                                                                        tcaatcaaCTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGa                                                              <  2:1965832/150‑1 (MQ=255)
                                                                                            tcaaCTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTa                                                          <  1:990460/150‑1 (MQ=255)
                                                                                               aCTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACa                                                       <  1:1149761/150‑1 (MQ=255)
                                                                                                cTCCCTGGTCCATGTCGATTTATCTACCACAATCGGCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAg                                                      >  1:248116/1‑150 (MQ=255)
                                                                                                cTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAg                                                      >  2:709833/1‑150 (MQ=255)
                                                                                                cTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAg                                                      <  2:1295610/150‑1 (MQ=255)
                                                                                                cTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAg                                                      >  1:2273632/1‑150 (MQ=255)
                                                                                                cTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAg                                                      >  2:1372561/1‑150 (MQ=255)
                                                                                                  cccTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGAc                                                    >  1:502630/1‑150 (MQ=255)
                                                                                                     tGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGAt                                                                                                     >  1:230211/1‑98 (MQ=255)
                                                                                                     tGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGAt                                                                                                     <  2:230211/98‑1 (MQ=255)
                                                                                                      ggTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCt                                                >  2:2191142/1‑150 (MQ=255)
                                                                                                         ccATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGc                                             <  2:732125/150‑1 (MQ=255)
                                                                                                           aTGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCAt                                           >  2:1100133/1‑150 (MQ=255)
                                                                                                            tGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATg                                          <  2:2114815/150‑1 (MQ=255)
                                                                                                               cGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGgtg                                       <  1:1100133/150‑1 (MQ=255)
                                                                                                               cGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGgtg                                       >  2:474860/1‑150 (MQ=255)
                                                                                                               cGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGgtg                                       <  2:1914640/150‑1 (MQ=255)
                                                                                                                gATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGgtgt                                      >  2:700038/1‑150 (MQ=255)
                                                                                                                gATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGgtgt                                      >  2:636490/1‑150 (MQ=255)
                                                                                                                 aTTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGGTGTa                                     <  2:224763/150‑1 (MQ=255)
                                                                                                                        tACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGGTGTAATTCTAt                              >  1:1119265/1‑150 (MQ=255)
                                                                                                                        tACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGGTGTAATTCTAt                              >  2:1772810/1‑150 (MQ=255)
                                                                                                                         aCCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGGTGTAATTCTATg                             >  1:458867/1‑150 (MQ=255)
                                                                                                                          ccACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGGTGTAATTCTATgg                            <  1:2220306/150‑1 (MQ=255)
                                                                                                                           cacaATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGGTGTAATTCTATggg                           <  1:363615/150‑1 (MQ=255)
                                                                                                                           cacaATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGGTGTAATTCTATggg                           <  1:2261538/150‑1 (MQ=255)
                                                                                                                           cacaATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGGTGTAATTCTATggg                           >  2:418523/1‑150 (MQ=255)
                                                                                                                            acaATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGGTGTAATTCTATgggg                          >  1:924204/1‑150 (MQ=255)
                                                                                                                               aTCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGGTGTAATTCTATGGGGATg                       <  2:1151705/150‑1 (MQ=255)
                                                                                                                                 cGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGGTGTAATTCTATGGGGATGaa                     <  2:420882/150‑1 (MQ=255)
                                                                                                                                   ccACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGGTGTAATTCTATGGGGATGAAta                   >  1:413534/1‑149 (MQ=255)
                                                                                                                                   ccACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGGTGTAATTCTATGGGGATGAATc                   <  1:2349107/150‑1 (MQ=255)
                                                                                                                                    tgccggTGGAGTTTGAACCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGGTGTAATTCTATGGGGATGAATCg                  >  1:1409641/6‑150 (MQ=255)
                                                                                                                                      ccggTGGAGTTTGAACCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGGTGTAATTCTATGGGGATGAATCGGc                >  2:2251936/4‑150 (MQ=255)
                                                                                                                                      ccggTGGAGTTTGAACCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGGTGTAATTCTATGGGGATGAATCGGc                >  2:1887478/4‑150 (MQ=255)
                                                                                                                                       cggTGGAGTTTGAACCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGGTGTAATTCTATGGGGATGAATCGGCt               >  2:1745706/3‑150 (MQ=255)
                                                                                                                                         gTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGGTGTAATTCTATGGGGATGAATCGGCTCa             <  1:2144403/150‑1 (MQ=255)
                                                                                                                                           ggAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGGTGTAATTCTATGGGGATGAATCGGCTCATg           <  2:144571/150‑1 (MQ=255)
                                                                                                                                              ggTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGGTGTAATTCTATGGGGATGAATCGGCTCATGGTa        >  2:1573181/1‑150 (MQ=255)
                                                                                                                                                ttCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGGTGTAATTCTATGGGGATGAATCGGCTCATGGTATc      >  2:386637/1‑150 (MQ=255)
                                                                                                                                                ttCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGGTGTAATTCTATGGGGATGAATCGGCTCATGGTATc      >  1:416536/1‑150 (MQ=255)
                                                                                                                                                ttCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGGTGTAATTCTATGGGGATGAATCGGCTCATGGTATc      >  1:273223/1‑150 (MQ=255)
                                                                                                                                                 tCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGGTGTAATTCTATGGGGATGAATCGGCTCATGGTATCg     >  1:1881116/1‑150 (MQ=255)
                                                                                                                                                 tCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGGTGTAATTCTATGGGGATGAATCGGCTCATGGTATCg     <  2:1325851/150‑1 (MQ=255)
                                                                                                                                                 tCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGGTGTAATTCTATGGGGATGAATCGGCTCATGGTATCg     <  1:1019678/150‑1 (MQ=255)
                                                                                                                                                   aaaCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGGTGTAATTCTATGGGGATGAATCGGCTCATGGTATCGcc   <  1:416755/150‑1 (MQ=255)
                                                                                                                                                    aaCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGGTGTAATTCTATGGGGATGAATCGGCTCATGGTATCGCCt  >  1:1283233/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TTCATCACGCCCCGACACACTTCCCACCGCCCACGTCCAGTTAGTTGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGGTGTAATTCTATGGGGATGAATCGGCTCATGGTATCGCCT  >  NC_000913/1533564‑1533861

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