Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,444,193 A→C 5.1% T28P (ACC→CCC)  tamB → translocation and assembly module for autotransporter export, inner membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,444,1930AC 5.1% 422.9 / 11.1 119T28P (ACC→CCC) tamBtranslocation and assembly module for autotransporter export, inner membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (59/53);  new base C (3/3);  total (63/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

Only 100 of 135 total aligned reads displayed.

TGCCGTACCGGTCGCGGATAAAGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATC  >  NC_000913/4444046‑4444340
                                                                                                                                                   |                                                                                                                                                   
tGCCGTACCGGTCGCGGATAAAGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGcac                                                                                                                                                    <  4:1557288/149‑1 (MQ=255)
tGCCGTACCGGTCGCGGAAAAAGACGAACACGGGGTACCGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAATAAAATCAGCCTCGGAGTGGTGATCGTTATCTTACTGTTGCCGGGATCGGTGGCGCTTCTGGTGGGcac                                                                                                                                                    <  2:1863544/149‑1 (MQ=255)
     tACCGGTCGCGGATAAAGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGTGCccacca                                                                                                                                               <  1:832570/149‑1 (MQ=255)
         ggTCGCGGATAAAGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGTGCCCACCAgaaa                                                                                                                                           <  4:682988/149‑4 (MQ=255)
           tCGCGGATAAAGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGcc                                                                                                                                         <  1:768194/149‑1 (MQ=255)
                 aTAAAGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGcacca                                                                                                                                                  <  3:1783151/134‑1 (MQ=255)
                 aTAAAGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGAGcacca                                                                                                                                                  >  4:1783151/1‑134 (MQ=255)
                      gACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGTGCCCACCAgaaa                                                                                                                                           >  3:1761862/1‑133 (MQ=255)
                      gACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTa                                                                                                                              <  2:1300488/149‑1 (MQ=255)
                         gAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGTGCCCACCAgaa                                                                                                                                            >  3:740239/1‑130 (MQ=255)
                            cacGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGTGCCCACCAgaaac                                                                                                                                          >  4:380777/1‑127 (MQ=255)
                               gggTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAgcg                                                                                                                     >  2:235585/1‑149 (MQ=255)
                                 gTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAgcggc                                                                                                                   <  3:275878/149‑1 (MQ=255)
                                      aGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATc                                                                                                              >  4:1317527/1‑149 (MQ=255)
                                       gTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAg                                                                                                                       >  1:359206/1‑139 (MQ=255)
                                       gTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAg                                                                                                                       <  2:359206/139‑1 (MQ=255)
                                         tttACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCt                                                                                                           <  2:1218678/149‑1 (MQ=255)
                                           tACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTgg                                                                                                         <  1:1854465/149‑1 (MQ=255)
                                           tACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTgg                                                                                                         <  1:235585/149‑1 (MQ=255)
                                           tACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCACGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCCgg                                                                                                         <  4:1265820/149‑1 (MQ=255)
                                              aTCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGc                                                                                                                                      >  2:1566158/1‑117 (MQ=255)
                                              aTCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGc                                                                                                                                      <  1:1566158/117‑1 (MQ=255)
                                                         gCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGa                                                                                           <  3:588693/149‑1 (MQ=255)
                                                         gCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGa                                                                                           >  3:1789432/1‑149 (MQ=255)
                                                          ccAGAATTATGAGTTTATGCAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGat                                                                                          >  3:1117608/1‑149 (MQ=255)
                                                             gAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAg                                                                                                  <  3:1435379/138‑1 (MQ=255)
                                                              aaTTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTg                                                                                      <  4:1554384/149‑1 (MQ=255)
                                                                ttATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGc                                                                                    <  4:1638422/149‑1 (MQ=255)
                                                                   tGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCTGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAg                                                                                 <  4:1789432/149‑1 (MQ=255)
                                                                       tttATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACTAGCGGCCTGCATCTgg                                                                                                                                >  1:1206061/1‑98 (MQ=255)
                                                                           tGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTgg                                                                                                         >  2:731709/1‑117 (MQ=255)
                                                                            ggAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTc                                                                              <  1:1734071/143‑1 (MQ=255)
                                                                            ggAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTc                                                                              >  2:1734071/1‑143 (MQ=255)
                                                                              aaaaaaaTCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTc                                                                              <  2:1826756/141‑1 (MQ=255)
                                                                              aaaaaaaTCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTc                                                                              >  1:1826756/1‑141 (MQ=255)
                                                                              aaaaaaaTCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACcggcgg                                                                      >  3:1495284/1‑149 (MQ=255)
                                                                               aaaaaaTCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCggcggc                                                                     >  1:657001/1‑149 (MQ=255)
                                                                                aaaaaTCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAgc                                                                                                                      >  3:324416/1‑99 (MQ=255)
                                                                                   aaTCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTg                                                                                                                                 >  3:163404/1‑85 (MQ=255)
                                                                                    aTCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGgcg                                                                >  1:1374829/1‑149 (MQ=255)
                                                                                              gCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCGCCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGt                                                                                                       >  2:486607/1‑100 (MQ=255)
                                                                                               cGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCt                                                     <  1:409612/149‑1 (MQ=255)
                                                                                                gTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCtt                                                    <  4:813177/149‑1 (MQ=255)
                                                                                                  ggTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGt                                                  <  2:1236361/149‑1 (MQ=255)
                                                                                                   gttatcgttatcTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTgg                                                                                                         <  2:1551833/93‑1 (MQ=255)
                                                                                                   gttatcgttatcTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTgg                                                                                                         >  1:1551833/1‑93 (MQ=255)
                                                                                                    ttatcgttatcTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCt                                                >  3:1116540/1‑149 (MQ=255)
                                                                                                      atcgttatcTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGa                                                                                                <  2:1472123/99‑1 (MQ=255)
                                                                                                      atcgttatcTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGa                                                                                                >  1:1472123/1‑99 (MQ=255)
                                                                                                        cgttatcTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGAcg                                            >  4:1540194/1‑149 (MQ=255)
                                                                                                          ttatcTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAg                                                                                 >  1:697402/1‑110 (MQ=255)
                                                                                                          ttatcTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGAcgtt                                          >  2:1634066/1‑149 (MQ=255)
                                                                                                             tcTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGAcgttcgt                                       <  4:158292/149‑1 (MQ=255)
                                                                                                             tcTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGAcgttcgt                                       >  3:158292/1‑149 (MQ=255)
                                                                                                              cTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGAcgttcgt                                       >  4:1672538/1‑148 (MQ=255)
                                                                                                              cTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGAcgttcgt                                       <  3:1672538/148‑1 (MQ=255)
                                                                                                               ttACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAg                                                                                 >  2:790505/1‑105 (MQ=255)
                                                                                                               ttACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAg                                                                                 <  1:790505/105‑1 (MQ=255)
                                                                                                               ttACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGt                                                                               >  1:928604/1‑107 (MQ=255)
                                                                                                               ttACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATAGTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTa                                     >  2:702320/1‑149 (MQ=255)
                                                                                                                     ttGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGAcgttc                                         <  3:846500/139‑1 (MQ=255)
                                                                                                                     ttGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGAcgttc                                         >  4:846500/1‑139 (MQ=255)
                                                                                                                     ttGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGagca                               <  2:657001/149‑1 (MQ=255)
                                                                                                                         tGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCa                           <  4:1495284/149‑1 (MQ=255)
                                                                                                                          gggATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTg                                                                                      >  3:678238/1‑89 (MQ=255)
                                                                                                                          gggATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTg                                                                                      <  4:678238/89‑1 (MQ=255)
                                                                                                                          gggATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGgc                                                                 >  1:1413939/1‑110 (MQ=255)
                                                                                                                          gggATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGgc                                                                 <  2:1413939/110‑1 (MQ=255)
                                                                                                                          gggATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAg                          >  4:1819096/1‑149 (MQ=255)
                                                                                                                             aTCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGAcgttc                                         >  3:521188/1‑131 (MQ=255)
                                                                                                                             aTCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGAcgttc                                         <  4:521188/131‑1 (MQ=255)
                                                                                                                             aTCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGc                        >  4:489380/1‑148 (MQ=255)
                                                                                                                             aTCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGc                        <  3:489380/148‑1 (MQ=255)
                                                                                                                             aTCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCg                       >  1:1777088/1‑149 (MQ=255)
                                                                                                                             aTCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCg                       >  1:723349/1‑149 (MQ=255)
                                                                                                                                  tGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTg                                                                                             <  4:663151/74‑1 (MQ=255)
                                                                                                                                  tGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTAt                                    >  1:1401651/1‑131 (MQ=255)
                                                                                                                                  tGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTg                     <  1:1365270/146‑1 (MQ=255)
                                                                                                                                  tGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCg                  >  3:1507082/1‑149 (MQ=255)
                                                                                                                                    gCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCaa                                                                                  <  1:1057269/83‑1 (MQ=255)
                                                                                                                                    gCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCaa                                                                                  >  2:1057269/1‑83 (MQ=255)
                                                                                                                                      gTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAgcgg                                                                                                                    <  2:1472414/47‑1 (MQ=255)
                                                                                                                                      gTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTaa              >  3:1144468/1‑149 (MQ=255)
                                                                                                                                       tttCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTaaa             >  3:1823718/1‑149 (MQ=255)
                                                                                                                                       tttCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTaaa             >  3:1761651/1‑149 (MQ=255)
                                                                                                                                           tggtggGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTAAAAGCg         <  4:1507082/149‑1 (MQ=255)
                                                                                                                                           tggtggGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGAGATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTAAAAGCg         >  4:1671299/1‑149 (MQ=255)
                                                                                                                                            ggtggGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTAAAAGCgg        <  2:1374829/149‑1 (MQ=255)
                                                                                                                                            ggtggGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTAAAAGCgg        <  4:1116540/149‑1 (MQ=255)
                                                                                                                                              tggGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGAcgttcg                                        <  3:724722/115‑1 (MQ=255)
                                                                                                                                                ggCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCaa    <  2:1393815/149‑1 (MQ=255)
                                                                                                                                                  caccacCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTAAAAGCgg        >  3:999733/1‑143 (MQ=255)
                                                                                                                                                  caccacCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATc  >  1:1564577/1‑149 (MQ=255)
                                                                                                                                                  caccacCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATc  >  3:251919/1‑149 (MQ=255)
                                                                                                                                                   accaccAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGa                                                                                           <  1:480802/59‑1 (MQ=255)
                                                                                                                                                   accaccAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGa                                                                                           >  2:480802/1‑59 (MQ=255)
                                                                                                                                                   accaccAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGt                      <  2:103327/128‑1 (MQ=255)
                                                                                                                                                   accaccAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGt                      >  1:103327/1‑128 (MQ=255)
                                                                                                                                                   accaccAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTAAAAGCgg        <  2:139061/142‑1 (MQ=255)
                                                                                                                                                   accaccAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTAAAAGCgg        >  1:139061/1‑142 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TGCCGTACCGGTCGCGGATAAAGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATC  >  NC_000913/4444046‑4444340

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: