Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,553,315 A→T 11.3% L552Q (CTG→CAG)  PP_4041 ← glycoside hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,553,3150AT11.3% 209.3 / 16.4 71L552Q (CTG→CAG) PP_4041glycoside hydrolase
Reads supporting (aligned to +/- strand):  ref base A (31/32);  new base T (5/3);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 7.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

ATTCGCTACCCACAACGCGCCCGGTCGCGCGCCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCAGGAAGGTGCCCTCTCCCGGCGTCAGGCCGTCGCTGCAG  >  NC_002947/4553168‑4553458
                                                                                                                                                    |                                                                                                                                               
aTTCGCTACCCACAACGCGCCCGGTCGCGCGCCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATAATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAgg                                                                                                                                               >  2:268691/1‑150 (MQ=255)
     cTACCCACAACGCGCCCGGTCGCGCGCCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGaa                                                                                                                                          <  2:85546/150‑2 (MQ=255)
     cTACCCACAACGCGCCCGGTCGCGCGCCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAg                                                                                                                                          <  2:526760/150‑1 (MQ=255)
      tACCCACAACGCGCCCGGTCGCGCGCCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAgg                                                                                                                                         <  2:310803/150‑1 (MQ=255)
        cccACAACGCGCCCGGTCGCGCGCCGGGCACTGTGCGCGGCGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTc                                                                                                                                       <  1:703665/150‑1 (MQ=255)
             acagtgcCCGGCGCGCGACCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGc                                                                                                                                  >  2:99246/6‑150 (MQ=255)
             acagtgcCCGGCGCGCGACCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGGCGttcg                                                                                                                                  >  1:1108948/6‑148 (MQ=255)
              ccgcgCCCGGTCGCCCGCCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGca                                                                                                                                 <  1:199402/149‑1 (MQ=255)
                                 cGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTg                                                                                                               <  1:481893/150‑1 (MQ=255)
                                          ctcgcgGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGACGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAgcgc                                                                                                      <  1:932935/148‑1 (MQ=255)
                                             gcgGTGCAGATTCAGCGCGGCGTTGATGATGCCTGTCTGGCTGAAGGCCTGGGGGAGGTTGCCGAGCATGCGCTGGCCGGCCGGCTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGGGGGTCAGGTGTTGGTACAGCGCtaa                                                                                                   >  2:932935/1‑148 (MQ=255)
                                                  gCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCtt                                                                                              <  2:628553/150‑1 (MQ=255)
                                                    aGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTg                                                                                            <  1:437397/150‑1 (MQ=255)
                                                      aTTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGc                                                                                          >  1:513751/1‑150 (MQ=255)
                                                           gcgcGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTctggcctggc                                                                                     >  2:402883/1‑150 (MQ=255)
                                                             gcgGTGTTGATAATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGgcgc                                                                                   <  1:268691/150‑1 (MQ=255)
                                                             gaggTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGgcgc                                                                                   >  2:481758/3‑150 (MQ=255)
                                                                gtgtTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGcccc                                                                                <  1:24263/150‑1 (MQ=255)
                                                                   ttGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCTGCTGGCCgcaca                                                                                                                                        >  2:237287/1‑87 (MQ=255)
                                                                   ttGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCTGCTGGCCGaaca                                                                                                                                        <  1:237287/91‑4 (MQ=255)
                                                                   ttGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGt                                                                             >  1:628990/1‑150 (MQ=255)
                                                                           gCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGc                                                                                                                                  <  1:423783/89‑1 (MQ=255)
                                                                           gCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGc                                                                                                                                  >  2:423783/1‑89 (MQ=255)
                                                                           gCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATa                                                                     >  2:630274/1‑150 (MQ=255)
                                                                              tATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAAc                                                                  <  2:377468/150‑1 (MQ=255)
                                                                                tCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTAAGGTGTTGGTACAGCGCTTCGGCttt                                                                                             <  2:1137675/121‑1 (MQ=255)
                                                                                tCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTAAGGTGTTGGTACAGCGCTTCGGCttt                                                                                             >  1:1137675/1‑121 (MQ=255)
                                                                                 cTGGCTGAAAGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTc                                                               <  1:121722/150‑1 (MQ=255)
                                                                                  tGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTca                                                              <  1:481758/150‑1 (MQ=255)
                                                                                     cTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGc                                                                                                                                              <  1:588378/67‑1 (MQ=255)
                                                                                     cTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGc                                                                                                                                              >  2:588378/1‑67 (MQ=255)
                                                                                      tGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAgg                                                                                                                                         <  2:1010374/71‑1 (MQ=255)
                                                                                      tGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAgg                                                                                                                                         >  1:1010374/1‑71 (MQ=255)
                                                                                      tGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGcac                                                                                                                                <  1:653129/80‑1 (MQ=255)
                                                                                      tGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGcac                                                                                                                                >  2:653129/1‑80 (MQ=255)
                                                                                      tGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTcagcca                                                          >  1:190517/1‑150 (MQ=255)
                                                                                       gAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTgg                                                                                      <  1:476345/121‑1 (MQ=255)
                                                                                       gAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTgg                                                                                      >  2:476345/1‑121 (MQ=255)
                                                                                          ggCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTCTGGCCTGGCGCCCCAGTAACACATAAACGTCAgccagcca                                                      <  1:402883/150‑1 (MQ=255)
                                                                                             cTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAgaa                                                   >  2:49594/1‑150 (MQ=255)
                                                                                               gggggAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCg                                                                                                                          <  2:598770/77‑1 (MQ=255)
                                                                                               gggggAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCg                                                                                                                          >  1:598770/1‑77 (MQ=255)
                                                                                                ggggAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAgaacga                                                >  2:77251/1‑150 (MQ=255)
                                                                                                 gggAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAgaacgaa                                               <  1:49594/150‑1 (MQ=255)
                                                                                                         gCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCAg                                       >  2:121171/1‑150 (MQ=255)
                                                                                                         gCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCAg                                       <  1:413556/150‑1 (MQ=255)
                                                                                                         gCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCAg                                       >  2:413556/1‑150 (MQ=255)
                                                                                                          ccGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCAgg                                      >  2:159109/1‑150 (MQ=255)
                                                                                                            gAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCAGGaa                                    >  1:520268/1‑150 (MQ=255)
                                                                                                                  gcgcTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAacc                                                                         <  1:173974/107‑2 (MQ=255)
                                                                                                                  gcgcTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAaca                                                                         >  2:173974/1‑107 (MQ=255)
                                                                                                                  gcgcTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCAGGAAGGTGcc                              <  2:190517/150‑1 (MQ=255)
                                                                                                                  gcgcTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCAGGAAGGTGcc                              <  2:625961/150‑1 (MQ=255)
                                                                                                                          cggccggGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCAGGAAGGTGCCCTCTCCCg                      >  1:940605/1‑150 (MQ=255)
                                                                                                                                 gTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCAGGAAGGTGCCCTCTCCCGGCGTCAg               >  2:1009884/1‑150 (MQ=255)
                                                                                                                                 gTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCAGGAAGGTGCCCTCTCCCGGCGTCAg               <  2:520268/150‑1 (MQ=255)
                                                                                                                                     tACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCAGGAAg                                   <  2:276163/126‑1 (MQ=255)
                                                                                                                                     tACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCAGGAAg                                   >  1:276163/1‑126 (MQ=255)
                                                                                                                                      cgacctCGGCCTGCTGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCAGGAAGGTGCCCTCTCCCGGCGTCAGGCCGt          >  1:77988/6‑150 (MQ=255)
                                                                                                                                         gTTCGGCCAGCAGGCCAAGGTAGTAGCACAGGACGGACAGATGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCAGGAAGGTGCCCTCTCCCGGCGTCAGGCCGTCGc       >  2:263932/1‑150 (MQ=255)
                                                                                                                                         cctCGGCCTGCTGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAgaa                                                   >  1:894889/3‑106 (MQ=255)
                                                                                                                                         cctCGGCCTGCTGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAgaa                                                   <  2:894889/104‑1 (MQ=255)
                                                                                                                                          ctCGGCCTGCTGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCAGGAAGGTGCCCTCTCCCGGCGTCAGGc             >  1:159908/2‑143 (MQ=255)
                                                                                                                                          ctCGGCCTGCTGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCAGGAAGGTGCCCTCTCCCGGCGTCAGGc             <  2:159908/142‑1 (MQ=255)
                                                                                                                                          ctCGGCCTGCTGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCAGGAAGGTGCCCTCTCCAGGCGTCAGGCCGTCGCt      >  1:976261/2‑150 (MQ=255)
                                                                                                                                          atCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAGCACACCAGGAAGGTGCCCTCTCCCGGCGTCAGGCCGTCGCt      >  2:826326/2‑150 (MQ=255)
                                                                                                                                              gCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCAGGAAGGTGCCCTCTCCCGGCGTCAGGCCGTCGCTGCAg  <  2:933325/150‑1 (MQ=255)
                                                                                                                                               ccAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCAGGAAGGTGCCCTCTCCCGGCGTCAg               <  1:1203961/136‑1 (MQ=255)
                                                                                                                                               ccAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCAGGAAGGTGCCCTCTCCCGGCGTCAg               <  1:1204309/136‑1 (MQ=255)
                                                                                                                                               ccAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCAGGAAGGTGCCCTCTCCCGGCGTCAg               >  2:1204309/1‑136 (MQ=255)
                                                                                                                                               ccAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCAGGAAGGTGCCCTCTCCCGGCGTCAg               >  2:1203961/1‑136 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
ATTCGCTACCCACAACGCGCCCGGTCGCGCGCCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGGTGTTGGTACAGCGCTTCGGCTTTGGCCTGGCGCCCCAGTAACACATAAACGTCAGCCAGCCAGAACGAACACACCAGGAAGGTGCCCTCTCCCGGCGTCAGGCCGTCGCTGCAG  >  NC_002947/4553168‑4553458

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