Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,189,322 G→C 8.5% H457Q (CAC→CAG PP_2799 ← class III aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,189,3220GC 8.5% 247.9 / 12.7 82H457Q (CAC→CAGPP_2799class III aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (37/38);  new base C (3/4);  total (40/42)
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.21e-01

GGCCATGATGGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATGCTTGTGAGC  >  NC_002947/3189179‑3189471
                                                                                                                                                |                                                                                                                                                     
ggCCATGATGGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATTACTGCGATGACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCtgt                                                                                                                                                 <  1:654605/150‑1 (MQ=255)
  ccATGATGGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATTAGATGCAGCCTTTGCTCACAGCACCTGCTGTGa                                                                                                                                               <  1:259641/150‑1 (MQ=255)
  ccATGATGGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGa                                                                                                                                               <  1:1088939/150‑1 (MQ=255)
     tgatgGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCACCTGCTGTGAg                                                                                                                                              <  2:266019/148‑1 (MQ=255)
     tgatgGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCACCTGCTGTGAg                                                                                                                                              >  1:266019/1‑148 (MQ=255)
      gatgGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGGCGCGCAGGGCAGTCTTTGCTCACAGCAGGTGCTGTGAgcgc                                                                                                                                           >  1:1330137/1‑150 (MQ=255)
       atgGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAgcgcg                                                                                                                                          >  1:601735/1‑150 (MQ=255)
        tgGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCGCACAGCAGGTGCTGTGAgcgcgc                                                                                                                                         >  2:2199228/1‑150 (MQ=255)
        tgGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGCGAGTAGCCATCAGATGCAGCCTTCGCTCACAGCACCTGCTGTGAgcaaag                                                                                                                                         <  1:1682438/150‑5 (MQ=255)
          gTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGccc                                                                                                                                       >  2:1089785/1‑150 (MQ=255)
               cATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCACCTGCTGTGa                                                                                                                                               <  1:1240081/137‑1 (MQ=255)
               cATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCACCTGCTGTGa                                                                                                                                               >  2:1240081/1‑137 (MQ=255)
                  ggCACCGAGCCCGATCACTGAAATTTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTAATTGATGCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTc                                                                                                                               <  1:1956123/150‑1 (MQ=255)
                    cACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTCCTGTGAGCGCGCCCAATCCGTcgc                                                                                                                             <  1:1530167/150‑1 (MQ=255)
                     aCCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTcgcg                                                                                                                            >  2:1628164/1‑150 (MQ=255)
                      ccGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCAGCAGCTGCACGCTCTTCGCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGt                                                                                                                           <  2:1892406/150‑1 (MQ=255)
                            ccGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTcagcca                                                                                                                     >  2:589153/1‑150 (MQ=255)
                             cGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAgcgcgc                                                                                                                                         >  2:2367576/1‑129 (MQ=255)
                             cGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAgcgcgc                                                                                                                                         <  1:2367576/129‑1 (MQ=255)
                                   cGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCg                                                                                                              >  2:1447972/1‑150 (MQ=255)
                                               aCTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCt                                                                                                  >  1:1394240/1‑150 (MQ=255)
                                                 tCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTc                                                                                                           <  1:1424184/139‑1 (MQ=255)
                                                 tCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTc                                                                                                           >  2:1424184/1‑139 (MQ=255)
                                                 tCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGa                                                                                                >  1:1781060/1‑150 (MQ=255)
                                                  cGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGaa                                                                                               <  2:1985495/150‑1 (MQ=255)
                                                  cGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGaa                                                                                               <  2:1394240/150‑1 (MQ=255)
                                                   gAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAg                                                                                              >  2:1905128/1‑150 (MQ=255)
                                                     aGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGAt                                                                                            <  2:2141477/150‑1 (MQ=255)
                                                      gttttCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAAGCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATc                                                                                           <  2:1568189/148‑1 (MQ=255)
                                                           catcaCTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTTTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTc                                                                                      <  1:1735916/150‑1 (MQ=255)
                                                                cTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGGGCGCGCCCAATCCGTCGCGTTGCGCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGatat                                                                                 >  1:24099/1‑150 (MQ=255)
                                                                    gtTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTc                                                                             <  1:589153/150‑1 (MQ=255)
                                                                        acgcacgGTCCACCCCTCATGATTATCGCCCGCCGTAACTCGAAACCAGATGCAGTCTTTGCTCACAGCAGGTGCTGTCAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTgg                                                                         <  2:1330137/150‑1 (MQ=255)
                                                                        acgcacgGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTgg                                                                         <  2:1381738/150‑1 (MQ=255)
                                                                          gcacgGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGggt                                                                       >  2:1757204/1‑150 (MQ=255)
                                                                          gcacgGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGATCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGggt                                                                       >  1:1782286/1‑150 (MQ=255)
                                                                              ggTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTa                                                                   >  1:791901/1‑150 (MQ=255)
                                                                               gTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTaa                                                                  <  1:1905128/150‑1 (MQ=255)
                                                                                 cAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATc                                                                <  1:1447972/150‑1 (MQ=255)
                                                                                    aCCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAAc                                                             >  1:938231/1‑150 (MQ=255)
                                                                                      ccTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACgg                                                           >  1:1163880/1‑150 (MQ=255)
                                                                                        tCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACggag                                                         <  2:1451289/150‑1 (MQ=255)
                                                                                        tCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACggag                                                         <  2:1847668/150‑1 (MQ=255)
                                                                                          aTGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACggagga                                                       >  1:12064/1‑150 (MQ=255)
                                                                                             aTTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGata                                                                                  >  1:1987905/1‑120 (MQ=255)
                                                                                             aTTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGata                                                                                  <  2:1987905/120‑1 (MQ=255)
                                                                                              ttgttgTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGt                                                                                                                                >  1:2170552/1‑73 (MQ=255)
                                                                                              ttgttgTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGt                                                                                                                                <  2:2170552/73‑1 (MQ=255)
                                                                                                gttgttTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTgg                                                                         >  2:1879988/1‑126 (MQ=255)
                                                                                                gttgttTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTgg                                                                         <  1:1879988/126‑1 (MQ=255)
                                                                                                gttgttTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAg                                                 <  2:1781060/150‑1 (MQ=255)
                                                                                                   gttTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCg                                              >  2:2167769/1‑150 (MQ=255)
                                                                                                    ttttACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGa                                             <  1:1628164/150‑1 (MQ=255)
                                                                                                          gTGAGTAGCCATCAGATGCAGCCTTTGCTCACTGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATc                                           <  1:1933780/146‑1 (MQ=255)
                                                                                                          gTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATc                                           >  2:1933780/1‑146 (MQ=255)
                                                                                                          gTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCgtgt                                       <  2:601735/150‑1 (MQ=255)
                                                                                                            gAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATc                                                                                                                                   <  2:1794751/56‑1 (MQ=255)
                                                                                                            gAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATc                                                                                                                                   >  1:1794751/1‑56 (MQ=255)
                                                                                                                 gCCATCAGATGCAGCCTTTGCTCACAGCAGGTTCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCACCGGAGGACAGAAGGCGATCGTGTCGTAGAt                                <  1:942863/150‑1 (MQ=255)
                                                                                                                 gCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGAt                                >  1:2353040/1‑150 (MQ=255)
                                                                                                                 gCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGAt                                >  1:1065198/1‑150 (MQ=255)
                                                                                                                     aTCAGATGCAGCCTTTGCTCCCAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAAc                                                             <  2:881494/118‑1 (MQ=255)
                                                                                                                     aTCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAAc                                                             >  1:881494/1‑118 (MQ=255)
                                                                                                                          aTGCAGCATTTGCTCACAGCAGGTGCTGTGAGCGAGCACAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGaa                        >  1:802763/1‑150 (MQ=255)
                                                                                                                           tGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAAt                       >  2:1572391/1‑150 (MQ=255)
                                                                                                                             cAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGa                     <  1:1938493/150‑1 (MQ=255)
                                                                                                                             cAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGa                     >  2:930829/1‑150 (MQ=255)
                                                                                                                               gCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATc                   <  1:1225018/150‑1 (MQ=255)
                                                                                                                               gCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATc                   >  2:1225018/1‑150 (MQ=255)
                                                                                                                                ccTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTgg                                                                         >  1:2221390/1‑95 (MQ=255)
                                                                                                                                ccTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTgg                                                                         <  2:2221390/95‑1 (MQ=255)
                                                                                                                                 cTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACgga                                                          >  1:2181883/1‑109 (MQ=255)
                                                                                                                                 cTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACgga                                                          <  2:2181883/109‑1 (MQ=255)
                                                                                                                                  tttGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGc                >  2:1090539/1‑150 (MQ=255)
                                                                                                                                     gCTCACAGCACCTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCa             >  2:641427/1‑150 (MQ=255)
                                                                                                                                       tCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACggag                                                         <  2:938821/104‑1 (MQ=255)
                                                                                                                                       tCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACggag                                                         >  1:938821/1‑104 (MQ=255)
                                                                                                                                         acaGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCCAACGAGAGCAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATGCt         >  1:1351763/1‑150 (MQ=255)
                                                                                                                                           agcagGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATGCTtg       <  2:12064/150‑1 (MQ=255)
                                                                                                                                           agcagGTGCTGTGAGCCCGCTCAATCCGTCGCGTCATCCAGTGTCCGCTCGGAGGCGCTGGAGATCGCCTCGATATCGTCTTGGGTGGTAATCAACGTAGGACAGAAGGCGATCGTGTCGTTGATGGCGCGATTGATCAGCCCGTGCTtg       <  2:1205817/150‑1 (MQ=255)
                                                                                                                                            gcagGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATGCTtgt      <  1:777048/150‑1 (MQ=255)
                                                                                                                                            gcagGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATGCTtgt      <  1:2167769/150‑1 (MQ=255)
                                                                                                                                                gTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATGCTTGTGAGc  <  2:938231/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
GGCCATGATGGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATGCTTGTGAGC  >  NC_002947/3189179‑3189471

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: