Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,439,588 G→C 15.0% intergenic (‑79/+28) yhdN ← / ← rplQ DUF1992 family protein/50S ribosomal subunit protein L17

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,439,5880GC15.0% 516.8 / 73.0 191intergenic (‑79/+28)yhdN/rplQDUF1992 family protein/50S ribosomal subunit protein L17
Reads supporting (aligned to +/- strand):  ref base G (85/76);  new base C (20/9);  total (105/86)
Fisher's exact test for biased strand distribution p-value = 1.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

Only 100 of 214 total aligned reads displayed.

CCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCC  >  NC_000913/3439489‑3439687
                                                                                                   |                                                                                                   
 cACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCgg                                                                                                    >  2:3442085/1‑100 (MQ=255)
 cACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCgg                                                                                                    >  1:180436/1‑100 (MQ=255)
  aCTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACccgccgag                                                                                                   <  1:2849344/100‑3 (MQ=255)
   cTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCgggg                                                                                                  >  2:1326832/1‑100 (MQ=255)
    tGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGc                                                                                                 <  1:762342/100‑1 (MQ=255)
    tGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGc                                                                                                 >  1:1883451/1‑100 (MQ=255)
     ggTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCg                                                                                                <  2:3240722/100‑1 (MQ=255)
       tCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCggg                                                                                              <  1:1434190/100‑1 (MQ=255)
        cAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCggg                                                                                              >  2:3736745/1‑99 (MQ=255)
        cAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGt                                                                                             >  1:981526/1‑100 (MQ=255)
          aGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCggg                                                                                              >  2:2619172/1‑97 (MQ=255)
           gTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGtttt                                                                                          >  1:2792170/1‑100 (MQ=255)
            tAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGttttt                                                                                         <  2:419804/100‑1 (MQ=255)
            tAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGttttt                                                                                         >  1:1120992/1‑100 (MQ=255)
             aaCCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATAGTACGGGTATAAAAAAACCCGCCGGGGCGGGtttttt                                                                                        >  1:2941510/1‑100 (MQ=255)
              aCCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGttttttt                                                                                       <  1:547582/100‑1 (MQ=255)
               ccACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTa                                                                                      <  2:2961363/100‑1 (MQ=255)
                cacaTAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTAc                                                                                     >  1:2154518/1‑100 (MQ=255)
                 acaTAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCCCGGCGGGTTTTTTTata                                                                                    <  1:384160/100‑3 (MQ=255)
                  caTAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGt                                                                                   <  2:64480/100‑1 (MQ=255)
                    tAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTg                                                                                 >  2:4086304/1‑100 (MQ=255)
                          ctacAGTGGATGAACAAAAACAGTACAGCTATTGTCGACAAGGTGGAATACTACGGGTATAACACAACCCTGCGGGTCGGGATTTTTTACGTCGTTTCa                                                                            >  2:2297988/4‑99 (MQ=255)
                          cTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAACAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAg                                                                           <  1:1046521/100‑1 (MQ=255)
                          cTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAg                                                                           >  2:1280902/1‑100 (MQ=255)
                          cTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAg                                                                           >  1:2370990/1‑100 (MQ=255)
                           tCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAACCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGa                                                                          <  1:1442747/100‑1 (MQ=255)
                           tCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGa                                                                          <  2:1605346/100‑1 (MQ=255)
                           tCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGa                                                                          <  2:4116664/100‑1 (MQ=255)
                            ccAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGAt                                                                         <  1:458133/100‑1 (MQ=255)
                            ccAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGATTTTTTACGTTGCTTCAGAt                                                                         >  2:1179327/1‑100 (MQ=255)
                             cAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGAtt                                                                        <  1:2423288/100‑1 (MQ=255)
                               ggggATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTAc                                                                      <  1:2905631/100‑1 (MQ=255)
                               ggggATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTAc                                                                      >  2:3212214/1‑100 (MQ=255)
                                  gATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTActct                                                                   >  2:4043013/1‑100 (MQ=255)
                                    tGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTgc                                                                 <  2:180436/100‑1 (MQ=255)
                                      taCAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAGACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTgcag                                                               <  1:2354794/99‑1 (MQ=255)
                                       aCAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTgcaga                                                              >  2:2354789/1‑99 (MQ=255)
                                          aaaaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCagcagc                                                           >  2:154157/1‑100 (MQ=255)
                                          aaaaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCagcagc                                                           <  2:2268501/100‑1 (MQ=255)
                                          aaaaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCagcagc                                                           >  1:4146431/1‑100 (MQ=255)
                                           aaaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCt                                                          <  2:560376/100‑1 (MQ=255)
                                              gagTACAGCTATTGTAGATACGTGGGGCTACTATGGATATAAAAATACCCACGGTGGCGGGTTTTTTTACGTTGGTGCAGAGTAATGTgaagcaacgtaa                                                       <  2:3629220/100‑7 (MQ=255)
                                              gagTACAGCTATTGTAGATAAGTGGCGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTGCAGAGTAATCTgaagcaacgtaa                                                       <  2:2702589/100‑7 (MQ=37)
                                                   cAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCttt                                                  >  1:3438515/1‑100 (MQ=255)
                                                    aGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTCCTTCAGATTACTCTGCAGCAGCTTCTGCTTtc                                                 <  1:3977312/100‑1 (MQ=255)
                                                    aGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCAGCCGGGGCGGGTGTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTtc                                                 >  1:275231/1‑100 (MQ=255)
                                                     gCTAATGTAGATAAGTGGTGATACTACTGGTAGAAAAATGCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCGGCAGCAGCTTCTGCTTtct                                                <  2:3818416/100‑1 (MQ=255)
                                                        aTTGTAGATAAGTGGGGATACTACGGGTCTAAAAAATCCCGCCGGGGCGGGGTTTTTTACGTTGCTTCAGATTACACTGCAGCAGCTTCTGCTTTgtct                                              >  2:3577554/1‑99 (MQ=255)
                                                        aTTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTg                                             <  1:1193469/100‑1 (MQ=255)
                                                         ttGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGa                                            >  1:2075770/1‑100 (MQ=255)
                                                          tGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGaa                                           <  1:1002671/100‑1 (MQ=255)
                                                           gTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAAc                                          <  1:3736766/100‑1 (MQ=255)
                                                           gTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAAc                                          <  1:4043040/100‑1 (MQ=255)
                                                           gTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAAc                                          <  2:804285/100‑1 (MQ=255)
                                                           gTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAAc                                          <  1:1669480/100‑1 (MQ=255)
                                                           gTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGGTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAAc                                          >  1:286475/1‑100 (MQ=255)
                                                            tAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACg                                         >  1:261019/1‑100 (MQ=255)
                                                            tAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACg                                         <  2:3891272/100‑1 (MQ=255)
                                                             aGATAAGTGGGCATACTACGGGTATAAAAAAACCCGGCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAAc                                          <  1:1961879/98‑1 (MQ=255)
                                                               aTAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATc                                      <  1:2703759/100‑1 (MQ=255)
                                                                 aaGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCaa                                    <  1:846337/100‑1 (MQ=255)
                                                                   gTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTCAAGCAGCTTCTGCTTTCTCTGAACGATCAAcc                                  >  1:934624/1‑100 (MQ=255)
                                                                    tGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCa                                 >  1:3021707/1‑100 (MQ=255)
                                                                     ggggATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAg                                >  1:4190889/1‑100 (MQ=255)
                                                                      gggATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGc                               >  1:699881/1‑100 (MQ=255)
                                                                       ggATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCt                              <  1:2253830/100‑1 (MQ=255)
                                                                        gATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTTAACGATCAACCAGCTc                             <  1:2411056/100‑1 (MQ=255)
                                                                         aTACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCGCTGAACGATCAACCAGCTCg                            <  2:2559057/100‑1 (MQ=255)
                                                                            ctacGGGTATAACAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATg                         <  2:1952080/100‑1 (MQ=255)
                                                                            ctacGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATg                         >  1:35607/1‑100 (MQ=255)
                                                                             aagcaacgtaAAAAAACCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGt                        >  1:1712720/9‑100 (MQ=37)
                                                                                gggTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAg                     <  1:3626245/100‑1 (MQ=255)
                                                                                 ggTATAAAAAAACCCCCCGGGGCGGGCTTTTTTACGTTGCTTCAGATTACCCTGCAGCAGCTTCTGCTTTCTCTGCACGATCAACCAGCTCGATGTaan                     <  2:1010747/99‑2 (MQ=255)
                                                                                  gTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGcc                   >  2:3268450/1‑100 (MQ=255)
                                                                                  gTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGcc                   >  1:1462372/1‑100 (MQ=255)
                                                                                    gtaAAAAAACCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATTTAAGCCAt                 >  2:1555690/2‑100 (MQ=38)
                                                                                    gtaAAAAAACCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCAt                 >  2:3085603/2‑100 (MQ=255)
                                                                                    ataAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCAt                 <  1:4086333/100‑1 (MQ=255)
                                                                                     taAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATc                >  2:1733130/1‑100 (MQ=255)
                                                                                     taAAAAAACCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCGGCAGCAGCTTCTGCTTTCTCTGAACGAGCAACCAGCTCGATGTAAGCCATc                >  2:186352/1‑100 (MQ=38)
                                                                                     taAAAAAACCCGCCCCGGCGGGTTTTTTTAAGTTGCTTTAAATTAATATGCAGCAGCTTTTGGTTTCTTTGAAAGATAAACCAGCTCGAAGTAAGCCAt                 >  2:28004/1‑99 (MQ=16)
                                                                                      aaaaaaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCg               >  2:587422/1‑100 (MQ=255)
                                                                                      aaaaaaaCCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCg               >  1:1319737/1‑100 (MQ=255)
                                                                                       aaaaaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCgg              >  1:3890420/1‑100 (MQ=255)
                                                                                       aaaaaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCgg              >  1:3536133/1‑100 (MQ=255)
                                                                                        aaaaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGgc             <  2:1069905/100‑1 (MQ=255)
                                                                                        aaaaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGgc             >  1:3369768/1‑100 (MQ=255)
                                                                                        aaaaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGgc             <  1:3212228/100‑1 (MQ=255)
                                                                                            aCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGt         <  1:2554875/100‑1 (MQ=255)
                                                                                             cccGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGtt        >  2:247895/1‑100 (MQ=255)
                                                                                             cccGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGtt        >  1:248149/1‑100 (MQ=255)
                                                                                              ccgccgGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTg       >  1:2634018/1‑100 (MQ=255)
                                                                                                gccgGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTc     >  1:3169231/1‑100 (MQ=255)
                                                                                                 ccgGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCGACCAGCTCGATGTAAGCCATCGGCGCGTTGTc     >  2:826889/1‑99 (MQ=255)
                                                                                                 ccgGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCg    <  2:3368715/100‑1 (MQ=255)
                                                                                                 ccgGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCg    <  2:3094476/100‑1 (MQ=255)
                                                                                                  cgGGGCGGGTTTTTTTACGTTGCTTCAGATTACCCGTCAGCAGGTTCTGCCGTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGt         <  2:2947148/94‑1 (MQ=255)
                                                                                                   ggggCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGcc  >  1:3593967/1‑100 (MQ=255)
                                                                                                   ggggCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGcc  <  2:1360065/100‑1 (MQ=255)
                                                                                                   ggggCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGcc  <  2:1497667/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
CCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCC  >  NC_000913/3439489‑3439687

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: