Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,145,433 C→T 7.3% intergenic (‑100/+41) hemO ← / ← PP_1006 heme oxygenase/heme receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,145,4330CT 7.3% 196.9 / 4.5 68intergenic (‑100/+41)hemO/PP_1006heme oxygenase/heme receptor
Reads supporting (aligned to +/- strand):  ref base C (32/31);  new base T (3/2);  total (35/33)
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.58e-01

GTGATCTGGTTCAGGCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTG  >  NC_002947/1145287‑1145579
                                                                                                                                                  |                                                                                                                                                  
gTGATCTGGTTCAGGCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTcgc                                                                                                                                                  <  2:165665/149‑1 (MQ=255)
        gTTCAGGCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGc                                                                                                                                          >  8:198866/1‑149 (MQ=255)
        gTTCAGGCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAACCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGc                                                                                                                                          >  3:41829/1‑149 (MQ=255)
               cgcTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCaaa                                                                                                                                   <  2:235483/149‑1 (MQ=255)
                    ggAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAagag                                                                                                                                <  5:146521/147‑1 (MQ=255)
                      aGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTTGCGGGCACGCCCGCGAAGAGGCCg                                                                                                                            <  8:257662/149‑1 (MQ=255)
                           aGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAaga                                                                                                                                 >  1:186381/1‑139 (MQ=255)
                           aGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAaga                                                                                                                                 <  2:186381/139‑1 (MQ=255)
                           aGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTa                                                                                                                         >  1:175682/1‑147 (MQ=255)
                                     ggtcggtcATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGc                                                                                                             >  3:67593/1‑149 (MQ=255)
                                         ggtcATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTc                                                                                                          <  2:167065/148‑1 (MQ=255)
                                         gggcATGGGGGAGTCCTTGGAAACTGGGCGTCTATATACAAGACGAAACAGCGGGTTCGATTCATTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTc                                                                                                          <  6:345164/145‑1 (MQ=255)
                                               gggggAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATc                                                                                                   >  2:221600/1‑149 (MQ=255)
                                                ggggAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATc                                                                                                   <  6:17659/148‑1 (MQ=255)
                                                 gggAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATccc                                                                                                 <  7:198866/149‑1 (MQ=255)
                                                      tCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAAc                                                                                             >  8:60130/1‑148 (MQ=255)
                                                         ttGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACcaca                                                                                         >  1:197480/1‑149 (MQ=255)
                                                         ttGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCGCAGATATGCCAAACcaca                                                                                         >  3:49/1‑149 (MQ=255)
                                                         ttGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGACGGTACAGACGTAACGCCTCAGATATCCCAAACcaca                                                                                         >  1:202691/1‑149 (MQ=255)
                                                         ttGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCAAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTAAAGACGTAACGCATCAGATATCCCAAACcac                                                                                          >  5:212444/1‑148 (MQ=255)
                                                              aaGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGa                                                                                     >  5:336560/1‑148 (MQ=255)
                                                                   ggCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCg                                                                                <  1:221600/148‑1 (MQ=255)
                                                                      gTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACAtt                                                                                       >  2:214949/1‑138 (MQ=255)
                                                                      gTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACAtt                                                                                       <  1:214949/138‑1 (MQ=255)
                                                                       tCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGTGCGTTCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACAGAAGtt                                                                           >  1:218295/1‑149 (MQ=255)
                                                                            aTACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTTGCGGGCACGCCCGCGAAGAGGCCGGTAct                                                                                                                       <  8:197237/100‑2 (MQ=25)
                                                                            aTACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTTGCGGGCACGCCCGCGAAGAGGCCGGTAct                                                                                                                       >  7:197237/1‑99 (MQ=25)
                                                                              aCAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACTGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGcccc                                                                    >  8:67310/1‑149 (MQ=255)
                                                                              aCAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGcccc                                                                    >  4:143221/1‑149 (MQ=255)
                                                                                aaGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCgg                                                                  >  4:51946/1‑149 (MQ=255)
                                                                                aaGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCgg                                                                  >  3:186719/1‑149 (MQ=255)
                                                                                 aGACGAAACAGCGGGATCGAGACAGTAAAAAAAACCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGc                                                                 <  4:41829/149‑1 (MQ=255)
                                                                                     gAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATgg                                                             <  3:143221/149‑1 (MQ=255)
                                                                                       aaCAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGCTACCGAAGTTGCGCCCCGGCATGGt                                                            >  6:221492/1‑148 (MQ=255)
                                                                                         cAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCg                                                                   >  3:289208/1‑139 (MQ=255)
                                                                                         cAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCg                                                                   <  4:289208/139‑1 (MQ=255)
                                                                                          aGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGGAAAGAGGCCTGTAGAGACGTATCTCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAgg                                                        <  5:241137/149‑1 (MQ=255)
                                                                                               gAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTc                                                   <  4:186719/149‑1 (MQ=255)
                                                                                                  cGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGAt                                                <  2:175682/149‑1 (MQ=255)
                                                                                                   gagaCAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATg                                               >  5:134043/1‑149 (MQ=255)
                                                                                                    agaCAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATAACCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGt                                              >  4:154448/1‑149 (MQ=255)
                                                                                                     gaCAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGtt                                             >  3:155127/1‑149 (MQ=255)
                                                                                                      aCAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTtg                                            <  3:51946/149‑1 (MQ=255)
                                                                                                       cAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCTGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTtg                                            <  2:197480/148‑1 (MQ=255)
                                                                                                       cAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTtgg                                           <  6:336560/149‑1 (MQ=255)
                                                                                                       cAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAAAGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTtg                                            <  2:202691/148‑1 (MQ=255)
                                                                                                         gTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTgc                                                                         >  4:36809/1‑117 (MQ=255)
                                                                                                         gTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTgc                                                                         <  3:36809/117‑1 (MQ=255)
                                                                                                              aaaaGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGgctggctg                                    <  3:154448/149‑1 (MQ=255)
                                                                                                               aaaGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGctggctgg                                   <  4:67593/149‑1 (MQ=255)
                                                                                                                 aGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGt                                 <  6:134043/149‑1 (MQ=255)
                                                                                                                     ggTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTca                               >  1:243938/1‑147 (MQ=255)
                                                                                                                              gTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGTGACCCAgccagccaa                    <  5:44605/149‑7 (MQ=255)
                                                                                                                                  ccAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACAt                                                                                        <  7:252881/77‑1 (MQ=255)
                                                                                                                                  ccAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACAt                                                                                        >  8:252881/1‑77 (MQ=255)
                                                                                                                                   cAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTggctggct                                     <  2:139640/127‑1 (MQ=255)
                                                                                                                                   cAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTggctggct                                     >  1:139640/1‑127 (MQ=255)
                                                                                                                                   cAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTg               <  7:60130/149‑1 (MQ=255)
                                                                                                                                    aGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGGTGGCTGGCTGGGTCACAGCCGCCTCACCCTgg              <  7:53023/149‑1 (MQ=255)
                                                                                                                                       aCCGGCCTCTTTGCGGGCACGCCCGCGAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCcg           >  6:351841/1‑149 (MQ=255)
                                                                                                                                         cGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTcac                              >  2:278274/1‑128 (MQ=255)
                                                                                                                                         cGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTcac                              <  1:278274/128‑1 (MQ=255)
                                                                                                                                         cGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCcgtc         >  1:231805/1‑149 (MQ=255)
                                                                                                                                            ccTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCcgtcg        >  5:217764/1‑147 (MQ=255)
                                                                                                                                            ccTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTa      <  4:155127/149‑1 (MQ=255)
                                                                                                                                              tctTTGCGGGCACGCCCGCGAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGc    >  7:263976/1‑149 (MQ=255)
                                                                                                                                              tctTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCcg           <  4:59184/142‑1 (MQ=255)
                                                                                                                                              tctTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCcg           >  3:59184/1‑142 (MQ=255)
                                                                                                                                                ttCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGcaa  >  1:124063/1‑147 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GTGATCTGGTTCAGGCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTG  >  NC_002947/1145287‑1145579

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