Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,999,755 (C)6→7 coding (137/756 nt) ygeR ← LysM domain‑containing putative peptidase lipoprotein YgeR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,999,7491.C97.0% 807.4 / 14.4 200coding (143/756 nt)ygeRLysM domain‑containing putative peptidase lipoprotein YgeR
Reads supporting (aligned to +/- strand):  ref base . (3/3);  new base C (110/84);  total (113/87)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 226 total aligned reads displayed.

TTTCTGACCAACTTCAATGGTGTAAGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTGCTACCCGAACAGCCCGCCAAAAGCAGTCCAACCGATAACAACA  >  NC_000913/2999649‑2999848
                                                                                                     |                                                                                                   
tttCTGACCAACTTCAATGGTGTAAGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTAt                                                                                                      <  2:4561158/101‑1 (MQ=255)
 ttCTGACCAACTTCAATGGTGTAAGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATc                                                                                                    <  2:4561322/101‑1 (MQ=255)
 ttCTGACCAACTTCAATGGTGTAAGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATc                                                                                                    >  1:4561558/1‑101 (MQ=255)
  tCTGACCAACTTCAATGGTGTAAGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATcc                                                                                                   <  1:4561170/101‑1 (MQ=255)
   cTGACCAACTTCAATGGTGTAAGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATccc                                                                                                  >  1:4561560/1‑101 (MQ=255)
       ccAACTTCAATGGTGTAAGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATccccccc                                                                                              >  1:4561564/1‑100 (MQ=255)
       ccAACTTCAATGGTGTAAGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATccccccc                                                                                              >  2:4561563/1‑100 (MQ=255)
        cAACTTCAATGGTGTAAGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATcccccccg                                                                                             <  1:4561437/101‑3 (MQ=255)
         aaCTTCAATGGTGTAAGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATcccccccgt                                                                                            <  1:4561120/101‑4 (MQ=255)
         aaCTTCAATGGTGTAAGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATcccccccgt                                                                                            >  2:4561568/1‑98 (MQ=255)
         aaCTTCAATGGTGTAAGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATcccccccgt                                                                                            >  1:4561566/1‑98 (MQ=255)
            ttCAATGGTGTAAGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTc                                                                                          >  1:4561574/1‑101 (MQ=255)
            ttCAATGGTGTAAGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTc                                                                                          >  1:4561572/1‑101 (MQ=255)
             tCAATGGTGTAAGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCa                                                                                         <  2:4561103/101‑1 (MQ=255)
             tCAATGGTGTAAGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCa                                                                                         >  1:4561577/1‑101 (MQ=255)
              cAATGGTGTAAGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCAc                                                                                        <  2:4561189/101‑1 (MQ=255)
              cAATGGTGTAAGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCAc                                                                                        >  2:4561580/1‑101 (MQ=255)
              cAATGGTGTAAGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCAc                                                                                        <  1:4561124/101‑1 (MQ=255)
               aaTGGTGTAAGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACg                                                                                       <  2:4561168/101‑1 (MQ=255)
                aTGGTGTAAGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACgg                                                                                      >  1:4561582/1‑101 (MQ=255)
                 tGGTGTAAGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGgt                                                                                     <  1:4561305/101‑1 (MQ=255)
                     gtAAGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTaa                                                                                 <  1:4561145/101‑1 (MQ=255)
                     gtAAGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTaa                                                                                 >  1:4561585/1‑101 (MQ=255)
                      tAAGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTaaa                                                                                <  2:4561180/101‑1 (MQ=255)
                      tAAGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTaaa                                                                                >  2:4561586/1‑101 (MQ=255)
                      tAAGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTGTAAAc                                                                               <  1:4561141/101‑1 (MQ=255)
                       aaGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAGATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAAc                                                                               >  1:4561590/1‑101 (MQ=255)
                        aGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACg                                                                              >  2:4561591/1‑101 (MQ=255)
                          ggggggAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGa                                                                            <  1:4561232/101‑1 (MQ=255)
                          ggggggAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGa                                                                            >  2:4561593/1‑101 (MQ=255)
                            ggggAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGc                                                                          >  1:4561596/1‑101 (MQ=255)
                            ggggAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGc                                                                          <  2:4561139/101‑1 (MQ=255)
                            ggggAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGc                                                                          >  2:4561597/1‑101 (MQ=255)
                              ggAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCg                                                                        >  1:4561601/1‑101 (MQ=255)
                               gAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCgg                                                                       <  1:4561354/101‑1 (MQ=255)
                               gAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCgg                                                                       >  2:4561602/1‑101 (MQ=255)
                                 aaTGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGaa                                                                     <  2:4561222/101‑1 (MQ=255)
                                  aTGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAAt                                                                    >  1:4561605/1‑101 (MQ=255)
                                  aTGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAAt                                                                    >  2:4561606/1‑101 (MQ=255)
                                   tGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATa                                                                   <  1:4561155/101‑1 (MQ=255)
                                     ccGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACg                                                                 <  2:4561042/101‑1 (MQ=255)
                                      cGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGt                                                                >  1:4561611/1‑101 (MQ=255)
                                      cGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGt                                                                >  1:4561609/1‑101 (MQ=255)
                                       gTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGtt                                                               <  1:4561478/101‑1 (MQ=255)
                                        ttCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTc                                                              >  2:4561614/1‑101 (MQ=255)
                                         tCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTcc                                                             >  1:4561615/1‑101 (MQ=255)
                                           aGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTg                                                           >  1:4561619/1‑101 (MQ=255)
                                           aGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTg                                                           >  2:4561618/1‑101 (MQ=255)
                                            gTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGt                                                          <  2:4561241/101‑1 (MQ=255)
                                            gTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGt                                                          <  2:4561547/101‑1 (MQ=255)
                                              cgcgCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTAt                                                        <  1:4561336/101‑1 (MQ=255)
                                               gcgcCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATc                                                       <  2:4561271/101‑1 (MQ=255)
                                               gcgcCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATc                                                       <  2:4561095/101‑1 (MQ=255)
                                                cgcCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCg                                                      <  1:4561053/101‑1 (MQ=255)
                                                cgcCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCg                                                      >  1:4561624/1‑101 (MQ=255)
                                                cgcCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCg                                                      <  1:4561244/101‑1 (MQ=255)
                                                  ccAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGga                                                    >  2:4561625/1‑101 (MQ=255)
                                                   cAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGgat                                                   <  2:4561211/101‑1 (MQ=255)
                                                    aGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGgatg                                                  >  1:4561628/1‑101 (MQ=255)
                                                     gTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGgatga                                                 >  1:4561631/1‑101 (MQ=255)
                                                     gTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGgatga                                                 >  2:4561630/1‑101 (MQ=255)
                                                       tCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGAtt                                               >  2:4561633/1‑101 (MQ=255)
                                                        cTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGAttt                                              <  2:4561255/101‑1 (MQ=255)
                                                         ttttACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTg                                             >  1:4561635/1‑101 (MQ=255)
                                                          tttACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGgatg                                                  >  2:4561640/1‑95 (MQ=255)
                                                          tttACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTGc                                            >  1:4561639/1‑101 (MQ=255)
                                                          tttACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTGc                                            <  2:4561280/101‑1 (MQ=255)
                                                           ttACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTGCt                                           <  1:4561262/101‑1 (MQ=255)
                                                               gCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTGCTAccc                                       >  2:4561644/1‑101 (MQ=255)
                                                                cTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTGCTACCCg                                      >  1:4561645/1‑101 (MQ=255)
                                                                 tGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTGCTACCCGa                                     >  1:4561647/1‑101 (MQ=255)
                                                                 tGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTGCTACCCGa                                     >  1:4561649/1‑101 (MQ=255)
                                                                  ggTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTTCTAccc                                       >  1:4561653/1‑98 (MQ=255)
                                                                   gTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTGCTACCCGAAc                                   >  2:4561654/1‑101 (MQ=255)
                                                                     tCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTGCTACCCGAACAg                                 >  2:4561657/1‑101 (MQ=255)
                                                                     tCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTGCTACCCGAACAg                                 <  2:4561327/101‑1 (MQ=255)
                                                                     tCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTGCTACCCGAACAg                                 <  2:4561237/101‑1 (MQ=255)
                                                                     tCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTGCTACCCGAACAg                                 <  1:4561182/101‑1 (MQ=255)
                                                                        cGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTGCTACCCGAACAGccc                              >  2:4561659/1‑101 (MQ=255)
                                                                        cGTGGTGCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTGCTACCCGAACAGccc                              <  2:4561226/101‑1 (MQ=255)
                                                                         gtggtgCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTGCTACCCGAACAGCCCg                             <  2:4561464/101‑1 (MQ=255)
                                                                         gtggtgCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTGCTACCCGAACAGCCCg                             >  1:4561663/1‑101 (MQ=255)
                                                                         gtggtgCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTGCTACCCGAACAGCCCg                             >  1:4561660/1‑101 (MQ=255)
                                                                            gtgCGCGAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTGCTACCCGAACAGCCCGCCa                          >  1:4561666/1‑101 (MQ=255)
                                                                               cgcgAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCgg                                                     >  1:4561669/1‑71 (MQ=255)
                                                                                gcgAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTGCTCCCCGAACAGCCCGCCAAAAg                      <  2:4561022/101‑1 (MQ=255)
                                                                                  gAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTGCTACCCGAACAGCCCGCCAAAAGCa                    >  2:4561670/1‑101 (MQ=255)
                                                                                  gAAATACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTGCTACCCGAACAGCCCGCCAAAAGCa                    <  2:4561062/101‑1 (MQ=255)
                                                                                   aaaTACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTGCTACCCGAACAGCCCGCCAAAAGCAg                   <  2:4561215/101‑1 (MQ=255)
                                                                                    aaTACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTGCTACCCGAACAGCCCGCCAAAAGCAGt                  >  1:4561672/1‑101 (MQ=255)
                                                                                    aaTACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTGCTACCCGAACAGCCCGCCAAAAGCAGt                  <  1:4561104/101‑1 (MQ=255)
                                                                                     aTACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCTGATGATTTGCt                                           >  2:4561675/1‑75 (MQ=255)
                                                                                     aTACGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTGCTACCCGAACAGCCCGCCAAAAGCAGTc                 <  1:4561380/101‑1 (MQ=255)
                                                                                       aCGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTGCTACCCGAACAGCCCGCCAAAAGCAGTCCa               >  1:4561676/1‑101 (MQ=255)
                                                                                        cGATATAGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTGCTACCCGAACAGCCCGCCAAAAGCAGTCCaa              <  2:4561235/101‑1 (MQ=255)
                                                                                           tataGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTGCTACCCGAACAGCCCGCCAAAAGCAGTCCAACCg           >  2:4561678/1‑101 (MQ=255)
                                                                                            ataGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTGCTACCCGAACAGCCCGCCAAAAGCAGTCCAACCGa          <  2:4561302/101‑1 (MQ=255)
                                                                                            ataGCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTGCTACCCGAACAGCCCGCCAAAAGCAGTCCAACCGa          >  1:4561682/1‑101 (MQ=255)
                                                                                               gCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTGCTACCCGAACAGCCCGCCAAAAGCAGTCCAACCGATaa       <  2:4561294/101‑1 (MQ=255)
                                                                                               gCGTATCCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTGCTACCCGAACAGCCCGCCAAAAGCAGTCCAACCGATaa       >  1:4561685/1‑101 (MQ=255)
                                                                                                     |                                                                                                   
TTTCTGACCAACTTCAATGGTGTAAGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTCCTGTATCGGATGATTTGCTACCCGAACAGCCCGCCAAAAGCAGTCCAACCGATAACAACA  >  NC_000913/2999649‑2999848

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