Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,533,707 G→T 5.3% P49T (CCT→ACT)  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,7070GT 5.3% 408.7 / 6.5 120P49T (CCT→ACT) yddHflavin reductase‑like protein YddH
Reads supporting (aligned to +/- strand):  ref base G (59/54);  new base T (4/3);  total (63/57)
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.85e-01

Only 100 of 122 total aligned reads displayed.

TATCTTCATCACGCCCCGACACACTTCCCACCGCCCACGTCCAGTTAGTTGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGGTGTAATTCTATGGGGATGAATCGGCTCATGGTA  >  NC_000913/1533560‑1533855
                                                                                                                                                   |                                                                                                                                                    
tATCTTCATCACGCCCCGACACACTTCCCACCGCCCACGTCCAGTTAGTTGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGtt                                                                                                                                                    <  1:1184802/150‑1 (MQ=255)
   cTTCATCACGCCCCGACACACTTCCCACCGCCCACGTCCAGTTAGTTGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCaa                                                                                                                                                 >  1:2114815/1‑150 (MQ=255)
   cTTCATCACGCCCCGACACACTTCCCACCGCCCACGTCCAGTTAGTTGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCaa                                                                                                                                                 >  2:2137998/1‑150 (MQ=255)
    ttCATCACGCCCCGACACACTTCCCACCGCCCACGTCCAGTTAGTTGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCaaa                                                                                                                                                >  1:367560/1‑150 (MQ=255)
        tcaCGCCCCGACACACTTCCCACCGCCCACGTCCAGTTAGTTGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTcc                                                                                                                                            <  2:659033/150‑1 (MQ=255)
           cGCCCCGACACACTTCCCACCGCCCACGTCCAGTTAGTTGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCAcc                                                                                                                                         <  1:578941/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:1711261/1‑150 (MQ=255)
                     cacTTCCCACCGCCCACGTCCAGTTAGTTGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGAcc                                                                                                                               <  1:841209/150‑1 (MQ=255)
                      acTTCCCACCGCCCACGTCCAGTTAGTTGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCACTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGAccc                                                                                                                              >  2:45369/1‑149 (MQ=255)
                        ttCCCACCGCCCACGTCCAGTTAGTTGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAgg                                                                                                                            <  1:995783/150‑1 (MQ=255)
                                    aCGTCCAGTTAGTTGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATaa                                                                                                                >  1:1216321/1‑150 (MQ=255)
                                        ccAGTTAGTTGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGtt                                                                                                            >  2:2337049/1‑150 (MQ=255)
                                        ccAGTTAGTTGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGtt                                                                                                            >  2:1654740/1‑150 (MQ=255)
                                        ccAGTTAGTTGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCAGCAGCCAAAATGtt                                                                                                            >  2:1663532/1‑150 (MQ=255)
                                         cAGTTAGTTGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTa                                                                                                           <  1:2179314/150‑1 (MQ=255)
                                          agttagttGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTAc                                                                                                          <  2:271518/150‑1 (MQ=255)
                                           gttagttGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTAcg                                                                                                         >  1:1192733/1‑150 (MQ=255)
                                           gttagttGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCAGGTCCAAGTCGATTTATATCCCACAATCGCCACGCGTGGAGGTACAAACTCCACCGGCATTGCCCAGGCTGACGCCATAATGTTAcg                                                                                                         >  2:1311391/1‑150 (MQ=255)
                                           gttagttGCTGCAACGACCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAAATCCACAGGCATTGACCAGGCTGCAGCCATAATGTTAcg                                                                                                         >  2:1319801/1‑150 (MQ=255)
                                            ttagttGCTTCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACTCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTAcgc                                                                                                        <  1:1656134/150‑1 (MQ=255)
                                             tagttGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAACCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTAcgcc                                                                                                       <  2:144091/150‑1 (MQ=255)
                                                 tgctgcAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTg                                                                                                   >  2:765829/1‑150 (MQ=255)
                                                  gctgcAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTgg                                                                                                  <  1:180015/150‑1 (MQ=255)
                                                   ctgcAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTggg                                                                                                 <  1:2396901/150‑1 (MQ=255)
                                                    tgcAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGa                                                                                                >  1:224763/1‑150 (MQ=255)
                                                    tgcAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGa                                                                                                >  2:1361189/1‑150 (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)
                                                                    aCGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTTACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTg                                                                                >  1:724437/1‑150 (MQ=255)
                                                                         gCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGAt                                                                           <  2:2485991/150‑1 (MQ=255)
                                                                           cAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCa                                                                         >  1:1742310/1‑150 (MQ=255)
                                                                           cAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCa                                                                         <  2:1145159/150‑1 (MQ=255)
                                                                           cAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCa                                                                         >  2:134665/1‑150 (MQ=255)
                                                                           cAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCa                                                                         <  1:1756726/150‑1 (MQ=255)
                                                                           cAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCGCATCGAAACTGGTGATCa                                                                         <  2:1036048/150‑1 (MQ=255)
                                                                             aaTTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTTATCGAAACTGGTGATCAAg                                                                       <  2:1468662/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)
                                                                                       tATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGc                                                             <  2:1880879/150‑1 (MQ=255)
                                                                                          gTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCAt                                                          >  1:266889/1‑150 (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)
                                                                                            tcaatcaaCTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGa                                                        <  2:1965832/150‑1 (MQ=255)
                                                                                            tcaatcaaCTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGa                                                        >  2:2220306/1‑150 (MQ=255)
                                                                                               atcaaCTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGAttt                                                     <  2:1216321/150‑1 (MQ=255)
                                                                                                tcaaCTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTa                                                    >  1:847231/1‑150 (MQ=255)
                                                                                                tcaaCTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTa                                                    <  1:990460/150‑1 (MQ=255)
                                                                                                   aCTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACa                                                 <  1:1149761/150‑1 (MQ=255)
                                                                                                    cTCCCTGGTCCATGTCGATTTATCTACCACAATCGGCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAg                                                >  1:248116/1‑150 (MQ=255)
                                                                                                    cTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAg                                                >  1:2273632/1‑150 (MQ=255)
                                                                                                    cTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAg                                                >  2:1372561/1‑150 (MQ=255)
                                                                                                    cTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAg                                                >  2:709833/1‑150 (MQ=255)
                                                                                                        cTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACgg                                            >  1:443625/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)
                                                                                                                tGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATg                                    <  2:2114815/150‑1 (MQ=255)
                                                                                                                   cGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGgtg                                 >  2:474860/1‑150 (MQ=255)
                                                                                                                   cGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGgtg                                 <  1:1100133/150‑1 (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)
                                                                                                                           cTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGCGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGGTGTAATTCTa                         >  1:1142502/1‑150 (MQ=255)
                                                                                                                            tACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGGTGTAATTCTAt                        >  1:1119265/1‑150 (MQ=255)
                                                                                                                            tACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGGTGTAATTCTAt                        >  2:1772810/1‑150 (MQ=255)
                                                                                                                            tACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGGTGTAATTCTAt                        >  2:1962019/1‑150 (MQ=255)
                                                                                                                             aCCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGGTGTAATTCTATg                       >  1:458867/1‑150 (MQ=255)
                                                                                                                              ccACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGGTGTAATTCTATgg                      <  1:2137998/150‑1 (MQ=255)
                                                                                                                               cacaATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGGTGTAATTCTATggg                     <  1:2261538/150‑1 (MQ=255)
                                                                                                                               cacaATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGGTGTAATTCTATggg                     <  1:363615/150‑1 (MQ=255)
                                                                                                                                acaATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGGTGTAATTCTATgggg                    >  1:924204/1‑150 (MQ=255)
                                                                                                                                   aTCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGGTGTAATTCTATGGGGATg                 <  2:106779/150‑1 (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: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)
                                                                                                                                                   |                                                                                                                                                    
TATCTTCATCACGCCCCGACACACTTCCCACCGCCCACGTCCAGTTAGTTGCTGCAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACGGCTGGCATGGTGTAATTCTATGGGGATGAATCGGCTCATGGTA  >  NC_000913/1533560‑1533855

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: