Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,139,727 T→C 5.1% intergenic (+242/+32) yegH → / ← asmA inner membrane protein/suppressor of OmpF assembly mutants; putative outer membrane protein assembly factor; inner membrane‑anchored periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,139,7270TC 5.1% 269.1 / 3.1 78intergenic (+242/+32)yegH/asmAinner membrane protein/suppressor of OmpF assembly mutants; putative outer membrane protein assembly factor; inner membrane‑anchored periplasmic protein
Reads supporting (aligned to +/- strand):  ref base T (39/35);  new base C (2/2);  total (41/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

TAAGGATGGCGAGATGGAGTACGAGGTGTGATGGTATCGGTCATTTTACCCAATAGCCGAAAGGAATAAGATGAGGGAACCGCGCCATCGGCGATTTACCGCCGGGAGCCCGGGTTGCCAGGGTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCAGTAACTGATCCACTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCGTTC  >  NC_000913/2139455‑2139975
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                        
tAAGGATGGCGAGATGGAGTACGAGGTGTGATGGTATCGGTCATTTTACCCAATAGCCGAAAGGAATAAGATGAGGGAACCGCGCCATCGGCGATTTACCGCCGGGAGCCCGGGTTGCCAGGGTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTa                                                                                                                                                                                                                                                         >  4:18854/1‑274 (MQ=255)
 aaGGATGGCGGGATGGAGTACGGGGGTGGATGGTATCGGTCATTTTACCCAATAGCCGAAAGGAATAAGATGAGGGAACCGCGCCAGCGGCGATTTACCGCCGGGAGCCCGGGTTGCCAGGGTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTAt                                                                                                                                                                                                                                                        <  4:124552/274‑1 (MQ=255)
   ggATGGCGAGATGGAGTACGAGGTGTGATGGTATCGGTCATTTTACCCAATAGCCGAAAGGAATAAGATGAGGGAACCGCGCCATCGGCGATTTACCGCCGGGAGCCCGGGTTGCCAGGGTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAg                                                                                                                                                                                                                                                     <  4:211572/275‑1 (MQ=255)
   ggATGGCGAGAGGGAGTACGAGTGGTGATGGTATCGGTCATTTTACCCAATAGCCGAAAGGAATAAGATGAGGGAACCGGGCCATCGGCGATTTACCGCCGGGAGCCCGGGTTGCCAGGGTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAg                                                                                                                                                                                                                                                     <  1:114422/275‑1 (MQ=255)
               tgggtgcgagtggtgACGGTATCGGTCATTTTACCCAATAACCGGAAGGAATAAGATGAGGGAACCGCGCCATCGGGGATTTACCGCCGGGAGGCCGGGTTGCCAGGGTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGAc                                                                                                                                                                                                                                          <  2:529541/262‑1 (MQ=255)
                      gAGGTGTGATGGTAGCGGTCATGTTACCCAATAGCCGAAAGGAATAAGATGAGGGAACCGCGCCATCGGCGATTTACCGCCGGGAGCCCGGGTTGCCAGGGTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCAc                                                                                                                                                                                                                                   <  2:532817/274‑1 (MQ=255)
                          tgtgATGGTATCGGTCATTTTACCCAATAGCCGAAAGGAATAAGATGAGGGAACCGCGCCATCGGCGATTTACCGCCGGGAGCCCGGGTTGCCAGGGTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTgg                                                                                                                                                                                                                               <  2:18854/274‑1 (MQ=255)
                           gtgATGGTATCGGTCATTTTACCCAATAGCCGAAAGGAATAAGATGAGGGAACCGCGCCATCGGCGATTTACCGCCGGGAGCCCGGGTTGCCAGGGTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCCGGCCTACAGGACTGGTCACCTGGt                                                                                                                                                                                                                              >  2:344882/1‑274 (MQ=255)
                           gtgATGGTATCGGTCATTTTACCCAATAGCCGAAAGGAATAAGATGAGGGAACCGCGCCATCGGCGATTTACCGCCGGGAGCCCGGGTTGCCAGGGTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGt                                                                                                                                                                                                                              <  2:169148/274‑1 (MQ=255)
                                                                                                                      aggggTTGTGGAGATTGGGCCACCCTGGNNNTTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGGGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCt                                                                                                                                   <  3:149225/272‑1 (MQ=255)
                                                                                                                              tggAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTa                                                                                                                                                                                                                                                         >  7:1001813/1‑148 (MQ=255)
                                                                                                                                 agaTTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCa                                                                                                                                                                                                                                                      <  8:1001813/148‑1 (MQ=255)
                                                                                                                                   aTTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCg                                                                                                                                                                  >  2:83959/1‑230 (MQ=255)
                                                                                                                                       aGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCg                                                                                                                                                              <  4:83958/230‑1 (MQ=255)
                                                                                                                                       aGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGttt                                                                                                                  >  2:533225/1‑274 (MQ=255)
                                                                                                                                           aCCCTGGCACGTGCACCGGCTTGGGCGTTCCCAAAGAGCGAGGAACATAAAGTGAACGGAAACGGCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCa                                                                                                             <  4:176492/275‑1 (MQ=255)
                                                                                                                                               tGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGc                                                                                                                                                                                                                                                   >  5:799071/1‑137 (MQ=255)
                                                                                                                                               tGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGc                                                                                                                                                                                                                                                   <  6:799071/137‑1 (MQ=255)
                                                                                                                                                      ttCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCa                                                                                                                                                                                                                                    <  8:276488/145‑1 (MQ=255)
                                                                                                                                                      ttCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCa                                                                                                                                                                                                                                    >  7:276488/1‑145 (MQ=255)
                                                                                                                                                           cGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTaa                                                                                                                                                                                                                            >  5:99006/1‑148 (MQ=255)
                                                                                                                                                                 gAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACa                                                                                                                                                                                                                      <  6:99006/148‑1 (MQ=255)
                                                                                                                                                                         cATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATc                                                                                                                                                                                                                                                       >  8:898280/1‑107 (MQ=255)
                                                                                                                                                                         cATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATc                                                                                                                                                                                                                                                       <  7:898280/107‑1 (MQ=255)
                                                                                                                                                                                 cGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTCCGCAGTAACTGATCCAATTGCAGGCTGTAATTTAGTtgctg                                                                        >  1:38952/1‑274 (MQ=255)
                                                                                                                                                                                   aGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTTGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCAGTAACTGATCCACTTGCAGGCTGTAATTGAGTTGCTGcc                                                                      >  2:253673/1‑274 (MQ=255)
                                                                                                                                                                                          tCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATAAGACGCGTAAAGCGTCGCATCAGGc                                                                                                                                                                                                                                                   <  6:156822/94‑1 (MQ=25)
                                                                                                                                                                                          tCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATAAGACGCGTAAAGCGTCGCATCAGGc                                                                                                                                                                                                                                                   >  5:156822/1‑94 (MQ=25)
                                                                                                                                                                                           cAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGTCCTGTAGGCCTgataag                                                                                                                                                                                                                                   <  6:408594/109‑6 (MQ=38)
                                                                                                                                                                                           cAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGTCCTGTAGGCCTgataag                                                                                                                                                                                                                                   >  5:408594/1‑104 (MQ=38)
                                                                                                                                                                                           cAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCt                                                                                                                                                                                            >  6:45569/1‑148 (MQ=255)
                                                                                                                                                                                           cAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCt                                                                                                                                                                                            <  5:45569/148‑1 (MQ=255)
                                                                                                                                                                                                gAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCa                                                                                                                                                                                                                                    <  6:706060/103‑1 (MQ=255)
                                                                                                                                                                                                gAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCa                                                                                                                                                                                                                                    >  5:706060/1‑103 (MQ=255)
                                                                                                                                                                                                     gAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTa                                                                                                                                                                                                                             >  5:896970/1‑105 (MQ=255)
                                                                                                                                                                                                     gAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTa                                                                                                                                                                                                                             <  6:896970/105‑1 (MQ=255)
                                                                                                                                                                                                        cAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCaa                                                                                                                                                                                                         <  7:445966/122‑1 (MQ=255)
                                                                                                                                                                                                        cAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCaa                                                                                                                                                                                                         >  8:445966/1‑122 (MQ=255)
                                                                                                                                                                                                            cTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTa                                                                                                                                                                                                                                                         <  5:282641/70‑1 (MQ=255)
                                                                                                                                                                                                            cTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTa                                                                                                                                                                                                                                                         >  6:282641/1‑70 (MQ=255)
                                                                                                                                                                                                            cTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACAtcttct                                                                                                                                                                                                                >  5:309547/1‑111 (MQ=255)
                                                                                                                                                                                                            cTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACAtcttct                                                                                                                                                                                                                <  6:309547/111‑1 (MQ=255)
                                                                                                                                                                                                              tCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTAc                                                                                                                                                                                                                                               <  5:684479/78‑1 (MQ=255)
                                                                                                                                                                                                              tCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTAc                                                                                                                                                                                                                                               >  6:684479/1‑78 (MQ=255)
                                                                                                                                                                                                              tCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACa                                                                                                                                                                                                                      <  8:104624/103‑1 (MQ=255)
                                                                                                                                                                                                              tCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACa                                                                                                                                                                                                                      >  7:104624/1‑103 (MQ=255)
                                                                                                                                                                                                                gATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCAGTAACTGATCCACTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCg                                        >  2:293749/1‑275 (MQ=255)
                                                                                                                                                                                                                          tttCACAAGATAAGGATCCCGCCCCGAACGGCGGATGCGACGCTTTACGCGTCTTATCAGGCCTACCGGACTGGTAACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCCCATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCAGTAACTGATCCACTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCa                               <  4:18657/273‑1 (MQ=255)
                                                                                                                                                                                                                            tCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTaa                                                                                                                                                                                                                            >  8:1074436/1‑83 (MQ=255)
                                                                                                                                                                                                                            tCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTaa                                                                                                                                                                                                                            <  7:1074436/83‑1 (MQ=255)
                                                                                                                                                                                                                            tCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACa                                                                                                                                                                                                                      <  6:134124/89‑1 (MQ=255)
                                                                                                                                                                                                                            tCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACa                                                                                                                                                                                                                      >  5:134124/1‑89 (MQ=255)
                                                                                                                                                                                                                            tCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCTACGTTTCGCTTCGTCTTGTAGTTGTTTGCTCAGTAACTGACCCATGTGCAGGCTGTAATTGAGTTGCTGCCAGTTGTCATAAACCCGAAGCGGTACTGGCGTTTTTTTCAGGa                            >  4:33939/1‑275 (MQ=255)
                                                                                                                                                                                                                            tCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCAGTAACTGATCCACTTGCAGGCTGTAAt                                                                                  >  2:478327/1‑221 (MQ=255)
                                                                                                                                                                                                                            tCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCAGTAACTGATCCACTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGtttct                                   >  2:18657/1‑268 (MQ=255)
                                                                                                                                                                                                                                aaTATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTTCGCAGTAACTGATCCACTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAaat                         >  2:254306/1‑274 (MQ=255)
                                                                                                                                                                                                                                aaTATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCAGTAACTGATCCACTTGCAGGCTGTAATTGAg                                                                              <  4:478326/221‑1 (MQ=255)
                                                                                                                                                                                                                                aaTATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCAGTAACTGATCCACTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTCCTTTCAGGAaat                         >  2:471150/1‑274 (MQ=255)
                                                                                                                                                                                                                                          tgCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCAGTAACTGATCCACTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAAATCAATCAGttt               >  2:524712/1‑274 (MQ=255)
                                                                                                                                                                                                                                          tgCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCAGTAACTGATCCACTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAAATCAATCAGTTTg              >  2:386840/1‑275 (MQ=255)
                                                                                                                                                                                                                                            ctgccCCGAAGCCCTGATGCGCGGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTTGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCAGTAACTGATCCACTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAAATCAATCAGTTTGCt            <  4:524711/273‑1 (MQ=255)
                                                                                                                                                                                                                                            ccgccACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCAGTAACTGATCCACTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAAATCAATCAGTTTGc             >  4:220238/1‑274 (MQ=255)
                                                                                                                                                                                                                                                   gAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCAGTAACTGATCCACTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCg     <  4:386839/275‑1 (MQ=255)
                                                                                                                                                                                                                                                    aaCGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCAGTAACTGATCCACTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCg     >  2:17790/1‑274 (MQ=255)
                                                                                                                                                                                                                                                       gCCTGAAGCGACGCTTTACGCGTCTTATCAGGCCTACAGGAATGGAAACCTAGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCAGTAACTGATCCACTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCGTTc  <  2:467599/274‑1 (MQ=255)
                                                                                                                                                                                                                                                        ccTGATAAGACGCGTAAAGCGTCGCATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGc                                                                                                                                    >  6:133450/1‑143 (MQ=255)
                                                                                                                                                                                                                                                        ccTGATAAGACGCGTAAAGCGTCGCATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGc                                                                                                                                    <  5:485251/143‑1 (MQ=255)
                                                                                                                                                                                                                                                          tGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAActtct                                                                                                                                                                                                    <  6:104768/77‑1 (MQ=255)
                                                                                                                                                                                                                                                          tGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAActtct                                                                                                                                                                                                    >  5:104768/1‑77 (MQ=255)
                                                                                                                                                                                                                                                                gACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCt                                                                                                                                                                                                                                 >  6:652141/1‑42 (MQ=255)
                                                                                                                                                                                                                                                                gACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCt                                                                                                                                                                                                                                 <  5:652141/42‑1 (MQ=255)
                                                                                                                                                                                                                                                                 aCGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCttt                                                                                                                                                                        >  7:752198/1‑98 (MQ=255)
                                                                                                                                                                                                                                                                 aCGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCttt                                                                                                                                                                        <  8:752198/98‑1 (MQ=255)
                                                                                                                                                                                                                                                                 aCGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCAGTAACTGATCCACt                                                                                               >  2:538420/1‑171 (MQ=255)
                                                                                                                                                                                                                                                                     tttACGCGTCTTATCAGGNCTACNGGACTGGTCACCTGGTAACCTACNTCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCAGTAACTGATCCACTTGCa                                                                                           <  4:538419/171‑1 (MQ=255)
                                                                                                                                                                                                                                                                      ttACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTg                                                                                                                 >  7:1047576/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                             tCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCAGTa                                                                                                          >  6:322313/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                              cTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGttt                                                                                                                                       >  2:566439/1‑118 (MQ=255)
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                        
TAAGGATGGCGAGATGGAGTACGAGGTGTGATGGTATCGGTCATTTTACCCAATAGCCGAAAGGAATAAGATGAGGGAACCGCGCCATCGGCGATTTACCGCCGGGAGCCCGGGTTGCCAGGGTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCAGTAACTGATCCACTTGCAGGCTGTAATTGAGTTGCTGCCAGTTGCCATAAACCCGCAGCGGTACTGGCGTTTCTTTCAGGAAATCAATCAGTTTGCTTTCCCCGTTC  >  NC_000913/2139455‑2139975

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