Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,148,503 G→T 5.1% intergenic (‑186/+29) yijO ← / ← eptC putative DNA‑binding transcriptional regulator YijO/phosphoethanolamine transferase EptC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,148,5030GT 5.1% 377.5 / 5.0 112intergenic (‑186/+29)yijO/eptCputative DNA‑binding transcriptional regulator YijO/phosphoethanolamine transferase EptC
Reads supporting (aligned to +/- strand):  ref base G (49/57);  new base T (2/4);  total (51/61)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.38e-01

Only 100 of 115 total aligned reads displayed.

CTGGCATTATTAAAAAAGGCCACCTGAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACG  >  NC_000913/4148354‑4148651
                                                                                                                                                     |                                                                                                                                                    
cTGGCATTATTAAAAAAGGCCACCTGAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCg                                                                                                                                                      >  2:176460/1‑150 (MQ=255)
 tGGCATTATTAAAAAAGGCCACCTGAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGc                                                                                                                                                     >  1:250632/1‑150 (MQ=255)
  ggCATTATTAAAAAAGGCCACCTGAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCa                                                                                                                                                    >  2:1356327/1‑150 (MQ=255)
   gCATTATTAAAAAAGGCCACCTGAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGGATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCAt                                                                                                                                                   <  1:680400/150‑1 (MQ=255)
     attattAAAAAAGGCCACCTGAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATcc                                                                                                                                                 <  1:57806/150‑1 (MQ=255)
       tattaAAAAAGGCCACCTGAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCgg                                                                                                                                               >  1:1415748/1‑150 (MQ=255)
       tattaAAAAAGGCCACCTGAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGc                                                                                                                                              <  1:78424/150‑1 (MQ=255)
        attaAAAAAGGCCACCTGAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCGGCGATCTGTAGGCCGGATGCGACGCTTGACGCGTCTTATCCGGc                                                                                                                                              <  2:4070/150‑1 (MQ=255)
        attaAAAAAGGCCACCTGAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGc                                                                                                                                              >  2:832496/1‑150 (MQ=255)
          taAAAAAGGCCACCTGAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATGCGACGCTTGACGCGTCTTATCCGGcct                                                                                                                                            <  1:1174870/150‑3 (MQ=255)
          taAAAAAGGCCACCTGAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATGCGACGCTTGACGCGTCTTATCCGGcct                                                                                                                                            <  1:1318961/150‑3 (MQ=255)
          taAAAAAGGCCACCTGAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCaa                                                                                                                                            <  1:213115/150‑1 (MQ=255)
           aaaaaaGGCCACCTGAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATGCGACGCTTGACGCGTCTTATcc                                                                                                                                                 >  1:717834/1‑144 (MQ=39)
           aaaaaaGGCCACCTGAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATGCGACGCTTGACGCGTCTTATcc                                                                                                                                                 <  2:717834/144‑1 (MQ=39)
            aaaaaGGCCACCTGAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAAt                                                                                                                                          >  2:1330598/1‑150 (MQ=255)
                 ggCCACCTGAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAgcgc                                                                                                                                     >  1:1564580/1‑150 (MQ=255)
                  gCCACCTGAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGcc                                                                                                                                    >  2:1558647/1‑150 (MQ=255)
                   ccACCTGAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCt                                                                                                                                   >  2:1464727/1‑150 (MQ=255)
                    cACCTGAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTg                                                                                                                                  <  1:1330951/150‑1 (MQ=255)
                     aCCTGAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTgg                                                                                                                                 >  1:33895/1‑150 (MQ=255)
                       cTGAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCTTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCt                                                                                                                               <  2:412821/150‑1 (MQ=255)
                       cTGAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCt                                                                                                                               >  2:43622/1‑150 (MQ=255)
                        tGAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTg                                                                                                                              <  2:440085/150‑1 (MQ=255)
                         gAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGa                                                                                                                             >  2:1054983/1‑150 (MQ=255)
                          aaaaGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGAt                                                                                                                            <  2:1415748/150‑1 (MQ=255)
                           aaaGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTGGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATa                                                                                                                           >  2:1058993/1‑150 (MQ=255)
                              gtgtGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAAtt                                                                                                                        <  2:1626567/150‑1 (MQ=255)
                               tgtgGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTa                                                                                                                       >  1:1353458/1‑150 (MQ=255)
                                gtgGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTAc                                                                                                                      >  1:1046859/1‑150 (MQ=255)
                                   gCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTTGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGa                                                                                                                   <  1:980203/150‑1 (MQ=255)
                                    ccTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGAt                                                                                                                  <  1:1216648/150‑1 (MQ=255)
                                    ccTGAGGGGAGTTAGATCTAAACATAGCGGCAGTTGTATTCCCGGATGTGGCGTAAACGCCTTATCCAGCATACCGCCGCGATCTGTACGCCGGATAACACGCGTCAACAGTCTCATCCGCCAAATAGCGCCTGGCTGAAAATTACTGAt                                                                                                                  >  2:1061003/1‑150 (MQ=255)
                                     cTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGAtt                                                                                                                 >  2:335987/1‑150 (MQ=255)
                                     cTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGAtt                                                                                                                 <  1:1127534/150‑1 (MQ=255)
                                           ggAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCAc                                                                                                           >  1:417113/1‑150 (MQ=255)
                                            gAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCAcc                                                                                                          <  2:1353458/150‑1 (MQ=255)
                                               ttCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGa                                                                                                       <  1:661983/150‑1 (MQ=255)
                                                  tatttAAACATTGCGTCAGTTGTATTGCCGGATGTGTCTTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGTCCGGATAAGACGCGTCAAGCGTCGCATCCGTCAAATAGCGCCTTGCTGATAATTACTGATTACCCACCTGATCg                                                                                                    <  1:499775/146‑1 (MQ=255)
                                                  gATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCg                                                                                                    >  2:135363/1‑150 (MQ=255)
                                                      tAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCAt                                                                                                <  1:1361648/150‑1 (MQ=255)
                                                       aaaCATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTGGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATa                                                                                               <  1:1058993/150‑1 (MQ=255)
                                                          cATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAggg                                                                                            >  1:326828/1‑150 (MQ=255)
                                                          aaTTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAggg                                                                                            <  1:627424/149‑1 (MQ=255)
                                                           aTTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGc                                                                                           <  2:828052/150‑1 (MQ=255)
                                                             tGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAg                                                                                         <  2:944203/150‑1 (MQ=255)
                                                             tGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAg                                                                                         <  1:939701/150‑1 (MQ=255)
                                                             tGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAg                                                                                         >  1:3587/1‑150 (MQ=255)
                                                              gCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAgt                                                                                        >  1:588461/1‑150 (MQ=255)
                                                               cGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAgtg                                                                                       <  2:1407551/150‑1 (MQ=255)
                                                               cGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAgtg                                                                                       <  1:500098/150‑1 (MQ=255)
                                                                  cAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAgtgtgt                                                                                    >  2:380830/1‑150 (MQ=255)
                                                                    gTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCg                                                                                  <  1:1558647/150‑1 (MQ=255)
                                                                     ttGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGt                                                                                 <  1:176460/150‑1 (MQ=255)
                                                                        tATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAAt                                                                              <  2:617873/150‑1 (MQ=255)
                                                                             ccGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATc                                                                         <  1:832496/150‑1 (MQ=255)
                                                                                     ggCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGtt                                                                 >  2:817777/1‑150 (MQ=255)
                                                                                      gCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTCCGttt                                                                >  2:1407576/1‑150 (MQ=255)
                                                                                       cGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGtttt                                                               >  1:1117311/1‑150 (MQ=255)
                                                                                         tAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCCAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCAGAGGGCAGTGTGTCGTAATCGCTCAGTGCGTGTTTc                                                             >  2:528463/1‑150 (MQ=255)
                                                                                         tAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTc                                                             >  1:1058821/1‑150 (MQ=255)
                                                                                         tAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTc                                                             <  1:1464727/150‑1 (MQ=255)
                                                                                         tAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTc                                                             <  2:1450170/150‑1 (MQ=255)
                                                                                          aaaCGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCCGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCt                                                            <  1:598242/150‑1 (MQ=255)
                                                                                           aaCGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCtt                                                           <  1:115055/150‑1 (MQ=255)
                                                                                              gCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATc                                                                             >  2:213940/1‑129 (MQ=255)
                                                                                              gCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATc                                                                             <  1:213940/129‑1 (MQ=255)
                                                                                                  tATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACggg                                                    <  2:831509/150‑1 (MQ=255)
                                                                                                   aTCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCATCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGt                                                   >  2:116545/1‑150 (MQ=255)
                                                                                                    tCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGtt                                                  <  1:886707/150‑1 (MQ=255)
                                                                                                       aGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTAcc                                               <  2:1564580/150‑1 (MQ=255)
                                                                                                         ccTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCaa                                             <  1:182501/150‑1 (MQ=255)
                                                                                                          cTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAAt                                            >  2:491919/1‑150 (MQ=255)
                                                                                                             ccgccgCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCa                                         <  2:306131/150‑1 (MQ=255)
                                                                                                             ccgccgCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCa                                         >  2:1015738/1‑150 (MQ=255)
                                                                                                                    gATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTa                                  >  2:176301/1‑150 (MQ=255)
                                                                                                                     aTCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAg                                 >  1:1130848/1‑150 (MQ=255)
                                                                                                                           aGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCtt                           <  2:3587/150‑1 (MQ=255)
                                                                                                                            ggCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCttt                          <  1:1172550/150‑1 (MQ=255)
                                                                                                                            ggCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCttt                          >  1:405913/1‑150 (MQ=255)
                                                                                                                             gCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTg                         <  1:135363/150‑1 (MQ=255)
                                                                                                                             gCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTg                         <  2:33895/150‑1 (MQ=255)
                                                                                                                              ccGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGa                        <  2:1515342/150‑1 (MQ=255)
                                                                                                                               cGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGaa                       <  1:1054983/150‑1 (MQ=255)
                                                                                                                               cGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGaa                       <  2:326828/150‑1 (MQ=255)
                                                                                                                                   tAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGc                   <  1:380830/150‑1 (MQ=255)
                                                                                                                                     aGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTta                                                          >  2:405009/1‑109 (MQ=255)
                                                                                                                                     aGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTta                                                          <  1:405009/109‑1 (MQ=255)
                                                                                                                                      gACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACTGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGa                >  2:1322010/1‑150 (MQ=255)
                                                                                                                                      gACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGa                >  2:209622/1‑150 (MQ=255)
                                                                                                                                       aCGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGAt               <  1:530430/150‑1 (MQ=255)
                                                                                                                                           gTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTcac           >  1:1213289/1‑150 (MQ=255)
                                                                                                                                            tCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTcacc          >  1:920936/1‑150 (MQ=255)
                                                                                                                                            tCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTcacc          <  1:56032/150‑1 (MQ=255)
                                                                                                                                               aGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACt       >  2:1205309/1‑150 (MQ=255)
                                                                                                                                                gCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTg      <  1:814166/150‑1 (MQ=255)
                                                                                                                                                  gTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTca            <  1:1246393/142‑1 (MQ=255)
                                                                                                                                                  gTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGaa    >  2:468174/1‑150 (MQ=255)
                                                                                                                                                   tCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAAc   >  1:606859/1‑150 (MQ=255)
                                                                                                                                                    cGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACg  <  2:1046859/150‑1 (MQ=255)
                                                                                                                                                     gCATCCGGCAAATAGCCCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACg  >  2:1483863/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTGGCATTATTAAAAAAGGCCACCTGAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACG  >  NC_000913/4148354‑4148651

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