Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,444,193 A→C 7.6% T28P (ACC→CCC)  tamB → translocation and assembly module subunit TamB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,444,1930AC 7.6% 413.4 / 22.2 143T28P (ACC→CCC) tamBtranslocation and assembly module subunit TamB
Reads supporting (aligned to +/- strand):  ref base A (65/67);  new base C (3/8);  total (68/75)
Fisher's exact test for biased strand distribution p-value = 2.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

Only 100 of 155 total aligned reads displayed.

TTTGCCGTACCGGTCGCGGATAAAGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATCTG  >  NC_000913/4444044‑4444342
                                                                                                                                                     |                                                                                                                                                     
tttGCCGTACGGGTCGCGGATAAAGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGca                                                                                                                                                       <  1:2374626/150‑1 (MQ=255)
  tGCCGTACCGGGCGCGGATAAAGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGGTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGGTATCTTACTGTTGCGGGGATCGGTGGCGTTTCTGGTGGGcacc                                                                                                                                                     <  2:209832/150‑1 (MQ=255)
   gCCGTACCGGTCGCGGATAAAGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGTAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGcacca                                                                                                                                                    <  2:1288265/150‑1 (MQ=255)
   gCCGTACCGGTCGCGGATAAAGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGcacca                                                                                                                                                    <  1:526319/150‑1 (MQ=255)
   gCCGTACCGGTCGCGGATAAAGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGcacca                                                                                                                                                    >  1:1240772/1‑150 (MQ=255)
    ccGTACCGGTCGCGGATAAAGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGcaccac                                                                                                                                                   >  1:2288349/1‑150 (MQ=255)
    ccGTACCGGTCGCGGATAAAGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGcaccac                                                                                                                                                   >  1:495314/1‑150 (MQ=255)
    ccGTACCGGTCGCGGATAAAGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGcaccac                                                                                                                                                   <  2:1588667/150‑1 (MQ=255)
      gTACCGGTCGCGGATAAAGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCAccacca                                                                                                                                                 >  2:1178559/1‑150 (MQ=255)
           ggTCGCGGATAAAGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGTGCCCACCAgaaac                                                                                                                                            <  2:1142568/150‑5 (MQ=255)
           ggTCGCGGATAAAGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGTGCCCACCAgaaac                                                                                                                                            <  2:2407477/150‑5 (MQ=255)
           ggTCGCGGATAAAGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGTGCCCACCAgaaac                                                                                                                                            <  1:1107854/150‑5 (MQ=255)
           ggTCGCGGATAAAGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGGCCCACCAgaaac                                                                                                                                            <  1:8039/150‑5 (MQ=255)
            gTCGCGGATAAAGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTGTATGGAAAAAAATCAGCCTCGCCGTGGTTTTCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGcc                                                                                                                                           >  2:1557224/1‑150 (MQ=255)
            gTCGCGGATAAAGACGAACACGGGTTACAGTTTTACATCGGGCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGcc                                                                                                                                           <  1:1003530/150‑1 (MQ=255)
                cgGATAAAGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCa                                                                                                                                       >  1:336707/1‑150 (MQ=255)
                      aaGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGt                                                                                                                                 >  1:402881/1‑150 (MQ=255)
                      aaGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGCGCACCACCAGCGGCCTGCATCTGGt                                                                                                                                 >  2:842411/1‑150 (MQ=255)
                         acaaaCACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGTGCCCACCAga                                                                                                                                               <  2:617459/130‑2 (MQ=255)
                          cGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTAttt                                                                                                                             >  2:1174559/1‑150 (MQ=255)
                           gAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTa                                                                                                                            <  2:2360633/150‑1 (MQ=255)
                           gAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTa                                                                                                                            <  1:860250/150‑1 (MQ=255)
                            aaaacGGGTTACAGTTTTACATCGGTCTGTGGCAAGAATTATGAGTTTATGGAAAAAAATAAGCATCGGCGTGTTTTTCGTTATCTTACTGTTTATGGTATCGGTGGCGTTTATGTTGGGAACAACCATCGGCCTGCATATCGTATTTaa                                                                                                                           >  2:2433021/4‑150 (MQ=255)
                             acacGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTaaa                                                                                                                          <  1:878684/150‑1 (MQ=255)
                              cacGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAg                                                                                                                         >  1:836722/1‑150 (MQ=255)
                               acGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGGTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAgc                                                                                                                        <  2:721283/150‑1 (MQ=255)
                                   gTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAgcggcg                                                                                                                    <  1:362977/150‑1 (MQ=255)
                                    ttACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGcggcgg                                                                                                                   >  2:2182714/1‑150 (MQ=255)
                                             tACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTggg                                                                                                          <  1:2187037/150‑1 (MQ=255)
                                             tACATCGGTCTGGGCCGAGAATTATGAGTTTATTGAAAAAAATCAGCCTCGGCGTGGTTTTCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGGATCTGGTATTTTAAGCGGCGGATCGCTggg                                                                                                          <  2:657302/150‑1 (MQ=255)
                                              aCATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGt                                                                                                         >  1:568048/1‑150 (MQ=255)
                                               cATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTg                                                                                                        <  1:1752334/150‑1 (MQ=255)
                                               cATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTg                                                                                                        >  2:1643896/1‑150 (MQ=255)
                                                  cGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCa                                                                                                     <  1:653068/150‑1 (MQ=255)
                                                   ggTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAg                                                                                                    >  2:2136907/1‑150 (MQ=255)
                                                          ggCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGa                                                                                             >  2:2016478/1‑150 (MQ=255)
                                                          ggCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGa                                                                                             >  1:1328784/1‑150 (MQ=255)
                                                           gCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGat                                                                                            <  2:1363870/150‑1 (MQ=255)
                                                           gCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGat                                                                                            >  2:2337518/1‑150 (MQ=255)
                                                           gCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGCCGTGGTTATCGTTATCTTACTGTAGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGCTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGat                                                                                            >  2:1908595/1‑150 (MQ=255)
                                                             cAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGatat                                                                                          <  2:1534858/150‑1 (MQ=255)
                                                               gAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTg                                                                                        <  2:501342/150‑1 (MQ=255)
                                                                aaTTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTgg                                                                                       <  2:1240772/150‑1 (MQ=255)
                                                                 aTTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGc                                                                                      >  2:144228/1‑150 (MQ=255)
                                                                  ttATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCa                                                                                     <  2:2071712/150‑1 (MQ=255)
                                                                     tGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAgg                                                                                  <  1:1892065/150‑1 (MQ=255)
                                                                        gTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGTTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCa                                                                               >  2:776429/1‑150 (MQ=255)
                                                                          ttATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCAcc                                                                             <  1:375772/150‑1 (MQ=255)
                                                                           tATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACcg                                                                            >  2:920787/1‑150 (MQ=255)
                                                                           tATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACcg                                                                            >  1:702688/1‑150 (MQ=255)
                                                                           tATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACcg                                                                            >  1:2061070/1‑150 (MQ=255)
                                                                            aTGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACcgg                                                                           <  1:1230807/150‑1 (MQ=255)
                                                                                aaaaaaaTCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCggcggc                                                                       >  1:217752/1‑150 (MQ=255)
                                                                                   aaaaTCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGTGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTgg                                                                    <  2:1239409/150‑1 (MQ=255)
                                                                                    aaaTCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGgc                                                                   >  1:116538/1‑150 (MQ=255)
                                                                                     aaTCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGAGCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGgcg                                                                  <  1:317253/150‑1 (MQ=255)
                                                                                     aaTCAGCCTCGGCGGGGGTAGCGTGAGCTTACTGTTGCGGGGATCGGTGGCGTTTCTGGTGGTCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGgcg                                                                  <  2:2398093/150‑1 (MQ=255)
                                                                                       tCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGcgcg                                                                <  1:1160735/150‑1 (MQ=255)
                                                                                        cAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGa                                                               <  1:1299619/150‑1 (MQ=255)
                                                                                          gCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGAtc                                                             >  1:1568946/1‑150 (MQ=255)
                                                                                          gCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGAtc                                                             <  1:24034/150‑1 (MQ=255)
                                                                                          gCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGAtc                                                             >  2:1828196/1‑150 (MQ=255)
                                                                                           ccTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGAtct                                                            <  1:425350/150‑1 (MQ=255)
                                                                                             tCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCa                                                          >  1:2308848/1‑150 (MQ=255)
                                                                                             tCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCa                                                          <  1:1178559/150‑1 (MQ=255)
                                                                                              cGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCAc                                                         <  2:370807/150‑1 (MQ=255)
                                                                                                gCGTGGTTATCGTTATCTTATTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCGGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCt                                                       <  2:1880142/150‑1 (MQ=255)
                                                                                                    ggTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTc                                                   >  2:2055378/1‑150 (MQ=255)
                                                                                                    ggTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATATCACCTTGTc                                                   <  1:2441255/150‑1 (MQ=255)
                                                                                                     gttatcgttatcTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCt                                                  <  1:160817/150‑1 (MQ=255)
                                                                                                       tatcgttatcTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGa                                                <  2:2058483/150‑1 (MQ=255)
                                                                                                       tatcgttatcTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGa                                                <  1:2354235/150‑1 (MQ=255)
                                                                                                         tcgttatcTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGAcg                                              >  1:919587/1‑150 (MQ=255)
                                                                                                          cgttatcTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGAcgt                                             >  1:461489/1‑150 (MQ=255)
                                                                                                              atcTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACGAGCGGCCTGCAGCTGGTATTTAAAGTTGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCAGCTGGCGCGATCTCACCTTGTCTGAcgttcgt                                         <  2:1112759/150‑1 (MQ=255)
                                                                                                               tcTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGAcgttcgtt                                        >  1:1997982/1‑150 (MQ=255)
                                                                                                                cTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTa                                       >  1:618115/1‑150 (MQ=255)
                                                                                                                    cTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGCTCTCACCTTGTCTGACGTTCGTTATGag                                   >  2:359007/1‑150 (MQ=255)
                                                                                                                      gTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGagca                                 >  1:2145061/1‑150 (MQ=255)
                                                                                                                      gTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTAGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGagca                                 >  2:2032960/1‑150 (MQ=255)
                                                                                                                        tGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGagcagc                               >  1:549833/1‑150 (MQ=255)
                                                                                                                         gCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGcc                              >  1:1422176/1‑150 (MQ=255)
                                                                                                                          cTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTGAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCa                             >  1:1476736/1‑150 (MQ=255)
                                                                                                                          cTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGATGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCa                             >  1:1449285/1‑150 (MQ=255)
                                                                                                                            gggATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAgg                           >  2:2429246/1‑150 (MQ=255)
                                                                                                                              gATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCg                         >  2:2389368/1‑150 (MQ=255)
                                                                                                                                 cGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTgg                      <  2:655768/150‑1 (MQ=255)
                                                                                                                                    gggCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCgg                   <  1:121565/149‑1 (MQ=255)
                                                                                                                                    gggCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCgg                   <  2:1943569/149‑1 (MQ=255)
                                                                                                                                      aggccgcTGGTGGTGCCCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTa                 >  1:676552/7‑150 (MQ=255)
                                                                                                                                       cGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGTCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCATCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTaa                >  1:2089202/1‑150 (MQ=255)
                                                                                                                                        gTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTaaa               <  2:85766/150‑1 (MQ=255)
                                                                                                                                         ccgcTGGTGGTGCCCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTaaaa              >  2:2299588/4‑150 (MQ=255)
                                                                                                                                            cTGGTGGGCACCACCAGGGGGCTGCATCTGTTATTTAAAGCGGCGGATCGCTGTGTGCGAGGAGTGGATATTGGCAAGGTGACGGGCGGCTGGCGCGATCTCACGTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTAAAAGCg           <  2:451784/150‑1 (MQ=255)
                                                                                                                                            cTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTAAAAGCg           >  2:2152096/1‑150 (MQ=255)
                                                                                                                                              ggtggGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTAAAAGCggg         >  1:1778701/1‑150 (MQ=255)
                                                                                                                                                 gggCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCaa      >  2:2116549/1‑150 (MQ=255)
                                                                                                                                                    caccacCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATCt   <  1:2055378/150‑1 (MQ=255)
                                                                                                                                                     accaccAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATCTg  <  1:2429246/150‑1 (MQ=255)
                                                                                                                                                     accaccAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATCTg  <  1:2136907/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTTGCCGTACCGGTCGCGGATAAAGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATCTG  >  NC_000913/4444044‑4444342

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: