Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,294,680 Δ1 bp 5.6% coding (225/1944 nt) ccmF ← heme lyase, CcmF subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,294,6800G. 5.6% 237.0 / 7.1 72coding (225/1944 nt)ccmFheme lyase, CcmF subunit
Reads supporting (aligned to +/- strand):  ref base G (35/32);  new base . (2/2);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.15e-01

CGGGCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGAAAGCATTGACCAGCACCAGAAATGCGCCAGCCACAGACATAAACAGCAGCCAGGCAAACAAGCGGGAAGACGCCATCATGCGCGCATCTCCGCGCGCC  >  NC_000913/2294534‑2294821
                                                                                                                                                  |                                                                                                                                             
cGGGCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGa                                                                                                                                             <  3:205151/149‑1 (MQ=255)
 gggCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCCCGCGATACCACACCGGAAGCTGGGTATTGGAg                                                                                                                                            >  7:248542/1‑149 (MQ=255)
        aaTATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTgg                                                                                                                                     >  8:341931/1‑149 (MQ=255)
          tatCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGc                                                                                                                                    >  8:98463/1‑148 (MQ=255)
               aGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACa                                                                                                                                >  7:152821/1‑147 (MQ=255)
                  ggAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGt                                                                                                                           >  5:240106/1‑149 (MQ=255)
                       aCGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGAGTTGCTGGCAACATAGGTGACGGt                                                                                                                     >  6:212925/1‑149 (MQ=255)
                       aCGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGAGTTGCTGGCAACATAGGTGACGGt                                                                                                                     >  6:207531/1‑149 (MQ=255)
                         gCTGACTAAAAATAGCCCCCGCAAAGGTCCAGCCGCTCATCAGCAGCCCCCCCAGCAGTAGCGAGCCTTAATGCGCGCACCAGGTAGCCGCCCCGCGATACCACACCGGAAGCTGGGTATGGGAGTTGATTGCAACATAGGTGACgtt                                                                                                                     >  6:312560/1‑146 (MQ=255)
                           tGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGACGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTg                                                                                                                    >  2:304933/1‑147 (MQ=255)
                             aCTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAg                                                                                                                 <  6:240106/148‑1 (MQ=255)
                             aCTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGt                                                                                                                >  1:66217/1‑149 (MQ=255)
                             aCTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGt                                                                                                                >  2:255915/1‑149 (MQ=255)
                             aCTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGt                                                                                                                >  3:314422/1‑149 (MQ=255)
                               tAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCg                                                                                                              >  7:46769/1‑149 (MQ=255)
                                    aTCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGAc                                                                                                         <  8:152821/149‑1 (MQ=255)
                                       gCCACCGCAAAGGTCCAGCCGCTCTTCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATAACACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCAc                                                                                                      >  1:278688/1‑149 (MQ=255)
                                       gCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCAc                                                                                                      <  7:341931/149‑1 (MQ=255)
                                        ccACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCaa                                                                                                                                  <  1:134676/120‑1 (MQ=255)
                                        ccACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCaa                                                                                                                                  >  2:134676/1‑120 (MQ=255)
                                        ccACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCAc                                                                                                      >  4:81415/1‑148 (MQ=255)
                                         cACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGa                                                                                                    <  4:314422/149‑1 (MQ=255)
                                             gCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGAAAg                                                                                                 >  4:105745/1‑148 (MQ=255)
                                             gCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGAAAGc                                                                                                <  8:327883/149‑1 (MQ=255)
                                                 aGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGAAAGCATTg                                                                                            >  5:282967/1‑149 (MQ=255)
                                                   gTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGAAAGCATTGAc                                                                                          <  1:362971/149‑1 (MQ=255)
                                                   gTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGAAAGCATTGAc                                                                                          <  8:1214/149‑1 (MQ=255)
                                                   gTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGAAAGCATTGAc                                                                                          <  3:139388/149‑1 (MQ=255)
                                                     ccAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGAAAGCATTGACCa                                                                                        <  2:186195/149‑1 (MQ=255)
                                                        gCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCCCGCGATACCACACCGGAAGCGGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCCCGAAAGCATTGACCAGCa                                                                                     >  1:232480/1‑149 (MQ=255)
                                                         ccGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGAAAGCATTGACCAGCAc                                                                                    <  3:105745/149‑1 (MQ=255)
                                                          cGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGAAAGCATTGACCAGCAc                                                                                    >  3:328188/1‑148 (MQ=255)
                                                          cGCTCATCAACAGCACCCACAGCAGTAGCGAGCCATCCTGCGCGACCCAGCCAGCCGCCACGCGATACAACACCGGCAGCGGGGAATAGGAGTTGCAGGCAACATAGGGGACGGTGAATTCGTTGACCACGAAACCATTGACCAGcccc                                                                                   <  1:318705/149‑1 (MQ=255)
                                                            cTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGAAAGCATTGACCAGCACCAg                                                                                 <  8:144285/149‑1 (MQ=255)
                                                               atcaGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGAAAGCATTGACCAGCACCAGaaa                                                                              <  5:207531/148‑1 (MQ=255)
                                                               atcaGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGAAAGCATTGACCAGCACCAGaaa                                                                              <  5:212925/148‑1 (MQ=255)
                                                                 cagcagCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTg                                                                                                           <  5:96811/118‑1 (MQ=255)
                                                                 cagcagCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTg                                                                                                           >  6:96811/1‑118 (MQ=255)
                                                                    cagcaCCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGAAAGCATTGACCAGCACCAGAAATgcgc                                                                         <  5:132895/149‑1 (MQ=255)
                                                                      gcaCCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGACAGCATTGCCCAGCACCAGAAATgcgc                                                                         >  1:168325/1‑147 (MQ=255)
                                                                      gcaCCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGAAAGCATTGACCAGCACCAGAAATGCgcca                                                                       <  1:255915/149‑1 (MQ=255)
                                                                        aCCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGAAAGCATTGACCAGCACCAGAAATGCgccagc                                                                     >  3:175575/1‑149 (MQ=255)
                                                                         cccACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGAAAGCATTGAcc                                                                                         >  4:218246/1‑128 (MQ=255)
                                                                         cccACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGAAAGCATTGAcc                                                                                         <  3:218246/128‑1 (MQ=255)
                                                                          ccACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGAAAGCATTGACCAGCACCAGAAATGCgccagcca                                                                   >  1:260832/1‑149 (MQ=255)
                                                                            acaGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGAAAGCATTGACCAGCACCAGAAATGCgccagcca                                                                   >  8:56008/1‑147 (MQ=255)
                                                                            acaGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGAAAGCATTGACCAGCACCAGAAATGCGCCAGCcaca                                                                 <  1:304933/149‑1 (MQ=255)
                                                                            acaGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGAAAGCATTGACCAGCACCAGAAATGCGCCAGCcaca                                                                 <  3:81415/149‑1 (MQ=255)
                                                                                cagGAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGGCACGGGATACCACACGGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGGGACGGTGAAGTCGTTGAGCACGAAAGCATTGACCAGCTCCAGAAATGCGCCAGCCACAGACa                                                             <  4:348911/149‑1 (MQ=255)
                                                                                 agTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGAAAGCATTGACCAGCACCAGAAATGCGCCAGCCACAGACAt                                                            <  1:128208/149‑1 (MQ=255)
                                                                                  gTAGCGAGCCTTCATGCGCGCCCCTGGTAGCCGCCACGCGATACGACCCCGGAAGCTGGGTTTTGGAGTTGCTGGCAACATAGGTGAGGGTGAAGTCGTTGAGGACGAAAGCATTGACCAGCACCAGAAATGCGCCAGCCCGAGACAt                                                            >  2:120991/1‑148 (MQ=255)
                                                                                     gCGAGCCTTCATGCGCGGCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCATGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGAAAGCATTGACCAGCACCAGAAATGCGCCAGCCACAGACATaaa                                                         <  3:332015/148‑1 (MQ=255)
                                                                                        aGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGAAAGCATTGACCAGCACCAGAAATGCGCCAGCCACAGACATAAAcag                                                      >  8:287570/1‑148 (MQ=255)
                                                                                         gCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGAAAGCATTGACCAGCACCAGAAATGCGCCAGCCACAGACATAAAcagca                                                    <  8:248542/149‑1 (MQ=255)
                                                                                            ttCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGAAAGCATTGACCAGCACCAGAAATGCGCCAGCCACAGACATAAACagcagc                                                  >  3:40854/1‑148 (MQ=255)
                                                                                               aTGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGAAAGCATTGACCAGCACCAGAAATGCGCCAGCCACAGACATAAACAGCAGCCAg                                               <  6:282967/148‑1 (MQ=255)
                                                                                               aTGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACACGAAAGCATTGACCAGCACCAGAAATGCGCCAGCCACAGACATAAACAGCAGCCAGGc                                             <  2:66217/149‑1 (MQ=255)
                                                                                               aTGCGCGCCCAAGGTAGCCGCCACGCGATACCACACAGGAAGATGGGTATTGGAGTTGCTGGCAACATAGTTGACGGTGAAGGCGTTGAAGACGAAAGCATTGATCATAACCAGAATTTCGCAAGCCACAGACATAAACATCAGCCAg                                               <  8:347519/148‑1 (MQ=255)
                                                                                                    cgcCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTTGAGTTGATGGCAACATAGGTGACGGTGAAGTCGTTGACCACGAAAGCATTGACCAGCACCAGAAATGCGCCAGCCACAGACATAAACAGCAGCCAGGCAAAc                                         <  6:20621/149‑1 (MQ=255)
                                                                                                        ccAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGAAAGCATTGACCAGCACCAGAAATGCGCCAGCCACAGACATAAACAGCAGCCAGGCAAACAAg                                      >  8:240797/1‑148 (MQ=255)
                                                                                                        ccAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGAAAGCATTGACCAGCACCAGAAATGCGCCAGCCACAGACATAAACAGCAGCCAGGCAAACAAg                                      <  7:240797/148‑1 (MQ=255)
                                                                                                         cAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTTAAGTCGTTGACCACGAAAGCATTGACCAGCACCAGAAATGCGCCAGCCACAGACATAAACAGCAGCCAGGCAAACAAGCg                                    <  4:40854/149‑1 (MQ=255)
                                                                                                         cAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGAAAGCATTGACCAGCACCAGAAATGCGCCAGCCACAGACATAAACAGCAGCCAGGCAAACAAGCg                                    <  8:46769/149‑1 (MQ=255)
                                                                                                             tAGCCGCCACGCGATACCACACCGGAAGCTGGGTAGTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGAAAGCATTGACCAGCACCAGAAATGCGCCAGCCACAGACATAAACAGCAGCCAGGCAAACAAGCGGGaa                                <  3:286232/149‑1 (MQ=255)
                                                                                                                         gATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGAAAGCATTGACCAGCACCAGAAATGCGCCAGCCACAGACATAAACAGCAGCCAGGCAAACAAGCGGGAAGACGCCATCATg                    >  7:315654/1‑149 (MQ=255)
                                                                                                                              cacaCCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGAAAGCATTGACCAGCACCAGAAATGCGCCAGCCACAGACATAAACAGCAGCCAGGCAAACAAGCGGGAAGACGCCATCATgcgca                >  3:220857/1‑147 (MQ=255)
                                                                                                                                 acCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGAAAGCATTGAc                                                                                          <  5:115309/71‑1 (MQ=255)
                                                                                                                                 acCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGAAAGCATTGAc                                                                                          >  6:115309/1‑71 (MQ=255)
                                                                                                                                   cGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGAAAGCATTGACCAGCACCAGAAATGCGCCAGCCACAGACATAAACAGCAGCCAGGCAAACAAGCGGGAAGACGCCATCATGCGCGCAtctc          >  3:184080/1‑149 (MQ=255)
                                                                                                                                   cGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGAAAGCATTGACCAGCACCAGAAATGCGCCAGCCACAGACATAAACAGCAGCCAGGCAAACAAGCGGGAAGACGCCATCATGCGCGCAtctc          >  7:60207/1‑149 (MQ=255)
                                                                                                                                       aGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGAAAGCATTGACCAGCACCAGAAATGCGCCAGCCACAGACATAAACAGCAGCCAGGCAAACAAGCGGGAAGACGCCATCATGCGCGCATCTCcgcg      >  1:49284/1‑149 (MQ=255)
                                                                                                                                           gggTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGAAAGCATTGACCAGCACCAGAAATGCGCCAGCCACAGCCATAAACAGCAGCCAGGCAAACAAGCGGGAAGACGCCATCATGCGCGCATCTCCGCGCGcc  >  6:376942/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
CGGGCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGAAAGCATTGACCAGCACCAGAAATGCGCCAGCCACAGACATAAACAGCAGCCAGGCAAACAAGCGGGAAGACGCCATCATGCGCGCATCTCCGCGCGCC  >  NC_000913/2294534‑2294821

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