Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,993,946 A→C 5.6% V254G (GTT→GGT)  uvrC ← excinuclease UvrABC, endonuclease subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,993,9460AC 5.6% 259.2 / 2.4 90V254G (GTT→GGT) uvrCexcinuclease UvrABC, endonuclease subunit
Reads supporting (aligned to +/- strand):  ref base A (47/37);  new base C (4/2);  total (51/39)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACTGGCCTACGAAGGTTTCTACCACCTCGCTCAGTTCCGTACCGCCAGGCACTTTCGGGAAATAGCTGCGGCTGCCGAGCACTTTGCCCTGACGAATGAACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCAGCTTCTTCAAACTCCAGATTCTGGCTGGCAGTTTCCATACGACTAATGAGTTGCGTAAGCACCTG  >  NC_000913/1993799‑1994094
                                                                                                                                                   |                                                                                                                                                    
aaCTGGCCTACGAAGGTTTCTACCACCTCGCTCAGTTCCGTACCGCCAGGCACTTTCGGGAAATAGCTGCGGCTGCCGAGCACTTTGCCCTGACGTATGAACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAAc                                                                                                                                                     >  1:154641/1‑149 (MQ=255)
     gCCTACGAAGGTTTCTACCACCTCGCTCAGTTCCGTACCGCCAGGCACTTTCGGGAAATAGCTGCGGCTGCCGAGCACTTTGCCCTGACGAATGAACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGa                                                                                                                                                <  1:390676/149‑1 (MQ=255)
     gCCTACGAAGGTTTCTACCACCTCGCTCAGTTCCGTACCGCCAGGCACTTTCGGGAAATAGCTGCGGCTGCCGAGCACTTTGCCCTGACGAATGAACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGa                                                                                                                                                >  2:1319275/1‑149 (MQ=255)
      ccTACGAAGGTTTCTACCACCTCGCTCAGTTCCGTACCGCCAGGCACTTTCGGGAAATAGCTGCGGCTGCCGAGCACTTTGCCCTGACGAATGAACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAg                                                                                                                                               <  1:529305/149‑1 (MQ=255)
       cTACGAAGGTTTCTACCACCTCGCTCAGTTCCGTACCGCCAGGCACTTTCGGGAAATAGCTGCGGCTGCCGAGCACTTTGCCCTGACGAATGAACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAgg                                                                                                                                              >  1:275465/1‑149 (MQ=255)
          cGAAGGTTTCTACCACCTCGCTCAGTTCCGTACCGCCAGGCACTTTCGGGAAATAGCTGCGGCTGCCGAGCACTTTGCCCTGACGAATGAACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGgtcg                                                                                                                                           >  2:1432295/1‑149 (MQ=255)
           gAAGGTTTCTACCACCTCGCTCAGTTCCGTACCGCCAGGCACTTTCGGGAAATAGCTGCGGCTGCCGAGCACTTTGCCCTGACGAATGAACAATACTTGGACACAAGCCATGCCCGCATCGCACGCCACACCAATAACGTCGAGgtcgt                                                                                                                                          >  1:216070/1‑149 (MQ=255)
           gAAGGTTTCTACCACCTCGCTCAGTTCCGTACCGCCAGGCACTTTCGGGAAATAGCTGCGGCTGCCGAGCACTTTGCCCTGACGAATGAACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGgtcgt                                                                                                                                          >  2:466276/1‑149 (MQ=255)
              ggTTTCTACCACCTCGCTCAGTTCCGTACCGCCAGGCACTTTCGGGAAATAGCTGCGGCTGCCGAGCACTTTGCCCTGACGAATGAACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGc                                                                                                                                       <  2:1305770/149‑1 (MQ=255)
               gTTTCTACCACCTCGCTCAGTTCCGTACCGCCAGGCACTTTCGGGAAATAGCTGCGGCTGCCGAGCACTTTGCCCTGACGAATGAACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGcc                                                                                                                                      <  1:1350507/149‑1 (MQ=255)
                tttCTACCACCTCGCTCAGTTCCGTACCGCCAGGCACTTTCGGGAAATAGCTGCGGCTGCCGAGCACTTTGCCCTGACGAATGAACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCg                                                                                                                                     <  1:572360/149‑1 (MQ=255)
                     accaccTCGCTCAGTTCCGTACCGCCAGGCTCTTTCGGGAAATAGCTGCGGCTGCCGAGCACTTTGCCCTGACGAATGAACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTAtt                                                                                                                                >  2:636597/1‑149 (MQ=255)
                     accaccTCGCTCAGTTCCCTACCGCCAGGCACTTTCGGGAAATACCTGCGGCTGCCGAGCACTTTGCCCTGACGAATGAACAATACGTTGACACAATCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTAtt                                                                                                                                <  2:220733/149‑1 (MQ=255)
                      ccaccTCGCTCAGTTCCGTACCGCCAGGCACTTTCGGGAAATAGCTGCGGCTGCCGAGCACTTTGCCCTGACGAATGAACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATACCGGCGACGACCTCGACGTTATTg                                                                                                                               <  2:642834/149‑1 (MQ=255)
                      ccaccTCGCTCAGTTCCGTACCGCCAGGCACTTTCGGGAAATAGCTGCGGCTGCCGAGCACTTTGCCCTGACGAATGAACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTg                                                                                                                               <  2:1167535/149‑1 (MQ=255)
                       caccTCGCTCAGTTCCGTACCGCCAGGCACTTTCGGGAAATAGCTGCGGCTGCCGAGCACTTTTCCCTGACGAATGAACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAAAGTCGAGGTAGTCGCCCGTATTgg                                                                                                                              >  1:800924/1‑149 (MQ=255)
                        accTCGCTCAGTTCCGTACCGCCAGGCACTTTCGGGAAATAGCTGCGGCTGCCGAGCACTTTGCCCTGACGAATGAACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGa                                                                                                                             >  1:1011905/1‑149 (MQ=255)
                         ccTCGCTCAGTTCCGTACCGCCAGGCACTTTCGGGAAATAGCTGCGGCTGCCGAGCACTTTGCCCTGACGAATGAACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGaa                                                                                                                            <  2:987426/149‑1 (MQ=255)
                          cTCGCTCAGTTCCGTACCGCCAGGCACTTTCGGGAAATAGCTGCGGCTGCCGAGCACTTTGCCCTGACGAATGAACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGaaa                                                                                                                           <  2:1336154/149‑1 (MQ=255)
                          cTCGCTCAGTTCCGTACCGCCAGGCACTTTCGGGAAATAGCTGCGGCTGCCGAGCACTTTGCCCTGACGAATGAACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGaaa                                                                                                                           <  2:84985/149‑1 (MQ=255)
                             gCTCAGTTCCGTACCGCCAGGCACTTTCGGGAAATAGCTGCGGCTGCCGAGCACTTTGCCCTGACGAATGAACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGa                                                                                                                        >  2:485851/1‑149 (MQ=255)
                               gcAGTTCCGTACCGCCAGGCACTTTCGGGAAATAGCTGCGGCTGCCGAGCACTTTGCCCTGACGAATGAACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAAt                                                                                                                      <  2:1461135/148‑1 (MQ=255)
                                cAGTTCCGTACCGCCAGGCACTTTCGGGAAATAGCTGCGGCTGCCGAGCACTTTGCCCTGACGAATGAACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAAtt                                                                                                                     <  2:417701/149‑1 (MQ=255)
                                  gTTCCGTACCGCCAGGCACTTTCGGGAAATAGCTGCTGCTGCCGAGCACTTTGCCCTGACGAATGAACAATACGTTGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGt                                                                                                                   <  1:466276/149‑1 (MQ=255)
                                   ttCCGTACCGCCAGGCACTTTCGGGAAATAGCTGCGGCTGCCGAGCACTTTGCCCTGACGAATGAACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGtt                                                                                                                  >  2:925174/1‑149 (MQ=255)
                                       gTACCGCCAGGCACTTTCGGGAAATAGCTGCGGCTGCCGAGCACTTTGCCCTGACGAATGAACAATACGTGGCCACACGCCATGCCCGCATCGAACGCCACCCCAATAACGTCGAGGTCGTCGCCGGTATTGGACACGAATTGTTTTTc                                                                                                              >  1:1282294/1‑149 (MQ=255)
                                         aCCGCCAGGCACTTTCGGGAAATAGCTGCGGCTGCCGAGCACTTTGCCCTGACGAATGAACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCgg                                                                                                            <  2:1437793/149‑1 (MQ=255)
                                          ccGCCAGGCACTTTCGGGAAATAGCTGCGGCTGCCGAGCACTTTGCCCTGACGAATGAACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGt                                                                                                           <  2:918448/149‑1 (MQ=255)
                                           cGCCAGGCACTTTCGGGAAATAGCTGCGGCTGCCGAGCACTTTGCCCTGACGAATGAACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTg                                                                                                          >  2:202195/1‑149 (MQ=255)
                                                         gggAAATAGCAGTGGCTGCCGAGCACTTTGCCCTGACGAATGAACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGc                                                                                            >  2:775777/1‑149 (MQ=255)
                                                          ggAAATAGCTGCGGCTGCCGAGCACTTTGCCCTGACGAATGAACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGcc                                                                                           >  1:1290219/1‑149 (MQ=255)
                                                          ggAAATAGCTGCGGCTGCCGAGCACTTTGCCCTGACGAATGAACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGcc                                                                                           >  1:1347743/1‑149 (MQ=255)
                                                          ggAAATAGCTGCGGCTGCCGAGCACTATGCCCTGACGAATGCACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATACCCTCGAGGTCCTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGcc                                                                                           >  1:1065947/1‑149 (MQ=255)
                                                              aTAGCTGCGGCTGCCGAGCACTTTGCCCTGACGAATGAACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGaa                                                                                       <  1:568443/149‑1 (MQ=255)
                                                              aTAGCTGCGGCTGCCGAGCACTTTGCCCTGACGAATGAACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGaa                                                                                       >  2:1223564/1‑149 (MQ=255)
                                                              aTAGCTGCGGCTGCCGAGCACTTTGCCCTGACGAATGAACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGGCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGaa                                                                                       <  2:54821/149‑1 (MQ=255)
                                                               tAGCTGCGGCTGCCGAGCACTTTGCCCGGACGAATGAACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAAt                                                                                      <  1:1498981/149‑1 (MQ=255)
                                                                  cTGCGGCTGCCGAGCACTTTGCCCTGCCGAATGAACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTg                                                                                   >  1:477219/1‑149 (MQ=255)
                                                                  cTGCGGCTGCCGAGCACTTTGCCCTGACGAATGAACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTg                                                                                   >  2:511555/1‑149 (MQ=255)
                                                                   tGCGGCTGCCGAGCACTTTGCCCTGACGAATGAACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTgg                                                                                  <  1:412181/149‑1 (MQ=255)
                                                                       gCTGCCGAGCACTTTGCCCTGACGAATGAACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTcgc                                                                              <  2:1506993/149‑1 (MQ=255)
                                                                                   tttGCCCTGACGAATGAACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCAg                                                                  >  1:557610/1‑149 (MQ=255)
                                                                                      gCCCTGACGAATGAACAATACGTGGACACAAGCCATGCCCGCTTCGACCGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGCAACGAATTGATTTACGGTGACGCGTCGCCCCGCCGGAAGGTGGTCGCGAAGACGTGCAGctt                                                               >  1:253779/1‑149 (MQ=255)
                                                                                       cccTGACGAATGAACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCAGcttc                                                              >  2:354369/1‑149 (MQ=255)
                                                                                        ccTGACGAATGAACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCTTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCAGcttct                                                             <  2:275465/149‑1 (MQ=255)
                                                                                         cTGACGAATGAACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCAGcttctt                                                            <  2:182010/149‑1 (MQ=255)
                                                                                             cGAATGAACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCAGCTTCTTCaaa                                                        <  1:677189/149‑1 (MQ=255)
                                                                                               aaTGAACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCAGCTTCTTCAAACt                                                      >  1:989465/1‑149 (MQ=255)
                                                                                                 tGAACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCAGCTTCTTCAAACTcc                                                    <  2:1011905/149‑1 (MQ=255)
                                                                                                  gAACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCAGCTTCTTCAAACTCCa                                                   >  2:1425617/1‑149 (MQ=255)
                                                                                                  gAACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCAGCTTCTTCAAACTCCa                                                   >  1:632847/1‑149 (MQ=255)
                                                                                                   aaCAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCAGCTTCTTCAAACTCCAg                                                  >  2:1435015/1‑149 (MQ=255)
                                                                                                       aTACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCAGCTTCTTCAAACTCCAGATTc                                              <  1:768487/149‑1 (MQ=255)
                                                                                                         aCGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCAGCTTCTTCAAACTCCAGATTctg                                            <  1:775777/149‑1 (MQ=255)
                                                                                                           gTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCAGCTTCTTCAAACTCCAGATTctggc                                          >  2:1211767/1‑149 (MQ=255)
                                                                                                            tGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGCCGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCAGCTTCTTCAAACTCCAGATTctggct                                         >  1:85350/1‑149 (MQ=255)
                                                                                                            tGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCAGCTTCTTCAAACTCCAGATTctggct                                         <  1:1447577/149‑1 (MQ=255)
                                                                                                            gggACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCAGCTTCTTCAAACTCCAGATTctggct                                         <  1:787478/148‑1 (MQ=255)
                                                                                                                cacaAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTCGCGAATCCGTGCAGCTTCTTCAAACTCCAGATTCTGGCTGGCa                                     >  1:1347766/1‑149 (MQ=255)
                                                                                                                cacaAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCAGCTTCTTCAAACTCCAGATTCTGGCTGGCa                                     >  2:1592351/1‑149 (MQ=255)
                                                                                                                   aaGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCAGCTTCTTCAAACTCCAGATTCTGGCTGGCAGtt                                  >  1:1155166/1‑149 (MQ=255)
                                                                                                                   aaGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATAGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTAGCGAATACGTGCAGCTTCTTCAAACTCCAGATTCTGGCTGGCAGtt                                  <  2:800924/149‑1 (MQ=255)
                                                                                                                         tGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCAGCTTCTTCAAACTCCAGATTCTGGCTGGCAGTTTCCATa                            <  1:202195/149‑1 (MQ=255)
                                                                                                                            ccGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCCGCTTCTTCAAACTCCATATTCTGGCTGGCAGTTTCCATACGa                         >  2:1193210/1‑149 (MQ=255)
                                                                                                                            ccGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCAGCTTCTTCAAACTCCAGATTCTGGCTGGCAGTTTCCATACgc                         >  2:920872/1‑148 (MQ=255)
                                                                                                                            ccGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCAGCTTCTTCAAACTCCAGATTCTGGCTGGCAGTTTCCATACGa                         >  1:1249276/1‑149 (MQ=255)
                                                                                                                             cGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCAGCTTCTTCAAACTCCAGATTCTGGCTGGCAGTTTCCATACGAc                        <  2:364613/149‑1 (MQ=255)
                                                                                                                              gCATCGAACGCCACACCAATACCGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGAGTCTCTCCGCCTGAATTTGGTCGCGAATACGTGCAGCTTCTTCAAACTCCAGAGTCTGGCTGGAACTTTACAGACGaat                       >  2:86041/1‑147 (MQ=255)
                                                                                                                              gCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCAGCTTCTTCAAACTCCAGATTCTGGCTGGCAGTTTCCATACGACt                       >  1:667319/1‑149 (MQ=255)
                                                                                                                              gCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCAGCTTCTTCAAACTCCAGATTCTGGCTGGCAGTTTCCATACGACt                       >  1:1318244/1‑149 (MQ=255)
                                                                                                                              gCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCAGCTTCTTCAAACTCCAGATTCTGGCTGGCAGTTTCCATACGACt                       >  1:1466315/1‑149 (MQ=255)
                                                                                                                                    aaagCCACACCACTAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCAGCTTCTTCAAACTCCAGATTCTGGCTGGCAGTTTCCATACGACTAATGAg                 <  2:1050502/146‑1 (MQ=255)
                                                                                                                                    aaCGCCACACCACTACCGTCGAGGTCGGCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCAGCTTCTTCAAACTCCAGATTCTGGCTGGCAGTTTCCATACGACTAATGAg                 <  2:1055854/149‑1 (MQ=255)
                                                                                                                                     aCGCCACACCAATAACGTCGAGGTCGGCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCAGCTTCTTCAAACTCCAGATTCTGGCTGGCAGTTTCCATACGACTAATGAGt                <  1:449073/149‑1 (MQ=255)
                                                                                                                                      cGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCAGCTTCTTCAAACTCCAGATTCTGGCTGGCAGTTTCCATACGACTAATGAGtt               <  1:925174/149‑1 (MQ=255)
                                                                                                                                      cGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGCGACGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCAGCTTCTTCAAACTCCAGATTCTGGCTGGCAGTTTCCATACGACTAATGAGtt               >  1:1565031/1‑149 (MQ=255)
                                                                                                                                        ccACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCAGCTTCTTCAAACTCCAGATTCTGGCTGGCAGTTTCCATACGACTAATGAGTTGc             <  2:632847/149‑1 (MQ=255)
                                                                                                                                         gtttccAATACCGGCGACGACCTCGACGTTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCAGCTTCTTCAAACTCCAGATTCTGGCTGGCAGTTTCCATACGACTAATGAGTTGCg            >  1:253867/5‑149 (MQ=38)
                                                                                                                                          acacCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCAGCTTCTTCAAACTCCAGATTCTGGCTGGCAGTTTCCATACGACTAATGAGTTGCGt           >  1:1347357/1‑149 (MQ=255)
                                                                                                                                          acacCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCAGCTTCTTCAAACTCCAGATTCTGGCTGGCAGTTTCCATACGACTAATGAGTTGCGt           >  1:1528448/1‑149 (MQ=255)
                                                                                                                                          acacCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCAGCTTCTTCAAACTCCAGATTCTGGCTGGCAGTTTCCATACGACTAATGAGTTGCGt           >  2:1556131/1‑149 (MQ=255)
                                                                                                                                            acCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCAGCTTCTTCAAACTCCAGATTCTGGCTGGCAGTTTCCATACGACTAATGAGTTGCGTaa         >  2:1589593/1‑149 (MQ=255)
                                                                                                                                            acCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCAGCTTCTTCAAACTCCAGATTCTGGCTGGCAGTTTCCATACGACTAATGAGTTGCGTaa         >  1:1186442/1‑149 (MQ=255)
                                                                                                                                              cAATACCGGCGACGACCTCGACGATATTGGAAACGAATTGTTTTTCTGTGACGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCAGCTTCTTCACACTCCAGATTCTGGCTGGCAGTTTCCATACGACTAATGAGTTGCGTACGc       >  2:1507806/1‑149 (MQ=255)
                                                                                                                                              cAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCAGCTTCTTCAAACTCCAGATTCTGGCTGGCAGTTTCCATACGACTAATGAGTTGCGTAAGc       >  1:346900/1‑149 (MQ=255)
                                                                                                                                              cAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCAGCTTCTTCAAACTCCAGATTCTGGCTGGCAGTTTCCATACGACTAATGAGTTGCGTAAGc       >  2:1393468/1‑149 (MQ=255)
                                                                                                                                               aaTAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCAGCTTCTTCAAACTCCAGATTCTGGCTGGCAGTTTCCATACGACTAATGAGTTGCGTAAGCa      >  1:828258/1‑149 (MQ=255)
                                                                                                                                                aTAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCAGCTTCTTCAAACTCCAGATTCTGGCTGGCAGTTTCCATACGACTAATGAGTTGCGTAAGCAc     >  1:1383087/1‑149 (MQ=255)
                                                                                                                                                 tAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCAGCTTCTTCAAACTCCAGATTCTGGCTGGCAGTTTCCATACGACTAATGAGTTGCGTAAGCAcc    <  1:1223564/149‑1 (MQ=255)
                                                                                                                                                 tAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCAGCTTCTTCAAACTCCAGATTCTGGCTGGCAGTTTCCATACGACTAATGAGTTGCGTAAGCAcc    <  2:154641/149‑1 (MQ=255)
                                                                                                                                                   aCGTCGAGGTCGTGGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCAGCTTCTTCAAACTCCAGATTCTGGCTGGCAGTTTCCATACGACTAATGAGTTTCGTAAGCACCTg  <  1:1432295/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AACTGGCCTACGAAGGTTTCTACCACCTCGCTCAGTTCCGTACCGCCAGGCACTTTCGGGAAATAGCTGCGGCTGCCGAGCACTTTGCCCTGACGAATGAACAATACGTGGACACAAGCCATGCCCGCATCGAACGCCACACCAATAACGTCGAGGTCGTCGCCGGTATTGGAAACGAATTGTTTTTCGGTGACGCGTCGCACCGCCTGAATTTGGTCGCGAATACGTGCAGCTTCTTCAAACTCCAGATTCTGGCTGGCAGTTTCCATACGACTAATGAGTTGCGTAAGCACCTG  >  NC_000913/1993799‑1994094

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