Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,885,298 C→A 5.3% R75S (CGT→AGT) ‡ yidC → membrane protein insertase YidC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,885,2980CA 5.3% 508.1 / 9.0 148R75S (CGT→AGT) ‡yidCmembrane protein insertase YidC
Reads supporting (aligned to +/- strand):  ref base C (62/78);  new base A (6/2);  total (68/80)
Fisher's exact test for biased strand distribution p-value = 1.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 152 total aligned reads displayed.

AGGATAAAAACCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGATAACCCGGC  >  NC_000913/3885149‑3885446
                                                                                                                                                     |                                                                                                                                                    
aGGATAAAAACCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGCCGTGCTTGATCTGACCATCAACAccc                                                                                                                                                      >  1:1048556/1‑150 (MQ=255)
 ggATAAAAACCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCg                                                                                                                                                     >  1:114413/1‑150 (MQ=255)
   aTAAAAACCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGTGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCgtg                                                                                                                                                   <  1:1975617/150‑1 (MQ=255)
    tAAAAACCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCgtgg                                                                                                                                                  >  2:2131106/1‑150 (MQ=255)
     aaaaaCCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCgtggt                                                                                                                                                 <  2:2391999/150‑1 (MQ=255)
      aaaaCCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCgtggtg                                                                                                                                                >  2:583039/1‑150 (MQ=255)
      aaaaCCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCgtggtg                                                                                                                                                >  1:792609/1‑150 (MQ=255)
          cccGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGa                                                                                                                                            >  1:1494176/1‑150 (MQ=255)
          cccGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGa                                                                                                                                            <  1:934308/150‑1 (MQ=255)
          cccGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGa                                                                                                                                            <  1:428231/150‑1 (MQ=255)
            cGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCACCACCACGGGTGTTGATg                                                                                                                                          <  1:564022/150‑1 (MQ=255)
            cGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATg                                                                                                                                          >  2:760757/1‑150 (MQ=255)
            cGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATg                                                                                                                                          <  1:19119/150‑1 (MQ=255)
             gCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGt                                                                                                                                         <  1:72648/150‑1 (MQ=255)
             gCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGt                                                                                                                                         <  2:2289095/150‑1 (MQ=255)
              cAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGtt                                                                                                                                        <  1:354969/150‑1 (MQ=255)
                aCCTCAGGCCCAACAGACCACGCAGACAACGAGCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGa                                                                                                                                      <  2:1784268/150‑1 (MQ=255)
                aCCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGa                                                                                                                                      >  1:977479/1‑150 (MQ=255)
                  cTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGc                                                                                                                                    <  2:2420110/150‑1 (MQ=255)
                  cTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGc                                                                                                                                    >  2:95324/1‑150 (MQ=255)
                   tCAGGCCTAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCa                                                                                                                                   >  2:1917912/1‑150 (MQ=255)
                   tCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCa                                                                                                                                   <  1:1509706/150‑1 (MQ=255)
                   tCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCa                                                                                                                                   >  2:2238676/1‑150 (MQ=255)
                   tCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGCCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCACCACCACGGGTGTTGATggtcagat                                                                                                                                   <  1:1380144/150‑8 (MQ=255)
                      ggCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGc                                                                                                                                >  1:1795808/1‑150 (MQ=255)
                          cAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTctg                                                                                                                            >  2:753392/1‑150 (MQ=255)
                          cAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTctg                                                                                                                            >  2:1400148/1‑150 (MQ=255)
                           aaCAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTctgc                                                                                                                           >  1:1984759/1‑150 (MQ=255)
                              aGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCtgctgc                                                                                                                        <  1:2010997/150‑1 (MQ=255)
                               taCCACGCAGACAAGGACCACCGCAGCGGGTAGGGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGATAATGATCTTTGTTATTAATGATGTGGTTTATCTGAAGATCAACACGCGTTGTGTTTATTTTGATTAAGCTCtgctgcg                                                                                                                       <  2:89470/149‑2 (MQ=255)
                               gACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGctgcc                                                                                                                       <  1:1758840/150‑1 (MQ=255)
                                aCCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGctgcct                                                                                                                      <  2:2217761/150‑1 (MQ=255)
                                aCCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGctgcct                                                                                                                      <  1:583039/150‑1 (MQ=255)
                                 ccccGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACAAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGctgcctg                                                                                                                     >  2:1771806/4‑150 (MQ=255)
                                  cACGCTGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGctgcctgc                                                                                                                    <  2:236223/150‑1 (MQ=255)
                                  cACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGctgcctgc                                                                                                                    <  1:272339/150‑1 (MQ=255)
                                  cACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGctgcctgc                                                                                                                    <  2:1661810/150‑1 (MQ=255)
                                   aCGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCt                                                                                                                   <  1:1668513/150‑1 (MQ=255)
                                    cGCAGACAACGCCCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCtt                                                                                                                  >  2:1393283/1‑150 (MQ=255)
                                    cGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCtt                                                                                                                  >  1:646825/1‑150 (MQ=255)
                                      cAGACAACGACCACCGCAGCGGGTAGCGCTGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTAc                                                                                                                <  1:168220/150‑1 (MQ=255)
                                      cAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTAc                                                                                                                >  1:124916/1‑150 (MQ=255)
                                       aGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTAcc                                                                                                               >  2:2132019/1‑150 (MQ=255)
                                          cAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGa                                                                                                            >  2:2355075/1‑150 (MQ=255)
                                             cGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAag                                                                                                         <  2:1491736/150‑1 (MQ=255)
                                                ccaccGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGc                                                                                                      <  1:2315381/150‑1 (MQ=255)
                                                 caccGCAGCTGGTAGCTCCTCCGACCTGGGCGTAGCGGCCAGTGGGCAGTGTAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATTAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCt                                                                                                     <  2:77735/150‑1 (MQ=255)
                                                  accGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTg                                                                                                    <  1:1617392/150‑1 (MQ=255)
                                                      cAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAAct                                                                                                <  2:1625131/150‑1 (MQ=255)
                                                        gCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAActct                                                                                              >  2:331249/1‑150 (MQ=255)
                                                            gTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTAccc                                                                                          <  2:2250086/150‑1 (MQ=255)
                                                                cgccgcCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGcc                                                                                      <  1:1194114/150‑1 (MQ=255)
                                                                  ccgccgACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGt                                                                                    >  2:33957/1‑150 (MQ=255)
                                                                     ccgACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTcc                                                                                 <  1:1147318/150‑1 (MQ=255)
                                                                      cgACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAAc                                                                                                 >  2:301244/1‑133 (MQ=255)
                                                                      cgACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCa                                                                                <  1:953944/150‑1 (MQ=255)
                                                                         ccAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCt                                                                             >  2:1546198/1‑150 (MQ=255)
                                                                         ccAGGGCGTACCGGCCAGTGGCCAGGGAAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCt                                                                             <  2:2161460/150‑1 (MQ=255)
                                                                            gggCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGtt                                                                          >  2:803830/1‑150 (MQ=255)
                                                                              gCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTgg                                                                        <  1:2053031/150‑1 (MQ=255)
                                                                                 tACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTgctca                                                                     >  2:1475727/1‑146 (MQ=255)
                                                                                   ccGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACt                                                                   >  2:1084625/1‑150 (MQ=255)
                                                                                      gCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGctgcctg                                                                                                                     <  1:773356/97‑1 (MQ=255)
                                                                                         aGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCg                                                             <  2:130613/150‑1 (MQ=255)
                                                                                              ccAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGtt                                                        >  1:127943/1‑150 (MQ=255)
                                                                                               cAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGttt                                                       <  1:692045/150‑1 (MQ=255)
                                                                                               cAGGCGAAAATGATCTCGGTAAAGACCGACGTGCTTGATCTGACAATCAACACCCGTGGTGGTGATGTTGATCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGttt                                                       >  2:697308/1‑150 (MQ=255)
                                                                                                 ggggAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGtttat                                                     >  1:2168901/1‑150 (MQ=255)
                                                                                                 ggggAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGtttat                                                     <  2:2026013/150‑1 (MQ=255)
                                                                                                  gggAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGtttatt                                                    <  2:1526338/150‑1 (MQ=255)
                                                                                                  gggAAACTGATCTCGGTTAAGACCGACGTGCGTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGtttatt                                                    <  1:201379/150‑1 (MQ=255)
                                                                                                    gAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGtttattta                                                  <  1:1076718/150‑1 (MQ=255)
                                                                                                     aaaCTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTttatttat                                                 >  1:269865/1‑150 (MQ=255)
                                                                                                      aaCTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATc                                                >  2:1218577/1‑150 (MQ=255)
                                                                                                        cTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAg                                              >  2:848401/1‑150 (MQ=255)
                                                                                                        cTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAg                                              <  2:1048556/150‑1 (MQ=255)
                                                                                                        cTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAg                                              >  2:151890/1‑150 (MQ=255)
                                                                                                        cTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAg                                              <  1:2092806/150‑1 (MQ=255)
                                                                                                         tGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAgg                                             >  1:709795/1‑150 (MQ=255)
                                                                                                                 gTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGc                                     <  2:269865/150‑1 (MQ=255)
                                                                                                                 gTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGc                                     >  1:191650/1‑150 (MQ=255)
                                                                                                                   tAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCgg                                   <  1:868596/150‑1 (MQ=255)
                                                                                                                     aGACCGACGTGCTTGATCTGACCATCACCACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTc                                 <  2:2422893/150‑1 (MQ=255)
                                                                                                                     aGACCGACGTGCTTGATATGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCCACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTc                                 >  2:2012837/1‑150 (MQ=255)
                                                                                                                      gACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCt                                >  2:23205/1‑150 (MQ=255)
                                                                                                                         cGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGAc                             <  1:542355/150‑1 (MQ=255)
                                                                                                                           aCGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCg                           >  1:2354082/1‑150 (MQ=255)
                                                                                                                            cGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCgg                          >  2:194697/1‑150 (MQ=255)
                                                                                                                            cGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCgg                          >  2:186099/1‑150 (MQ=255)
                                                                                                                                      ctcaacATCACCACCACGGGTGTTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGc                >  1:2218794/6‑150 (MQ=255)
                                                                                                                                      ctcaacATCACCACCACGGGTGTTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGc                >  2:10977/6‑150 (MQ=255)
                                                                                                                                      ctcaacATCACCACCACGGGTGTTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGc                >  2:1954875/6‑150 (MQ=255)
                                                                                                                                        gACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGccc              >  1:2480526/1‑150 (MQ=255)
                                                                                                                                        gACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGccc              >  1:554764/1‑150 (MQ=255)
                                                                                                                                         aacATCACCACCACGGGTGTTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCg             >  1:1438574/3‑150 (MQ=255)
                                                                                                                                            aTCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGAt          <  1:1771806/150‑1 (MQ=255)
                                                                                                                                                 cacCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGATAAccc     >  2:2375525/1‑150 (MQ=255)
                                                                                                                                                   cccGTGGAGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGATAACCCgg   >  2:2220819/1‑150 (MQ=255)
                                                                                                                                                    ccGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGATAACCCGGc  <  1:1400148/150‑1 (MQ=255)
                                                                                                                                                    ccGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGATAACCCGGc  <  2:1759942/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGGATAAAAACCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGATAACCCGGC  >  NC_000913/3885149‑3885446

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: