Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 36,172 A→T 6.0% intergenic (+57/+26) PP_0032 → / ← PP_0033 hypothetical protein/undecaprenyl‑glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294736,1720AT 6.0% 230.9 / 2.4 77intergenic (+57/+26)PP_0032/PP_0033hypothetical protein/undecaprenyl‑glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (37/35);  new base T (2/3);  total (39/38)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

CGGACCAGAACAAACAACAAGTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGCACAGCAGGAGCTGGCTGGCGTTCGCTTCGGTGATCGATG  >  NC_002947/36039‑36319
                                                                                                                                      |                                                                                                                                                       
cGGACCAGAACAAACAACAAGTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCtttttt                                                                                                                                           <  1:200985/148‑1 (MQ=255)
cGGACCAGAACAAACAACAAGTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCtttttt                                                                                                                                           <  7:142790/148‑1 (MQ=255)
  gACCAGAACAAACAACAAGTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGTTTTACCCTGGCTTTTTtgtt                                                                                                                                        <  1:11319/149‑2 (MQ=255)
   aCCAGAACAAACAACAAGTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGc                                                                                                                                       >  5:318406/1‑149 (MQ=255)
     cAGAACAAACAACAAGTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCtt                                                                                                                                     <  5:297703/149‑1 (MQ=255)
     cAGAACAAACAACAAGTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCtt                                                                                                                                     >  6:297703/1‑149 (MQ=255)
      aGAACAAACAACAAGTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGACAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCtt                                                                                                                                     <  1:62433/149‑1 (MQ=255)
       gAACAAACAACAAGTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCt                                                                                                                                      >  4:36150/1‑146 (MQ=255)
       gAACAAACAACAAGTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCt                                                                                                                                      <  3:36150/146‑1 (MQ=255)
       gAACAAACAACAAGTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTg                                                                                                                                   >  2:168201/1‑148 (MQ=255)
        aacaaacaACAAGTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTgg                                                                                                                                  >  6:34534/1‑148 (MQ=255)
           aaacaaCAAGTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAAtt                                                                                                                             >  5:152692/1‑149 (MQ=255)
           aaacaaCAAGTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAAtt                                                                                                                             >  3:41800/1‑149 (MQ=255)
           aaacaaCAAGTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAAtt                                                                                                                             >  7:250446/1‑149 (MQ=255)
              caacaaGTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCtttttt                                                                                                                                           >  8:120905/1‑134 (MQ=255)
              caacaaGTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCtttttt                                                                                                                                           <  7:120905/134‑1 (MQ=255)
                 caaGTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGATTTTTTTGGCTGGGAATTTTGCCt                                                                                                                       >  7:146228/1‑149 (MQ=255)
                   aGTCGCTCAAGAGGAAATGAATGTAAGCAAGGATGGTATGTATCAAGCCGACAATTAATCGGCGGTTAATTAGGACTTCAAATCAGGAATGAAATCGTCAAAAAAGCGATCGGATAAGCTGGCTTTTTTGTGGTTGTAATTTTGCCTTc                                                                                                                     <  8:267084/149‑1 (MQ=255)
                         tCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTg                                                                                                                                   >  3:102339/1‑130 (MQ=255)
                         tCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTg                                                                                                                                   <  4:102339/130‑1 (MQ=255)
                          cAAGATCAAAAGAAAGCAACCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTg                                                                                                              <  8:146228/149‑1 (MQ=255)
                                cAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGc                                                                                                        >  8:250328/1‑149 (MQ=255)
                                  aaagaaagaaagCAAGGATGGTAAGGAACAAGCCGACAATTAAGCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGAGAGGCGt                                                                                                      >  3:17619/1‑149 (MQ=255)
                                   aagaaagaaagCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTtgt                                                                                                                                         >  7:256275/1‑115 (MQ=255)
                                   aagaaagaaagCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTtgt                                                                                                                                         <  8:256275/115‑1 (MQ=255)
                                   aagaaagaaagCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGtt                                                                                                     >  2:44729/1‑149 (MQ=255)
                                      aaagaaagCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGAAGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGtttac                                                                                                 >  1:40463/1‑147 (MQ=255)
                                         gaaagCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCa                                                                                             <  4:146654/149‑1 (MQ=255)
                                          aaagCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCAt                                                                                            <  8:250446/149‑1 (MQ=255)
                                          aaagCAACGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCAt                                                                                            <  8:262891/149‑1 (MQ=255)
                                           aagCAAGGATTGTAATGAACAAGCAGAAAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCAtt                                                                                           <  6:273091/149‑1 (MQ=255)
                                           aagCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCAt                                                                                            <  6:318406/148‑1 (MQ=255)
                                              cAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTg                                                                                          >  7:294388/1‑147 (MQ=255)
                                              cAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTg                                                                                          <  8:294388/147‑1 (MQ=255)
                                              cAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCt                                                                                        <  1:168201/149‑1 (MQ=255)
                                                     ggTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCg                                                                                 >  8:301130/1‑149 (MQ=255)
                                                            aaCAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGAc                                                                          >  8:187735/1‑149 (MQ=255)
                                                                 gCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCg                                                                     <  7:187735/149‑1 (MQ=255)
                                                                 gCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGAGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCg                                                                     >  8:186698/1‑149 (MQ=255)
                                                                   cGACAATTAATCCGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGCGTAAAACCTGGCTTTTATGTCCATGGAATTTAGCCTTCGTACTTGAGAGGCGTTTCCGCTTAGCTTTCTCCGGCCAGACCCACgcc                                                                   >  6:94655/1‑147 (MQ=255)
                                                                         ttAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTtg                                                             >  4:246538/1‑149 (MQ=255)
                                                                             tCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCg                                                         >  1:35464/1‑149 (MQ=255)
                                                                               ggctgTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGTGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCggg                                                       <  6:140885/149‑1 (MQ=255)
                                                                                  ggTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCg                                                    >  6:5610/1‑149 (MQ=255)
                                                                                      aaaTAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGttt                                                 <  1:44729/148‑1 (MQ=255)
                                                                                        aTAGGACTTCAAATGACGTAAGAAATCGACAAAAAAGCCATGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCAGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTcc                                               <  1:237428/149‑1 (MQ=255)
                                                                                         tAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCaa                                             <  4:41800/149‑1 (MQ=255)
                                                                                            gACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATgg                                          <  3:246538/149‑1 (MQ=255)
                                                                                              cTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGca                                        <  3:230389/149‑1 (MQ=255)
                                                                                               ttCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGc                                                     >  8:157703/1‑135 (MQ=255)
                                                                                               ttCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGc                                                     <  7:157703/135‑1 (MQ=255)
                                                                                                  aaaTGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTCTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCgggg                                                      <  1:67599/132‑1 (MQ=255)
                                                                                                  aaaTGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCgggg                                                      >  2:67599/1‑131 (MQ=255)
                                                                                                  aaaTGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCAAATGGCAcag                                     >  6:114872/1‑148 (MQ=255)
                                                                                                       aCGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTAGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGCAcagcag                                  >  1:79898/1‑147 (MQ=255)
                                                                                                         gaaagaaaTCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATgg                                          <  2:311618/136‑1 (MQ=255)
                                                                                                         gaaagaaaTCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATgg                                          >  1:311618/1‑136 (MQ=255)
                                                                                                         gaaagaaaTCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGCACAGCAGGAg                               >  5:54218/1‑147 (MQ=255)
                                                                                                          aaagaaaTCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGCAcagca                                   >  1:164828/1‑142 (MQ=255)
                                                                                                          aaagaaaTCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGCAcagca                                   <  2:164828/142‑1 (MQ=255)
                                                                                                            agaaaTCGACAAAAAAGCCAGCGTAAAACCTCCCTTTATTGTGCTTGGAATAGTAGCCTTCGTACTTCCGAGGCCTTTCCGCAATGCTTACTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGCACAGCAGGAgctcg                           >  1:84756/1‑147 (MQ=255)
                                                                                                              aaaTCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGcaca                                      <  7:21438/135‑1 (MQ=255)
                                                                                                              aaaTCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGcaca                                      >  8:21438/1‑135 (MQ=255)
                                                                                                                  agcaCAAAAAAGCCAGGTTTTACCCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGCACAGCAGGAGCTGGCTGGCg                     >  3:120008/4‑149 (MQ=39)
                                                                                                                       aaaaacGCCAGGGTAAAACCTGGCTTTTTTGTGCTTTGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCACATTGCTTTCTCCGGCCAGACCCTCGGCACCTtgtg                                                           <  6:212729/106‑1 (MQ=255)
                                                                                                                       aaaaaaGCCAGGTTTTACCCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATgg                                          >  1:70194/1‑123 (MQ=255)
                                                                                                                       aaaaaaGCCAGGTTTTACCCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATgg                                          <  2:70194/123‑1 (MQ=255)
                                                                                                                       aaaaaaGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGCACAGCAGGAGCTGGCTGGCGTTCGc                >  8:156986/1‑149 (MQ=255)
                                                                                                                          aaaGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGCACAGCAGGAgctggctg                        >  2:275634/1‑138 (MQ=255)
                                                                                                                            aGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGCACAGCAGGAGCTGGCTGGCGTTCGCTTCgg           <  7:250328/149‑1 (MQ=255)
                                                                                                                            aGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGCACAGCAGGAGCTGGCTGGCGTTCGCTTCgg           <  6:152692/149‑1 (MQ=255)
                                                                                                                                aGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGCACAGCAGGAGCTGGCTGGCGTTCGCTTCGGTGAt       <  2:79898/149‑1 (MQ=255)
                                                                                                                                  ggTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCTTTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATg                                           <  1:241779/111‑1 (MQ=255)
                                                                                                                                  ggTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATg                                           >  2:241779/1‑111 (MQ=255)
                                                                                                                                    tAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGCACAGCAGGAGCTGGCTGGCGTTCGCTTCGGTGATCGAt   <  3:96495/149‑1 (MQ=255)
                                                                                                                                    tAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGCACAGCAGGAGCTGGCTGGCGTTCGCTTCGGTGATCGAt   >  5:217955/1‑149 (MQ=255)
                                                                                                                                    tAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGCACAGCAGGAGCTGGCTGGCGTTCGCTTCGGTGATCGAt   <  6:54218/149‑1 (MQ=255)
                                                                                                                                     aaagCCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGCACAGCAGGAGCTGGCTGGCGTTCGCTTCGGTGATCGAtg  >  2:305388/1‑149 (MQ=255)
                                                                                                                                      |                                                                                                                                                       
CGGACCAGAACAAACAACAAGTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGCACAGCAGGAGCTGGCTGGCGTTCGCTTCGGTGATCGATG  >  NC_002947/36039‑36319

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: