Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,401,582 A→T 6.3% L584Q (CTG→CAG)  csrD ← targeting factor for csrBC sRNA degradation

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,401,5820AT 6.3% 541.0 / 19.1 157L584Q (CTG→CAG) csrDtargeting factor for csrBC sRNA degradation
Reads supporting (aligned to +/- strand):  ref base A (76/70);  new base T (4/6);  total (81/76)
Fisher's exact test for biased strand distribution p-value = 5.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.42e-01

Only 100 of 164 total aligned reads displayed.

AAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTACCCGTACCCCTAAAGCATTCACTAA  >  NC_000913/3401434‑3401730
                                                                                                                                                    |                                                                                                                                                    
aaGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACcagca                                                                                                                                                      >  1:634454/1‑149 (MQ=255)
  gTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCagcagc                                                                                                                                                    <  2:259930/149‑1 (MQ=255)
   tGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCt                                                                                                                                                   <  1:1088596/149‑1 (MQ=255)
   tGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCt                                                                                                                                                   <  3:87750/149‑1 (MQ=255)
     gCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTgg                                                                                                                                                 <  3:711809/149‑1 (MQ=255)
       tgtgAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGtt                                                                                                                                               >  1:917521/1‑149 (MQ=255)
           aGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCTGCTGGTtcaaa                                                                                                                                           <  1:1015629/149‑4 (MQ=255)
           aGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCTGCTGGTtcaaa                                                                                                                                           <  4:199071/149‑4 (MQ=255)
                 cgAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGc                                                                                                                                     <  3:849551/149‑1 (MQ=255)
                  gAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCt                                                                                                                                    <  2:1378551/149‑1 (MQ=255)
                  gAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCt                                                                                                                                    >  2:258736/1‑149 (MQ=255)
                  gAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCt                                                                                                                                    >  4:510145/1‑149 (MQ=255)
                   aaaaaaTCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCtt                                                                                                                                   >  3:307482/1‑149 (MQ=255)
                   aaaaaaTCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCtt                                                                                                                                   >  4:954534/1‑149 (MQ=255)
                      aaaTCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTtctc                                                                                                                                <  1:285743/149‑1 (MQ=255)
                         tCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAAt                                                                                                                             <  1:154582/149‑1 (MQ=255)
                            ccTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGtt                                                                                                                          >  2:1332752/1‑149 (MQ=255)
                            ccTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGtt                                                                                                                          >  3:1681287/1‑149 (MQ=255)
                                       tGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGcc                                                                                                               >  4:82770/1‑149 (MQ=255)
                                        gTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGccc                                                                                                              >  4:1037783/1‑149 (MQ=255)
                                          aaCACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCg                                                                                                            >  2:1278416/1‑149 (MQ=255)
                                            cacCGCGCTGAATCAGGGTGTGCCACTCGCTTCGCGAACGCACGCAGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCAGACCAGCCCCGGa                                                                                                          <  4:426427/149‑1 (MQ=255)
                                            cacCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGa                                                                                                          <  3:1581857/149‑1 (MQ=255)
                                             acCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGAt                                                                                                         >  3:60745/1‑149 (MQ=255)
                                             acCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGAt                                                                                                         >  3:817678/1‑149 (MQ=255)
                                                  gcTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGc                                                                                                    >  2:18955/1‑149 (MQ=255)
                                                    tGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATgg                                                                                                       <  1:1388621/144‑1 (MQ=255)
                                                    tGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATgg                                                                                                       >  2:1388621/1‑144 (MQ=255)
                                                    tGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCtt                                                                                                  <  3:510145/149‑1 (MQ=255)
                                                    tGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCtt                                                                                                  <  2:1623021/149‑1 (MQ=255)
                                                      aaTCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACGAGCTGCTGGTtcaaa                                                                                                                                           >  4:1447887/1‑103 (MQ=255)
                                                       aTCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAg                                                                                               >  3:1322944/1‑149 (MQ=255)
                                                       aTCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAg                                                                                               >  1:1134168/1‑149 (MQ=255)
                                                        tCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGt                                                                                              >  1:478434/1‑149 (MQ=255)
                                                                tGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAAc                                                                                      <  2:1636799/149‑1 (MQ=255)
                                                                 gCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACa                                                                                     <  4:713526/149‑1 (MQ=255)
                                                                  ccACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGa                                                                                                                     <  4:1587947/116‑1 (MQ=255)
                                                                  ccACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACAt                                                                                    >  3:847322/1‑149 (MQ=255)
                                                                      tCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTGTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAg                                                                                <  4:1001182/149‑1 (MQ=255)
                                                                            cgcgAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCttt                                                                          <  4:1681287/149‑1 (MQ=255)
                                                                             gcgAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTg                                                                         >  4:258601/1‑149 (MQ=255)
                                                                              cgAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGa                                                                        <  1:1278416/149‑1 (MQ=255)
                                                                              cgAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGa                                                                        <  2:188973/149‑1 (MQ=255)
                                                                               gAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGAt                                                                       >  2:965403/1‑149 (MQ=255)
                                                                                     acgcCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAg                                                                                               <  1:1414854/119‑1 (MQ=255)
                                                                                     acgcCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGa                                                                        >  3:507866/1‑142 (MQ=255)
                                                                                     acgcCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACt                                                                 >  1:834346/1‑149 (MQ=255)
                                                                                      cgcCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTg                                                                <  2:1178563/149‑1 (MQ=255)
                                                                                         cGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTCCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTa                                                             <  1:394317/149‑1 (MQ=255)
                                                                                              gCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTAc                                                        >  2:1507422/1‑149 (MQ=255)
                                                                                              gCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTAc                                                        >  4:619339/1‑149 (MQ=255)
                                                                                                gTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCa                                                      <  2:964839/149‑1 (MQ=255)
                                                                                                ctAAACCAGGCTGCTGGTCCCGGAGCAGGCTGCCACCTCGCTTTGAACCAGCTTGTGGTTGTTGGTTCGCTTCTCAATGTCTCTGTCCCGCCCCGGATGGAGCTTGAGTTAATCAACATCTAGTTCTTTGATCCCTCTCGTACttgcta                                                      >  4:1326773/2‑145 (MQ=255)
                                                                                                   aaCCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCg                                                   >  2:504895/1‑149 (MQ=255)
                                                                                                     ccTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCGACTGGTACTTACCAGCGTc                                                 <  4:57257/149‑1 (MQ=255)
                                                                                                             cTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGc                                         <  3:177481/149‑1 (MQ=255)
                                                                                                             cTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGc                                         <  4:847322/149‑1 (MQ=255)
                                                                                                                gTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCtt                                                                           >  4:1660496/1‑111 (MQ=255)
                                                                                                                 tcacGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCTTTCGCTTCTCAATGTTTCTGACCAGCCCCGGACAGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCTAACTGGTCCTTACCAGCGTCAAAACAGCCTgg                                     >  4:1298196/4‑149 (MQ=255)
                                                                                                                   ccGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGtt                                   <  4:481337/149‑1 (MQ=255)
                                                                                                                   ccGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGtt                                   <  2:1134168/149‑1 (MQ=255)
                                                                                                                     ggAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGa                                 >  4:358859/1‑149 (MQ=255)
                                                                                                                     ggAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGa                                 >  3:1552718/1‑149 (MQ=255)
                                                                                                                      gAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGAc                                <  1:258736/149‑1 (MQ=255)
                                                                                                                        gCAGGCTTCCACCAGGCTTTGAACCAGCGGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACgg                              >  3:1447871/1‑149 (MQ=255)
                                                                                                                        gCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGa                                                                                                >  4:680612/1‑83 (MQ=255)
                                                                                                                        gCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACGCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAACCAGGCTGGtttgacgc                             <  3:1371558/149‑7 (MQ=255)
                                                                                                                            gCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTAc                          >  3:19568/1‑149 (MQ=255)
                                                                                                                              ttCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGa                                                                                                                     <  3:1328947/56‑1 (MQ=255)
                                                                                                                              ttCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGa                                                                                                                     >  4:1328947/1‑56 (MQ=255)
                                                                                                                              ttCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTAccc                        >  2:1155246/1‑149 (MQ=255)
                                                                                                                              ttCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTAccc                        >  4:1369596/1‑149 (MQ=255)
                                                                                                                                 caccaGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTACCCGTa                     >  4:551896/1‑149 (MQ=255)
                                                                                                                                  accaGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTACCCGTAc                    >  4:1223808/1‑149 (MQ=255)
                                                                                                                                  accaGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTACCCGTAc                    <  1:1155246/149‑1 (MQ=255)
                                                                                                                                   ccaGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTACCCGTAcc                   >  3:784642/1‑149 (MQ=255)
                                                                                                                                    caGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAAc                                             <  3:481168/122‑1 (MQ=255)
                                                                                                                                     aGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCaa                                                                   >  4:1501614/1‑99 (MQ=255)
                                                                                                                                     aGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTg                                      >  2:543304/1‑128 (MQ=255)
                                                                                                                                     aGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTg                                      <  1:543304/128‑1 (MQ=255)
                                                                                                                                     aGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTACCCGTAcccc                 >  2:498825/1‑149 (MQ=255)
                                                                                                                                     aGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTACCCGTAcccc                 >  4:433529/1‑149 (MQ=255)
                                                                                                                                      ggCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTACCCGTAccc                  >  1:1709962/1‑147 (MQ=255)
                                                                                                                                      ggCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTACCCGTAccc                  <  2:1709962/147‑1 (MQ=255)
                                                                                                                                      ggCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTACCCGTACCCCt                >  1:1502959/1‑149 (MQ=255)
                                                                                                                                         tttGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGTTGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTAAAACCAGCCTGGTTGACGGCTACCCGTACCCCTaaa             >  2:726413/1‑149 (MQ=255)
                                                                                                                                         tttGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTACCCGTACCCCTaaa             >  1:85192/1‑149 (MQ=255)
                                                                                                                                         tttGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTACCCGTACCCCTaaa             <  2:10779/149‑1 (MQ=255)
                                                                                                                                         tttGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTACCCGTACCCCTaaa             <  3:1476778/149‑1 (MQ=255)
                                                                                                                                           agAACCAGCTGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTACCCGTACCCCTaaa             >  2:1280179/2‑147 (MQ=255)
                                                                                                                                           agAACCAGCTGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTACCCGTACCCCTaaa             <  1:1280179/146‑1 (MQ=255)
                                                                                                                                            gAACCAGCAGCTGGTTCTCTGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCttt                                                                          <  3:1641932/85‑1 (MQ=255)
                                                                                                                                              aCCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCaa                                                                   <  2:1073992/90‑1 (MQ=255)
                                                                                                                                              aCCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCaa                                                                   >  1:1073992/1‑90 (MQ=255)
                                                                                                                                              aCCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACt                                                                 >  4:1291691/1‑92 (MQ=255)
                                                                                                                                              aCCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGATGGCTACCCGTACCCCTAAAGCAtt        <  2:272806/149‑1 (MQ=255)
                                                                                                                                              aCCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTACCCGTACCCCTAAAGCAtt        >  3:1149981/1‑149 (MQ=255)
                                                                                                                                               ccAGCAGGTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTACCCGTACCCCTAAAGCATTc       >  2:933398/1‑149 (MQ=255)
                                                                                                                                               ccAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTACCCGTACCCCTAAAGCATTc       >  1:271799/1‑149 (MQ=255)
                                                                                                                                                cagcagCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTACCCGTACCCCTAAAGCATTCa      >  2:70282/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTACCCGTACCCCTAAAGCATTCACTAA  >  NC_000913/3401434‑3401730

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