Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,349,078 T→G 6.7% intergenic (‑399/+75) 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,0780TG 6.7% 212.2 / 3.3 73intergenic (‑399/+75)PP_3820/galUgroup II intron‑encoding maturase/UTP‑glucose‑1‑phosphate uridylyltransferase
Reads supporting (aligned to +/- strand):  ref base T (35/33);  new base G (2/3);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGAAGTTGGTCGCGTCGATGTAGCCTTCGGCGCCAC  >  NC_002947/4348929‑4349226
                                                                                                                                                     |                                                                                                                                                    
gTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCt                                                                                                                                                      >  2:359156/1‑150 (MQ=32)
 tCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACTCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCtt                                                                                                                                                     >  1:511360/1‑150 (MQ=17)
     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                                                                                                                         >  2:232024/1‑150 (MQ=32)
                             ttcAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGcaca                                                                                                                         >  1:1202592/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                                                                                                   <  2:641858/150‑1 (MQ=255)
                                                   acacCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGAt                                                                                                   <  1:201969/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                                                                                >  1:254079/1‑142 (MQ=255)
                                                                              ccTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAggg                                                                                <  2:254079/142‑1 (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                                                                                >  2:545807/1‑128 (MQ=255)
                                                                                            tAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAggg                                                                                <  1:545807/128‑1 (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                                                                                        >  1:243623/1‑105 (MQ=255)
                                                                                                           tCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGc                                                                                        <  2:243623/105‑1 (MQ=255)
                                                                                                             cAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGa                                         >  1:1001477/1‑150 (MQ=255)
                                                                                                              aaGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTGAGCCTTGCCAGTCTTGTAGTAGTTCTCGaa                                        >  2:483057/1‑150 (MQ=255)
                                                                                                                    cccACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGAt                                                                                                   >  2:713399/1‑85 (MQ=255)
                                                                                                                    cccACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGAt                                                                                                   <  1:713399/85‑1 (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                        >  2:343363/1‑150 (MQ=255)
                                                                                                                              caccacAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGAAGTTGGTcgcg                        >  1:269333/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)
                                                                                                                                                     |                                                                                                                                                    
GTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGAAGTTGGTCGCGTCGATGTAGCCTTCGGCGCCAC  >  NC_002947/4348929‑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: