Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,853,603 A→G 6.2% A17A (GCA→GCG tisB → toxic membrane persister formation peptide, LexA‑regulated

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,853,6030AG 6.2% 599.0 / 20.7 186A17A (GCA→GCGtisBtoxic membrane persister formation peptide, LexA‑regulated
Reads supporting (aligned to +/- strand):  ref base A (82/91);  new base G (5/7);  total (88/98)
Fisher's exact test for biased strand distribution p-value = 7.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.81e-01

Only 100 of 194 total aligned reads displayed.

TCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCACTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAGCACAACGTTTCTCCTTGAGATACCGCGTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTAAACAATTCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGG  >  NC_000913/3853332‑3853873
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                              
tCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCACTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAGCACAACGTTTCTCCTTGAGATACCGCGTGCACAACAGCTGGCAACAGCCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTTGGTATCGCCATTCTTATCCCCAAACTCATTGTTCCAGCACt                                                                                                                                                                                                                                                                              >  4:469298/1‑274 (MQ=255)
      tACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCACTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAGCACAACGTTTCCCCTTGAGATACCGCGTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCAAAGGAGAAGCGTATGAACCTGGTGGATAACGCCATTCTTATCCTCAAACTCATTGTTGCAGACCTTCAACTg                                                                                                                                                                                                                                                                       >  4:434242/1‑275 (MQ=255)
             tgtgTTATGCGGGGGCCGCATCGTTACCCGGCGCACTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAGCACAACGTTTCTCCTTGAGATACCGCGTGCACAACAGCTGGCAACCGGCAGCGGCAAGGTACGTCCGCTGGCCATTCCCCTGAACCACAGGCACCGCATATGACCCTGGGTGATTATCCCATTCTTACCCTCAAACTCATTGTTGAAACACCGCCACTTCTTggt                                                                                                                                                                                                                                                                 >  4:6865/1‑272 (MQ=255)
             tgtgTTATGCGGGGGCCGCATCGTTACCCGGCGCACTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAGCACAACGTTTCTCCTTGAGATACCGCGTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATg                                                                                                                                                                                                                                                                <  4:334687/275‑1 (MQ=255)
             gtggtgaTGCGGGGGCCGCATCGGTACCCGGCGCACTAAGTCCTGTGTGAACCGGGTGGGGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAGCACAACGTTTCTCCTTGAGATACCGCGTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATg                                                                                                                                                                                                                                                                <  4:485872/269‑1 (MQ=255)
                         gggCCGCATCGTTACCCGGCGCACTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAGCACAACGTTTCTCCTTGAGATACCGCGTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAAt                                                                                                                                                                                                                                                    >  2:449739/1‑275 (MQ=255)
                                            cgcACTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAGCACAACGTTTCTCCTTGAGATACCGCGTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAATTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGAtt                                                                                                                                                                                                                                 >  1:123172/1‑275 (MQ=255)
                                                           cTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAGCACAACGTTTCTCCTTGAGATACCGCGTGCACAACCGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCAATTCTTATCCCAAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTTAAATACCTGAAGTAATTCACATTCAAGTCGTCCCaaa                                                                                                                                                                                                                   >  4:530958/1‑274 (MQ=255)
                                                           cTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAGCACAACGTTTCTCCTTGAGATACCGCGTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTc                                                                                                                                                                                                                                >  4:173913/1‑261 (MQ=255)
                                                                          ccGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAGCACAACGTTTCTCCTTGAGATACCGCGTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCg                                                                                                                                                                                                   >  2:341843/1‑275 (MQ=255)
                                                                             tCAGCGCCTTAACCCCGCTTGAGCACACTGTGTTATGTCAACAAGCACAACGTTTCTCCTTGAGATACCGCGTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTc                                                                                                                                                                                                <  2:202899/275‑1 (MQ=255)
                                                                               aGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAGCACAACGTTTCTCCTTGAGATACCGCGTTCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTccc                                                                                                                                                                                              <  1:72617/275‑1 (MQ=255)
                                                                                              cgTGAGCACACTGTGTTATGTCAACAAGCACAACGTTTCTCCTTGAGATACCGCGTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCt                                                                                                                                                                                                                                                >  2:215110/1‑210 (MQ=255)
                                                                                              cgTGAGCACACTGTGTTATGTCAACAAGCACAACGTTTCTCCTTGAGATACCGCGTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGACATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGACTCAAGTCGCACCAAAGGGGAGCGGGAAACCGATCCCCTTTTATATTTAGCg                                                                                                                                                                               >  2:514641/1‑275 (MQ=255)
                                                                                                     acacTGTGTTATGTCAACAAGCACAACGTTTCTCCTTGAGATACCGCGTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGt                                                                                                                                                                        <  4:341842/275‑1 (MQ=255)
                                                                                                                     acaaGCACAGCGTTTCTCCTTGAGATACCGCGTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACCCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGgatgcgat                                                                                                                                                        >  4:120326/1‑275 (MQ=255)
                                                                                                                          cccaACGTTTCTCCTTGAGATACCGCGTGCACAACAGCTGGAAACAGGCAGCGGAAAGGTACGTCACCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGa                                                                                                                                                   <  2:485977/273‑1 (MQ=255)
                                                                                                                                     tcCTTGAGATACCGCGTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTg                                                                                                                                                                                                                                                                       <  6:183385/148‑1 (MQ=255)
                                                                                                                                        ttGAGATACCGCGTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCtt                                                                                                                                                                                                                                                                    >  8:1048993/1‑148 (MQ=255)
                                                                                                                                        ttGAGATACCGCGTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCtt                                                                                                                                                                                                                                                                    >  7:752478/1‑148 (MQ=255)
                                                                                                                                        ttGAGATACCGCGTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCCCAAACTCATTGTTGCAGCACTGCa                                                                                                                                                                                                                                                                           <  6:1019729/141‑1 (MQ=255)
                                                                                                                                         tGAGATACCGCGTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTg                                                                                                                                                                                                                                                                   <  6:77958/148‑1 (MQ=255)
                                                                                                                                              tACCGCGTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGa                                                                                                                                                                                                                                                                  <  6:137894/144‑1 (MQ=255)
                                                                                                                                                 cgcgTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGtt                                                                                                                                                                                                                                                           >  8:241473/1‑148 (MQ=255)
                                                                                                                                                 cgcgTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCAtc                                                                                                                            >  2:481462/1‑275 (MQ=255)
                                                                                                                                                  gcgTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTc                                                                                                                                                                                                                                                          <  6:948558/148‑1 (MQ=255)
                                                                                                                                                  gcgTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTc                                                                                                                                                                                                                                                          <  7:247774/148‑1 (MQ=255)
                                                                                                                                                      gCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGaa                                                                                                                                                                                                                                                      <  5:951016/148‑1 (MQ=255)
                                                                                                                                                       cacaACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTCGCAGCACTGCAACTGCTTGATGCTGTTc                                                                                                                                                                                                                                                          >  8:591770/1‑143 (MQ=255)
                                                                                                                                                         caacaGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGt                                                                                                                                                                                                                                                            >  8:633982/1‑139 (MQ=255)
                                                                                                                                                         caacaGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGaa                                                                                                                                                                                                                                                      <  6:247217/145‑1 (MQ=255)
                                                                                                                                                         caacaGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGaa                                                                                                                                                                                                                                                      >  5:247217/1‑145 (MQ=255)
                                                                                                                                                          aacaGCTGGCAACAGGCAGCGGAAATGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTg                                                                                                                                                                                                                                                                   >  5:1013585/1‑131 (MQ=255)
                                                                                                                                                          aacaGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGACCCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGCTGCTGACGCATCTTATCCGCCCTTACATCTCTccctg                                                                                                                    >  1:113985/1‑272 (MQ=255)
                                                                                                                                                          aacaGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATAcc                                                                                                                                                                                                                                                 <  7:87175/148‑1 (MQ=255)
                                                                                                                                                            caGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATa                                                                                                                                                                                                                                                   <  7:1071375/145‑1 (MQ=255)
                                                                                                                                                                  gCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTa                                                                                                           >  2:62666/1‑275 (MQ=255)
                                                                                                                                                                     aCAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGAt                                                                                                                                                                                                                                                                 <  8:845500/122‑1 (MQ=255)
                                                                                                                                                                     aCAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGAt                                                                                                                                                                                                                                                                 >  7:845500/1‑122 (MQ=255)
                                                                                                                                                                      cAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTc                                                                                                                                                                                                                                      <  8:1008035/148‑1 (MQ=255)
                                                                                                                                                                          cgggtgcaaGGTACGTCAGATGGCATTGGTCCTGCGCCCCCGGCGACGCGTATGAACCGGGGGGTTTTCGCCTTTTTNTTCCTCCAGTTCTTTGTTGCAGCGTTGCAACTGCTTGAAGCTGTTCTCAAAGACTTGAAATCATTCAGATTCAAGTCGCTCCAAATGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCTGGTTGTTGGCGGATGGGATGCTGACGCATCTTATCCGCCGTACCATCTCTCCCGGCAACAGTTATTGCCGCt                                                                                                   <  3:113985/268‑1 (MQ=255)
                                                                                                                                                                          cAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGTGCTGCAAcaa                                                                                                                                                                                                                                                                       >  5:936336/1‑109 (MQ=255)
                                                                                                                                                                          cAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGTGCTGCAAcaa                                                                                                                                                                                                                                                                       >  5:791543/1‑109 (MQ=255)
                                                                                                                                                                             cGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCa                                                                                                                                                                                                                               <  8:176351/148‑1 (MQ=255)
                                                                                                                                                                             cGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAAcc                                                                                                                                                                                                    >  4:350806/1‑175 (MQ=255)
                                                                                                                                                                               gAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTg                                                                                                                                                                                                                                                                       >  7:606488/1‑106 (MQ=255)
                                                                                                                                                                                aaaGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGTGCTGCAAca                                                                                                                                                                                                                                                                        <  5:570895/104‑2 (MQ=255)
                                                                                                                                                                                aaaGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGa                                                                                                                                                                                                                                              <  7:939242/130‑1 (MQ=255)
                                                                                                                                                                                aaaGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTg                                                                                                                                                                             <  2:158239/195‑1 (MQ=255)
                                                                                                                                                                                 aaGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTc                                                                                                                                                                                                <  2:350807/175‑1 (MQ=255)
                                                                                                                                                                                     tACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTaa                                                                                                                                                                                                                                         >  5:586536/1‑130 (MQ=255)
                                                                                                                                                                                         tCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGaa                                                                                                                                                                                                                                                      <  6:520372/113‑1 (MQ=255)
                                                                                                                                                                                          cAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTCTCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAg                                                                                                                                                                                                                  <  5:216523/148‑1 (MQ=255)
                                                                                                                                                                                          cAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAg                                                                                                                                                                                                                  >  6:216523/1‑148 (MQ=255)
                                                                                                                                                                                              tGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGTGCTGCa                                                                                                                                                                                                                                                                           >  6:932560/1‑87 (MQ=255)
                                                                                                                                                                                              tGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGTGCTGCa                                                                                                                                                                                                                                                                           <  5:932560/87‑1 (MQ=255)
                                                                                                                                                                                                   gTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTc                                                                                                                                                                                                                                      >  6:835019/1‑119 (MQ=255)
                                                                                                                                                                                                      cTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATAc                                                                                                                                                                                                                                                  <  7:690359/104‑1 (MQ=255)
                                                                                                                                                                                                       tCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAATTACTTCAggta                                                                                                                                                                                                                                 <  7:588341/120‑4 (MQ=255)
                                                                                                                                                                                                        ccTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAAt                                                                                                                                                                                                                                        <  5:1035303/112‑1 (MQ=255)
                                                                                                                                                                                                        ccTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGaaa                                                                                                                                                                                                      <  7:1057299/146‑1 (MQ=255)
                                                                                                                                                                                                         cTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATg                                                                                                                                                                                                                                                                <  7:53660/87‑1 (MQ=255)
                                                                                                                                                                                                         cTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATg                                                                                                                                                                                                                                                                >  8:53660/1‑87 (MQ=255)
                                                                                                                                                                                                         cTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCaaa                                                                                                                                                                                                                   <  8:1057775/132‑1 (MQ=255)
                                                                                                                                                                                                            aaCCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTAAACAATTccc                                                                  <  2:335783/274‑1 (MQ=255)
                                                                                                                                                                                                                acaGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGaaa                                                                                                                                                                                                                                                     >  4:104607/1‑91 (MQ=255)
                                                                                                                                                                                                                acaGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATAcc                                                                                                                                                                                                                                                 >  5:921547/1‑95 (MQ=255)
                                                                                                                                                                                                                acaGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATAcc                                                                                                                                                                                                                                                 <  6:921547/95‑1 (MQ=255)
                                                                                                                                                                                                                acaGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTa                                                                                                                                                                                                                                          >  8:718090/1‑102 (MQ=255)
                                                                                                                                                                                                                    gagaCGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATAcc                                                                                                                                                                                                                                                 <  2:104608/91‑1 (MQ=255)
                                                                                                                                                                                                                        cgcgTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTaa                                                                                                                                                                                                                                         >  6:498663/1‑95 (MQ=255)
                                                                                                                                                                                                                            gtcgtaaccgtggggATATCTCCATTTTTATCCTCAAACTACTGGTTGCAGCACTGGAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGAAGCAAGTCGCACCAAATGGGAGCGGGAAACCGCTCCCTTTGTATATTTATCGTGCGGGTTGGTTTCGGATGCGATGCTGATGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTAAACAATTCCCCATTCCTGGCGTatat                                                  <  4:514640/262‑1 (MQ=255)
                                                                                                                                                                                                                             aTGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTg                                                                                                                                                                                                                                                                       >  8:215043/1‑60 (MQ=255)
                                                                                                                                                                                                                             aTGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGaa                                                                                                                                                                                                                                                      <  5:511958/77‑1 (MQ=255)
                                                                                                                                                                                                                             aTGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGa                                                                                                                                                                                                                                              >  5:484220/1‑85 (MQ=255)
                                                                                                                                                                                                                             aTGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGc                                                                                                                                                                                                                         <  5:1108677/106‑1 (MQ=255)
                                                                                                                                                                                                                             aTGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCaa                                                                                                                                                                                                                    <  7:228972/111‑1 (MQ=255)
                                                                                                                                                                                                                                aaCCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGaaa                                                                                                                                                                                                      <  5:225271/122‑1 (MQ=255)
                                                                                                                                                                                                                                aaCCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGaaa                                                                                                                                                                                                      >  6:225271/1‑122 (MQ=255)
                                                                                                                                                                                                                                  ccTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAgg                                                                                                                                                                                                                 <  5:140715/109‑1 (MQ=255)
                                                                                                                                                                                                                                    tggtggATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAg                                                                                                                                                                                                                                    >  6:1094208/1‑88 (MQ=255)
                                                                                                                                                                                                                                      ggtggATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTg                                                                                                                                                                                                                                                                   <  7:174064/56‑1 (MQ=255)
                                                                                                                                                                                                                                          gATATCGCCATTCTTATCCTCAAACTCATTGTTgcagca                                                                                                                                                                                                                                                                                <  4:362530/39‑1 (MQ=255)
                                                                                                                                                                                                                                          gATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCGCTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAgg                                                                                                                                                                                                                 >  6:363591/1‑102 (MQ=255)
                                                                                                                                                                                                                                          gATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCGCTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAgg                                                                                                                                                                                                                 <  5:363591/102‑1 (MQ=255)
                                                                                                                                                                                                                                           atatCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAg                                                                                                                                                                                                                  >  6:63982/1‑100 (MQ=255)
                                                                                                                                                                                                                                            tatCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCa                                                                                                                              <  3:69893/183‑1 (MQ=255)
                                                                                                                                                                                                                                              tCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAAATGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGAGGGAAAc                                                                                                                                                                                                     >  6:576780/1‑110 (MQ=255)
                                                                                                                                                                                                                                               cGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAgg                                                                                                                                                                                                                 <  8:1096250/97‑1 (MQ=255)
                                                                                                                                                                                                                                                         aTCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAg                                                                                                                                                                                                                  >  6:767855/1‑86 (MQ=255)
                                                                                                                                                                                                                                                         aTCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTAAACAATTCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGAtt                     >  4:498928/1‑275 (MQ=255)
                                                                                                                                                                                                                                                              cAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCaa                                                                                                                                                                                                                    >  6:450697/1‑79 (MQ=255)
                                                                                                                                                                                                                                                              cAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCaa                                                                                                                                                                                                                    <  5:450697/79‑1 (MQ=255)
                                                                                                                                                                                                                                                              cAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGc                                                                                                                                                                                                  <  8:88942/97‑1 (MQ=255)
                                                                                                                                                                                                                                                               aaaCTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTAAACAATTCCCCATTCCTGGCGTATATCTGGCTCACATTCATCAATGTGATAGATTCCtctc               >  2:426240/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                 attCACTGTTGAAGCACGTCAACTGCTTGAGGCTGTACTGAAATACCTGAAGTAATGTAGCTTCAAGTCGCACCAGAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTGGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGATAAACAATTCCCCATGCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTccc             <  4:13099/273‑1 (MQ=255)
                                                                                                                                                                                                                                                                            cagcaCTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTAAACAATTCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATggg  <  2:120327/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                              gcaCTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAgg                                                                                                                                                                                                                 >  5:76313/1‑66 (MQ=255)
                                                                                                                                                                                                                                                                               caCTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAg                                                                                                                                                                                                             <  7:204648/69‑1 (MQ=255)
                                                                                                                                                                                                                                                                                aCTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCg                                                                                                                                                                                                           >  7:1041858/1‑70 (MQ=255)
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                              
TCACTGTACGCAATGTGTTATGCGGGGGCCGCATCGTTACCCGGCGCACTAAGTCCTGGCTGAAACGGGTGGTGCCGTCAGCGCCTTAACCCCGCGTGAGCACACTGTGTTATGTCAACAAGCACAACGTTTCTCCTTGAGATACCGCGTGCACAACAGCTGGCAACAGGCAGCGGAAAGGTACGTCAGCTGGCAGTGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTAAACAATTCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGG  >  NC_000913/3853332‑3853873

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