Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,349,081 C→A 5.7% intergenic (‑402/+72) PP_3820 ← / ← galU group II intron‑encoding maturase/UTP‑glucose‑1‑phosphate uridylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,349,0810CA 5.7% 221.8 / 2.6 72intergenic (‑402/+72)PP_3820/galUgroup II intron‑encoding maturase/UTP‑glucose‑1‑phosphate uridylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (34/32);  new base A (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 = 7.12e-01

GCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGAAGTTGGTCGCGTCGATGTAGCCTTCGGCGCCAC  >  NC_002947/4348934‑4349226
                                                                                                                                                   |                                                                                                                                                 
gCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTAcgc                                                                                                                                                 <  1:658744/150‑1 (MQ=32)
             cTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATg                                                                                                                                    >  2:345112/1‑150 (MQ=32)
              tCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGc                                                                                                                                   <  2:1175170/150‑1 (MQ=32)
              tCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCGTAAGCGCAAccc                                                                                                                                          <  1:578245/143‑3 (MQ=12)
              tCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCGTAAGCGCAAccc                                                                                                                                          >  2:578245/1‑141 (MQ=12)
                        ttcAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCCTGCCTGGCGcaca                                                                                                                         >  2:426482/1‑150 (MQ=17)
                        ttcAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGcaca                                                                                                                         >  1:1202592/1‑150 (MQ=32)
                        ttcAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGcaca                                                                                                                         >  2:232024/1‑150 (MQ=32)
                                 ccgatctcccGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGc                                                                                                                <  1:1075823/143‑1 (MQ=17)
                                      cACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGGGTATCAGGCGTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGAGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCAccc                                                                                                           >  1:688414/1‑150 (MQ=11)
                                       aCCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACCTGGCCACCACAACCACGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCAcccc                                                                                                          <  2:856302/150‑1 (MQ=25)
                                         ccGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGa                                                                                                        >  1:562814/1‑150 (MQ=38)
                                              acacCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGAt                                                                                                   <  1:201969/150‑1 (MQ=255)
                                              acacCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGAt                                                                                                   <  2:641858/150‑1 (MQ=255)
                                                 cccTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGc                                                                                                >  1:847481/1‑150 (MQ=255)
                                                  ccTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCt                                                                                               <  2:1202592/150‑1 (MQ=255)
                                                      gCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGa                                                                                           <  2:562814/150‑1 (MQ=255)
                                                       ccTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAg                                                                                          <  1:359156/150‑1 (MQ=255)
                                                         tCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGc                                                                                        >  2:661783/1‑150 (MQ=255)
                                                                  cACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGc                                                                               >  1:960395/1‑150 (MQ=255)
                                                                    cTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCtt                                                                             >  1:167677/1‑150 (MQ=255)
                                                                      tCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTTAGGCCAGCAGGGCTTaa                                                                           >  2:970777/1‑149 (MQ=255)
                                                                      tCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTAt                                                                           <  2:511360/150‑1 (MQ=255)
                                                                       ccccTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATc                                                                          >  2:639760/1‑150 (MQ=255)
                                                                         ccTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAggg                                                                                <  2:254079/142‑1 (MQ=255)
                                                                         ccTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAggg                                                                                >  1:254079/1‑142 (MQ=255)
                                                                          cTTGCCGGGTGTGTAGAGGACTTTGACCTCCAAGTCACCAACTTGGACACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCACCATTCCTGGCGTACACAAAAAAGCCACACCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGt                                                                       >  2:440092/1‑150 (MQ=18)
                                                                           ttGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCAGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTa                                                                      >  2:1194676/1‑150 (MQ=255)
                                                                            tGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTaa                                                                     <  2:847481/150‑1 (MQ=255)
                                                                               cGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGcc                                                                  <  1:661783/150‑1 (MQ=255)
                                                                                gggTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCt                                                                 >  2:195253/1‑150 (MQ=255)
                                                                                   tgtgtAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGc                                                              <  1:232024/150‑1 (MQ=255)
                                                                                       tAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAggg                                                                                <  1:545807/128‑1 (MQ=255)
                                                                                       tAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAggg                                                                                >  2:545807/1‑128 (MQ=255)
                                                                                       tAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGt                                                          >  1:733400/1‑150 (MQ=255)
                                                                                        agagGACTTTCACCTCCAAGTCACCCACTTGGCCACCACACCCAAGCGGGTTGCGCTTCGGCTCACCGCGCCATGCCTGGAGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTACTCAGTAAGCCTTGCCAGTc                                                         >  2:548490/1‑150 (MQ=21)
                                                                                        agagGACTCTCACCTCCAAGTCACCCAGTTGGCCACCACACCCAAGCGGGGGGGGGGGGGGGGCAAGGCGCCATGCCCGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTc                                                         <  2:167677/150‑1 (MQ=21)
                                                                                             aCTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTgta                                                    <  2:463371/150‑1 (MQ=255)
                                                                                              cTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTgtag                                                   <  1:782391/150‑1 (MQ=255)
                                                                                               tttCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTc                                                         >  1:973870/1‑143 (MQ=255)
                                                                                               tttCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTc                                                         <  2:973870/143‑1 (MQ=255)
                                                                                                ttCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTgtagta                                                 <  1:639760/150‑1 (MQ=255)
                                                                                                  cACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTagtagt                                               >  2:1090585/1‑150 (MQ=255)
                                                                                                     cTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTaa                                                                     >  2:387892/1‑125 (MQ=255)
                                                                                                     cTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTaa                                                                     <  1:387892/125‑1 (MQ=255)
                                                                                                     cTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTc                                                         <  2:1017373/137‑1 (MQ=255)
                                                                                                     cTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTc                                                         >  1:1017373/1‑137 (MQ=255)
                                                                                                     cTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTtct                                            <  1:195253/150‑1 (MQ=255)
                                                                                                      tCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGc                                                                                        <  2:243623/105‑1 (MQ=255)
                                                                                                      tCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGc                                                                                        >  1:243623/1‑105 (MQ=255)
                                                                                                        cAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGa                                         >  1:1001477/1‑150 (MQ=255)
                                                                                                         aaGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTGAGCCTTGCCAGTCTTGTAGTAGTTCTCGaa                                        >  2:483057/1‑150 (MQ=255)
                                                                                                               cccACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGAt                                                                                                   <  1:713399/85‑1 (MQ=255)
                                                                                                               cccACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGAt                                                                                                   >  2:713399/1‑85 (MQ=255)
                                                                                                                ccACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGAAg                                 >  1:525753/1‑150 (MQ=255)
                                                                                                                 cACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCAcc                                                                                                            >  2:1196941/1‑74 (MQ=34)
                                                                                                                 cACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCAcc                                                                                                            <  1:1196941/74‑1 (MQ=34)
                                                                                                                  aCTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAGAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTGTCAGTAAGCCTTGCTAGTCTTGCAGTAGTGCTCGAAGCAGAAGtt                               >  2:514759/1‑150 (MQ=25)
                                                                                                                       gCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGAAGTTGGTcg                          >  2:664816/1‑150 (MQ=255)
                                                                                                                         caccacAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGAAGTTGGTcgcg                        >  1:269333/1‑150 (MQ=255)
                                                                                                                         caccacAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGAAGTTGGTcgcg                        >  2:343363/1‑150 (MQ=255)
                                                                                                                              caACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGAAGTTGGTCGCGTCGAt                   <  2:723759/150‑1 (MQ=255)
                                                                                                                              caACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGAAGTTGGTCGCGTCGAt                   <  2:960395/150‑1 (MQ=255)
                                                                                                                                 ccAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGAAGTTGGTCGCGTCGATGTa                <  1:514759/150‑1 (MQ=255)
                                                                                                                                   aaGTGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTACTCGAAGCAGAAGTTGGTCGCGTCGATGTAGc              <  1:1009375/150‑1 (MQ=255)
                                                                                                                                   aaGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGAAGTTGGTCGCGTCGATGTAGc              <  2:1001477/150‑1 (MQ=255)
                                                                                                                                    aGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGc                                                                                        >  2:778086/1‑75 (MQ=255)
                                                                                                                                    aGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGc                                                                                        <  1:778086/75‑1 (MQ=255)
                                                                                                                                       gggTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCa                                                                         <  1:538244/87‑1 (MQ=255)
                                                                                                                                       gggTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCa                                                                         >  2:538244/1‑87 (MQ=255)
                                                                                                                                          ttGCGCGTAAGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGaa                                  <  2:779225/123‑1 (MQ=255)
                                                                                                                                          ttGCGCGTAAGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGaa                                  >  1:779225/1‑123 (MQ=255)
                                                                                                                                               cTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGAAGTTGGTCGCGTCGATGTAGCCTTCGGCGCCAc  <  1:345112/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGAAGTTGGTCGCGTCGATGTAGCCTTCGGCGCCAC  >  NC_002947/4348934‑4349226

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: