Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,830,727 C→T 6.1% Q91* (CAA→TAA)  yicH → putative inner membrane‑anchored periplasmic AsmA family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,830,7270CT 6.1% 209.8 / 2.6 68Q91* (CAA→TAA) yicHputative inner membrane‑anchored periplasmic AsmA family protein
Reads supporting (aligned to +/- strand):  ref base C (33/29);  new base T (2/2);  total (36/32)
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.22e-01

CGACTATCATCTGGCCTTCGGGGCGATGGATCACCGTTTTTCCGCGCCATCTCATATCGTGCTGGAGAACGTCACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAA  >  NC_000913/3830585‑3830873
                                                                                                                                              |                                                                                                                                                  
cGACTATCATCTGGCCTTCGGGGCGATGGATCACCGTTTTTCCGCGCCATCTCATATCGTGCTGGAGAACGTCACGTTTGGTCGTGATGGCCAGCCAGCGACCCTGGTGGACAAAAGTGTCGACATTGCGCTAAGCAGTCTGCAACTGa                                                                                                                                              >  7:301000/1‑149 (MQ=255)
    tATCATCTGGCCTTCGGGGCGATGGATCACCGTTTTTCCGCGCCATCTCATATCGTGCTGGAGAACGTCACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCg                                                                                                                                           >  8:354372/1‑148 (MQ=255)
      tcatcTGGCCTTCGGGGCGATGGATCACCGTTTTTCCGCGCCATCTCACATCGTGCTGGAGAACGTCACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGaa                                                                                                                                         <  5:36847/148‑1 (MQ=255)
        atcTGGCCTTCGGGGCGATGGATCACCGTTTTTCCGCGCCATCTCATATCGTGCTGGAGAACGTCACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCa                                                                                                                                      >  6:137064/1‑149 (MQ=255)
        atcTGGCCTTCGGGGCGATGGATCACCGTTTTTCCGCGCCAGCTCATATCGTGCTGGAGAACGTCACGTTTGGTCGTGATGGCCAGCCCGCGACCATGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCa                                                                                                                                      <  3:206128/149‑1 (MQ=255)
         tcTGGCCTTCCGGGCCATCGATCACCGATTTTACGCGCCATCTCATATCGTCCTGGAGAACGACACGTTTGGTCGTGATAGCCAGCCAGCGACCCTGCAGGCCAAAACAGTCGACATAGCGCTAAGCAGTCGGCAACTGAACAAACCAc                                                                                                                                     >  4:263038/1‑149 (MQ=255)
            ggCCTTCGGGGCGATGGATCACCGTTTTTCCGCGCCATCTCATATCGTGCTGGAGAACGTCACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGc                                                                                                                                   >  4:111397/1‑148 (MQ=255)
             ggcTTCGGGGCGTTGGATCACCGTTTTTCCGCGGCATCTCATATCGTGCTGGAGAACGGCACGTTTGGTCGTGTTGGCCCGCCCGCGACCCTGGTGGCCAAAAGTGTCGTCATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGcc                                                                                                                                  <  6:359188/146‑1 (MQ=255)
             gCCTTCGGGGCGATGGATCACCGTTTTTCCGCGCCATCTCATATCGTGCTGGAGAACGTCACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGc                                                                                                                                   <  7:354372/147‑1 (MQ=255)
                   ggggCGATGGATCACCGTTTTTCCGCGCCATCTCATATCGTGCTGGAGAACGTCACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGa                                                                                                                           <  6:81298/149‑1 (MQ=255)
                     ggCGATGGATCACCGTTTTTCCGCGCCATCTCATATCGTGCTGGAGAACGTCACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATa                                                                                                                         <  8:153849/149‑1 (MQ=255)
                       cGATGGATCACCGTTTTTCCGCGCCATCTCATATCGTGCTGGAGAACGTCACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATAcc                                                                                                                       <  4:327804/149‑1 (MQ=255)
                        gATGGATCACCGTTTTTCCGCGCCATCTCATATCGTGCTGGAGAACGTCACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCa                                                                                                                      >  6:149219/1‑149 (MQ=255)
                           ggATCACCGTTTTTCCGCGCCATCTCATATCGTGCTGGAGAACGTCACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGATAAGCAGTCGGCAACTTCACGAACCACGCCATGTCGAGACCATcc                                                                                                                   >  2:352936/1‑149 (MQ=255)
                               cACCGTTTTTCCGCGCCATCTCATATCGTGCTGGAGAACGTCACGTTTGGTCGTGATGGCCCGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTACGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCctgct                                                                                                               >  2:257679/1‑149 (MQ=255)
                                aCCGTTTTTCCGCGCCATCTCATATCGTGCTGGAGAACGTCACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCa                                                                                                                                 >  6:76146/1‑130 (MQ=255)
                                 ccGTTTTTCCGCGCCATCTCATATCGTGCTGGAGAACGTCACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTgg                                                                                                             <  3:111397/149‑1 (MQ=255)
                                        tCCGCGCCATCTCATATCGTGCTGGAGAACGTCACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACg                                                                                                       <  3:263038/148‑1 (MQ=255)
                                          cgcgCCATCTCATATCGTGCTGGAGAACGTCACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGa                                                                                                    <  1:46445/149‑1 (MQ=255)
                                                 tctcATATCGTGCTGGAGAACGTCACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGaa                                                                                             >  7:77221/1‑149 (MQ=255)
                                                   tcATATCGTGCTGGAGAACGTCACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAAtc                                                                                           >  7:52875/1‑149 (MQ=255)
                                                       atCGTGCTGGAGAACGTCACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCAc                                                                                       >  6:280127/1‑149 (MQ=255)
                                                        tCGTGCTGGAGAACGTCACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCAcc                                                                                      <  5:280127/149‑1 (MQ=255)
                                                        tCGTGCTGGACAACGTCACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGGGGCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCAcc                                                                                      <  1:21626/149‑1 (MQ=255)
                                                           tGCTGGAGAACGTCACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGAc                                                                                   <  7:131740/149‑1 (MQ=255)
                                                              tGGAGAACGTCACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCAc                                                                                       >  2:348492/1‑142 (MQ=255)
                                                              tGGAGAACGTCACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCAc                                                                                       <  1:348492/142‑1 (MQ=255)
                                                              tGGAGAACGTCACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGAc                                                                                   >  5:72596/1‑146 (MQ=255)
                                                                gagaACGTCACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGTGCTAAGCAGTCGGCAACTGACCGAACCGCGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAgcc                                                                              >  8:160825/1‑147 (MQ=255)
                                                                gagaACGGCACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGAc                                                                              <  2:66934/149‑1 (MQ=255)
                                                                          cGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTa                                                                    <  5:137064/149‑1 (MQ=255)
                                                                           gTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAAACACGCCATGTCGATACAATCCTGCTGGAAATCGGGACGCTGAATCTCACCGACCAGAACGCGCCGCTAc                                                                   >  7:379227/1‑149 (MQ=255)
                                                                                 tCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTc                                                              >  2:238351/1‑148 (MQ=255)
                                                                                 tCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTc                                                              >  8:82226/1‑148 (MQ=255)
                                                                                  cGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCa                                                             >  3:266795/1‑148 (MQ=255)
                                                                                       tGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCAAAGCCg                                                       <  1:238351/149‑1 (MQ=255)
                                                                                           cAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCACAGCCGATCg                                                   >  3:73952/1‑149 (MQ=255)
                                                                                              cccGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCAGGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCAAAGCCGATCGTCt                                                <  4:326400/149‑1 (MQ=255)
                                                                                                cgcgACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCCCGCCATGTCGATACCATCCTGCTGGAACACGGGACGCTGACTCTCACCGACCAGCCCGCGCCGCTACCGTTCAACGCCGATCGTCTg                                               >  4:94295/1‑148 (MQ=255)
                                                                                                   gACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCAAc                                           <  1:244153/149‑1 (MQ=255)
                                                                                                       cTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCg                                       >  7:7593/1‑149 (MQ=255)
                                                                                                        tggtggCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGATGAATCTAAACTACCAGAACTCGCCGCTACCTTTCAAAGCCGATCGTCTGCAACTGCGt                                      >  4:246811/1‑149 (MQ=255)
                                                                                                             gCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCaa                                            >  8:349749/1‑138 (MQ=255)
                                                                                                             gCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCaa                                            <  7:349749/138‑1 (MQ=255)
                                                                                                             gCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGata                                  <  8:77221/148‑1 (MQ=255)
                                                                                                              ccAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATg                                >  6:144548/1‑149 (MQ=255)
                                                                                                               cAAAAGTGTCGACATTGCGCTAAGCAGTCGGGAACTGACCGATCCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATGGACCAGACCAGACCGCTCCGCTCCCGTTCAACGCCGATCGTCTGCCACTGCGTGatattg                               >  6:106856/1‑147 (MQ=255)
                                                                                                               cAAAAGTGTCGACATTGCGCTAAGCAGACGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACAGCGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATgg                               >  4:195507/1‑149 (MQ=255)
                                                                                                                    gtgtCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTTCGTGATATGGAGttt                          <  6:72596/149‑1 (MQ=255)
                                                                                                                      gtCGACATTGCGCTAAGCAGTCGGAAACTGACCGAACCACGCAATGTCGATCCCATCCTGCTGGAAAACGGTACGATGTATCTCACCGACCAGACCGCTCCGCTACCGTTCAAAGCCGATCTTCTGCAACTGCGTGATATGGCGTTTaa                        <  8:301000/149‑1 (MQ=255)
                                                                                                                       tCGACATTGCGCTAAGCAGTCGGCAACTGACCGACCCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCACCTGCGTGATATGGCGTTTaa                        >  1:371849/1‑148 (MQ=255)
                                                                                                                       tCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTaa                        >  8:46501/1‑148 (MQ=255)
                                                                                                                       tCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAAt                       <  5:149219/149‑1 (MQ=255)
                                                                                                                       tCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATACTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGACGCAACCGTGCAAAGCCGATCGTCTGCAACGGAGTGATATGGCGTTTaa                        >  5:150390/1‑148 (MQ=255)
                                                                                                                        cGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATa                      <  3:289317/149‑1 (MQ=255)
                                                                                                                         ccgtttGCGCTAAGCAGTCGGTAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAg                     >  4:303023/5‑149 (MQ=255)
                                                                                                                         ccgtttGCGCTAAGCAGTCGGTAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAg                     >  2:222298/5‑149 (MQ=255)
                                                                                                                                 gcTAAGCAGTCGGTAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATa             <  3:303023/149‑1 (MQ=255)
                                                                                                                                 gcTAAGCAGTCGGTAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATa             <  1:222298/149‑1 (MQ=255)
                                                                                                                                     aGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGa         <  1:352936/149‑1 (MQ=255)
                                                                                                                                       cAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCa                                                                                 >  7:254162/1‑75 (MQ=255)
                                                                                                                                       cAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCa                                                                                 <  8:254162/75‑1 (MQ=255)
                                                                                                                                       cAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAAt       >  1:292077/1‑149 (MQ=255)
                                                                                                                                         gTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGCCCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATg      >  8:353501/1‑148 (MQ=255)
                                                                                                                                          tCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGa    <  3:110085/149‑1 (MQ=255)
                                                                                                                                          tCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGAGTGATATGGCGTTTAATAGCCCGAATAGCGAATgg     >  2:173592/1‑148 (MQ=255)
                                                                                                                                            ggCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGaaa  <  7:82226/149‑1 (MQ=255)
                                                                                                                                            ggCAACTGACCGAACCACCCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGaaa  <  7:172535/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
CGACTATCATCTGGCCTTCGGGGCGATGGATCACCGTTTTTCCGCGCCATCTCATATCGTGCTGGAGAACGTCACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAA  >  NC_000913/3830585‑3830873

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑