Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 3,566,823 +A pseudogene (1174/1275 nt) PA14_RS16215 → fatty acid desaturase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084633,566,8221.A98.8% 2077.7 / 11.7 500pseudogene (1173/1275 nt)PA14_RS16215fatty acid desaturase
Reads supporting (aligned to +/- strand):  ref base . (5/1);  new base A (228/266);  total (233/267)
Fisher's exact test for biased strand distribution p-value = 1.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 536 total aligned reads displayed.

CGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACAAGGCGCGCCCCTCAAAGTTCGGCCCGATGACCAAAACCTGGGGCATGTGGGAAAACTTCCTGATCCTGCGCAACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGATTGCCCAGATCACC  >  NC_008463/3566672‑3566967
                                                                                                                                                       |                                                                                                                                                  
 gggggTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGa                                                                                                                                                   >  1:6294750/1‑151 (MQ=255)
   gggTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGaat                                                                                                                                                 >  2:9080998/1‑149 (MQ=255)
   gggTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGaat                                                                                                                                                 <  2:5493276/151‑3 (MQ=255)
   gggTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGaat                                                                                                                                                 >  2:3261907/1‑149 (MQ=255)
   gggTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGaat                                                                                                                                                 <  2:1847577/151‑3 (MQ=255)
      tGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGaatcgt                                                                                                                                              <  1:10998033/151‑6 (MQ=255)
      tGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGaatcgt                                                                                                                                              >  1:612436/1‑146 (MQ=255)
       gAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAAtcgtc                                                                                                                                             >  2:2079623/1‑151 (MQ=255)
       gAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAAtcgtc                                                                                                                                             <  2:10809352/151‑1 (MQ=255)
        aaCTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAAtcgtcg                                                                                                                                            <  2:242888/151‑1 (MQ=255)
          cTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGAc                                                                                                                                          <  1:4607749/151‑1 (MQ=255)
            tCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTg                                                                                                                                        <  1:3297364/151‑1 (MQ=255)
            tCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTg                                                                                                                                        <  1:6933918/151‑1 (MQ=255)
            tCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTg                                                                                                                                        >  1:2003189/1‑151 (MQ=255)
             cTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTgg                                                                                                                                       <  1:8417370/151‑1 (MQ=255)
             cTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGAGCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGACTCGTCGACTgg                                                                                                                                       >  1:1638656/1‑151 (MQ=255)
               ggCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTa                                                                                                                                     <  2:8638733/151‑1 (MQ=255)
                gcgcATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTAc                                                                                                                                    <  1:7586734/151‑1 (MQ=255)
                  gcATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACcg                                                                                                                                  <  2:1250285/151‑1 (MQ=255)
                  gcATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACcg                                                                                                                                  <  1:9831713/151‑1 (MQ=255)
                  gcATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACcg                                                                                                                                  <  1:3813633/151‑1 (MQ=255)
                   cATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACcgc                                                                                                                                 <  1:3043673/151‑1 (MQ=255)
                    aTCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACcgcg                                                                                                                                <  1:5394183/151‑1 (MQ=255)
                      cgcgCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCgag                                                                                                                              >  1:911702/1‑151 (MQ=255)
                         gcCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAAc                                                                                                                           >  1:10060887/1‑151 (MQ=255)
                           cGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCa                                                                                                                         >  2:7174422/1‑151 (MQ=255)
                            ggAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCAt                                                                                                                        >  1:9514642/1‑151 (MQ=255)
                             gAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATg                                                                                                                       <  1:9579297/151‑1 (MQ=255)
                             gAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATg                                                                                                                       >  1:1182981/1‑151 (MQ=255)
                               gcgcgcGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGcc                                                                                                                     <  2:10997717/151‑1 (MQ=255)
                                   gcgGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGgc                                                                                                                 <  2:9645368/151‑1 (MQ=255)
                                    cgGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGgct                                                                                                                <  1:10944613/151‑1 (MQ=255)
                                     ggCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGgctc                                                                                                               <  1:7814931/151‑1 (MQ=255)
                                        caacgaGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGgctcgct                                                                                                            <  1:8367221/151‑1 (MQ=255)
                                            gaGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCgg                                                                                                        >  1:11212921/2‑151 (MQ=255)
                                            gaGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCgg                                                                                                        <  1:10388045/150‑1 (MQ=255)
                                              ggCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGgag                                                                                                      <  1:40650/151‑1 (MQ=255)
                                               gcgcgcCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGgaga                                                                                                     <  2:6302654/151‑1 (MQ=255)
                                               gcgcgcCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGgaga                                                                                                     <  2:11334895/151‑1 (MQ=255)
                                                   gcCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGa                                                                                                 >  2:7642614/1‑151 (MQ=255)
                                                    ccccTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGAt                                                                                                >  2:8966351/1‑151 (MQ=255)
                                                      ccTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATgc                                                                                              >  2:8170667/1‑151 (MQ=255)
                                                       cTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATgcc                                                                                             >  1:552689/1‑151 (MQ=255)
                                                        tcgaaGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATgccg                                                                                            >  1:5150471/4‑151 (MQ=255)
                                                          gaaGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGACCCATGCCCGGCTCGCTCCGGAGACCGATgccgcc                                                                                          >  1:1836937/2‑151 (MQ=255)
                                                           aaGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCa                                                                                         >  2:3122285/1‑151 (MQ=255)
                                                                cGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAg                                                                                    <  1:8020361/151‑1 (MQ=255)
                                                                cGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAg                                                                                    <  1:3912745/151‑1 (MQ=255)
                                                                 ggCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGc                                                                                   >  1:7630628/1‑151 (MQ=255)
                                                                   cccGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGccc                                                                                 >  2:1689518/1‑151 (MQ=255)
                                                                   cccGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGccc                                                                                 <  1:4833774/151‑1 (MQ=255)
                                                                     cGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCca                                                                               <  2:3957862/151‑1 (MQ=255)
                                                                              agaCCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCCGCCCCACaccaacca                                                                      >  2:5122560/3‑151 (MQ=255)
                                                                              agaCCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACaccatcca                                                                      >  1:3944205/3‑151 (MQ=255)
                                                                              agaCCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACaccaacca                                                                      <  2:1024926/149‑1 (MQ=255)
                                                                               gaCCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACAccaaccaa                                                                     <  1:3134341/150‑1 (MQ=255)
                                                                                aCCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACGGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACcaaccaac                                                                    <  1:10549215/151‑1 (MQ=255)
                                                                                 ccTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCaaccaacc                                                                   <  2:1012017/151‑1 (MQ=255)
                                                                                 ccTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCaaccaacc                                                                   <  2:11516977/151‑1 (MQ=255)
                                                                                   tGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCAGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCg                                                                 >  2:9264331/1‑151 (MQ=255)
                                                                                       gCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTCCCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATc                                                             >  1:4172525/1‑151 (MQ=255)
                                                                                        cATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCACCCCGCATcc                                                            <  2:1464977/151‑1 (MQ=255)
                                                                                          tgtgGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCa                                                          >  1:2243631/1‑151 (MQ=255)
                                                                                          tgtgGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCa                                                          >  2:2968713/1‑151 (MQ=255)
                                                                                               gagaaCTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCCACCAACCCGCATCCCAGGCtc                                                     >  2:7781500/4‑150 (MQ=255)
                                                                                               gagaaCTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTa                                                     >  1:11080449/4‑151 (MQ=255)
                                                                                               gagaaCTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTa                                                     <  1:2723370/148‑1 (MQ=255)
                                                                                                       ccTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGa                                             <  2:5738617/151‑1 (MQ=255)
                                                                                                        cTGATCCTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGAt                                            <  2:275544/151‑1 (MQ=255)
                                                                                                             ccTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATcacac                                       >  2:4092922/1‑151 (MQ=255)
                                                                                                             ccTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATcacac                                       >  1:2536090/1‑151 (MQ=255)
                                                                                                              cTGCGCGACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATcacaca                                      <  1:10018812/151‑1 (MQ=255)
                                                                                                                   cgaccaccTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGaaa                                 <  2:4474164/149‑1 (MQ=255)
                                                                                                                   cgaccaccTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGaaa                                 <  2:8092466/149‑1 (MQ=255)
                                                                                                                    gaccaccTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAAc                                <  1:8811316/150‑1 (MQ=255)
                                                                                                                    gaccaccTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAAc                                <  1:10833529/150‑1 (MQ=255)
                                                                                                                     accaccTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAAcc                               <  2:3752498/151‑1 (MQ=255)
                                                                                                                     accaccTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAAcc                               >  1:114929/1‑151 (MQ=255)
                                                                                                                        accTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCa                            <  2:9221321/151‑1 (MQ=255)
                                                                                                                         ccTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCAt                           <  2:10533523/151‑1 (MQ=255)
                                                                                                                         ccTCTGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCAt                           <  2:6924425/151‑1 (MQ=255)
                                                                                                                          ctctGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATg                          <  1:4257104/151‑1 (MQ=255)
                                                                                                                          ctctGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATg                          >  2:2614636/1‑151 (MQ=255)
                                                                                                                           tctGCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGa                         <  1:2341875/151‑1 (MQ=255)
                                                                                                                              gCCGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATa                      >  2:676129/1‑151 (MQ=255)
                                                                                                                              gCCGCGGCCAGTCGGACCCGTTCGGAATAGTCGACTGGTACCGCGAGAAACATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACAAACCAACACGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATa                      <  1:4417231/151‑1 (MQ=255)
                                                                                                                               ccGCGGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATAc                     <  2:3616305/151‑1 (MQ=255)
                                                                                                                                 gcgGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTc                   <  1:5905917/151‑1 (MQ=255)
                                                                                                                                  cgGCCAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCa                  >  1:3444464/1‑151 (MQ=255)
                                                                                                                                      cAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGAtt              <  2:3496015/151‑1 (MQ=255)
                                                                                                                                      cAGTCGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGAtt              <  1:7142881/151‑1 (MQ=255)
                                                                                                                                          cGGACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACCGGAAACCCCATGAATACTCAGATTGccc          >  2:6439722/1‑151 (MQ=255)
                                                                                                                                           ggACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGATTGCCCa         >  1:9393214/1‑151 (MQ=255)
                                                                                                                                           ggACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGATTGCCCa         >  1:65183/1‑151 (MQ=255)
                                                                                                                                            gACCCGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGATTGCCCAg        >  2:1942813/1‑151 (MQ=255)
                                                                                                                                              cccGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTCCCGCCTGATCACACAGGACACCCCATGAATACTCAGATTGCCCCGAt      >  2:6179440/1‑151 (MQ=255)
                                                                                                                                                cGTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGATTGCCCAGATCa    >  2:4781843/1‑151 (MQ=255)
                                                                                                                                                 gTTCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGATTGCCCAGATCAc   <  1:7423642/151‑1 (MQ=255)
                                                                                                                                                  ttCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGATTGCCCAGATCAcc  >  1:83030/1‑151 (MQ=255)
                                                                                                                                                  ttCGGAATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGATTGCCCAGATCAcc  <  2:5941208/151‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                  
CGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACAAGGCGCGCCCCTCAAAGTTCGGCCCGATGACCAAAACCTGGGGCATGTGGGAAAACTTCCTGATCCTGCGCAACCACCTCTGCCGCGGCCAGTCGGACCCGTTCGGATCGTCGACTGGTACCGCGAGAACCATGCCCGGCTCGCTCCGGAGACCGATGCCGCCAACCAGCCCCACACCAACCAACCCGCATCCCAGGCTACCGCCTGATCACACAGGAAACCCCATGAATACTCAGATTGCCCAGATCACC  >  NC_008463/3566672‑3566967

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