Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 996,354 C→A 7.7% Q415H (CAG→CAT) ‡ PP_0860 ← sulfite reductase flavoprotein subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947996,3540CA 7.7% 286.4 / 9.6 92Q415H (CAG→CAT) ‡PP_0860sulfite reductase flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base C (43/41);  new base A (2/5);  total (45/47)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

CGAGTACCTTGCGCTCGAAGGCGCGAGCGCTGTCGGGCGCTTCGCCGTCGCCGAAGGTACTCACCACGAAGAGCGCGCGATGGGCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCACTCTGGCTGGCAAAGCCGATCAGCCAGCCGTCGCCATTACCGCTGTTGGCGACGTTGCCACGGGCTGCGCGCACTTGACGTTTCTTGCGGCGGCGGTCGAG  >  NC_002947/996206‑996500
                                                                                                                                                      |                                                                                                                                                   
cGAGTACCTTGCGCTCGAAGGCGCGAGCGCTGTCGGGCGCTTCGCCGTCGCCGAAGGTACTCACCACGAAGAGCGCGCGATGGGCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGctgcct                                                                                                                                                    >  2:160804/1‑150 (MQ=255)
    tACCTTGCGCTCGAAGGCGCGAGCGCTGTCGGGCGCTTCGCCGTCGCCGAAGGTACTCACCACGAAGAGCGCGCGATGGGCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCAGCtggcctgccc                                                                                                                                                <  1:51112/150‑3 (MQ=255)
    tACCTTGCGCTCGAAGGCGCGAGCGCTGTCGGGCGCTTCGCCGTCGCCGAAGGTACTCACCACGAAGAGCGCGCGATCGACGCCACGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAg                                                                                                                                                >  1:150677/1‑150 (MQ=255)
       cTTGCGCTCGAAGGCGCGCGCGCTGTCGGCCGCTTCGCCTTCGGCGACGGAACTCACCGCGTCCGATCTCGGCGGTGCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAgctg                                                                                                                                             <  1:2223279/150‑1 (MQ=255)
       cTTGCGCTCGAAGGCGCGAGCGCTGTCGGGCGCTTCGCCGTCGCCGAAGGTACTCACCACGAAGAGCGCGCGATGGGCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAgctg                                                                                                                                             >  1:1292813/1‑150 (MQ=255)
       cTTGCGCTCGAAGGCGCGAGCCCTGTCGCGCGCTAAGACAGCCTTGGAGATGCTCACCACGAAGAGCGCGCCATATACCTTGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAgctg                                                                                                                                             <  1:2340737/150‑1 (MQ=255)
        ttGCGCTCGACGATGCGAGCGCTGTCGGGCGCTTCGCCGTCGCCGAAGGTACTCACCACGAAGAGCGCGCGATGGGCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCAGGCAgctgg                                                                                                                                            <  1:1098690/150‑2 (MQ=255)
        ttGCGCTCGAAGGCGCGAGCGCTGTCGGGCGCTTCGCCGTCGCCGAAGGTACTCACCACGAAGAGCGCGCGATGGGCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAgctgc                                                                                                                                            >  1:1935088/1‑150 (MQ=255)
        ttGCGCTCGAAGGCGCGAGCGCTGTCGGGCGCTTCGCCGTCGCCGAAGGTACTCACCACGAAGAGCGCGCGATGGGCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCAGGCAgctgg                                                                                                                                            <  2:2050275/150‑2 (MQ=255)
         tGCGCTCGAAGGCGCGAGCGCTGTCGGGCGCTTCGCCGTCGCCGAAGGTACTCACCACGAAGAGCGCGCGATGGGCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCAGGCAgctggc                                                                                                                                           <  1:492252/150‑3 (MQ=255)
          gcgcTCGCAGGCGCGAGCGATATCGGACGATTTACCGTCGCCCAACGTACTCACCACGAAGAGCGCGCGATGGGCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTc                                                                                                                                          <  2:264815/150‑1 (MQ=255)
          gcgcTCGAAGGCGCGAGCGCTGTCGGGCGCTTCGCCGTCGCCGAAGGTACTCACCACGAAGAGCGCGCGATGGGCCTGGCGCAGGTCGTTTTCACACAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCAGGCAgctggcc                                                                                                                                          <  1:1430865/150‑4 (MQ=255)
            gcTCGAAGGCGCGAGCGCTGTCGGGCGCTTCGCCGTCGCCGAAGGTACTCACCACGAAGAGCGCGCGATGGGCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCgg                                                                                                                                        >  2:1744653/1‑150 (MQ=255)
              tCGAAGGCGCGAGCGCTGTCGGGCGCTTCGCCGTCGCCGAAGGTACTCACCACGAAGAGCGCGCGATGGGCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGgcg                                                                                                                                      >  1:130212/1‑150 (MQ=255)
              tCGAAGGCGCGAGCGCTGTCGGGCGCTTCGCCGTCGCCGAAGGTACTCACCACGAAGAGCGCGCGATGGGCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGgcg                                                                                                                                      >  2:1114824/1‑150 (MQ=255)
              tCGAAGGCGCGAGCGCTGTCGGGCGCTTCGCCGTCGCCGAAGGTACTCACCACGAAGAGCGCGCGATGGGCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGgcg                                                                                                                                      >  2:1236296/1‑150 (MQ=255)
                 aaGGCGCGAGCGCTGTCGGGCGCTTCGCCGTCGCCGAAGGTACTCACCACGAAGAGCGCGCGATGGGCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGctc                                                                                                                                   >  1:542926/1‑150 (MQ=255)
                    gcgaGAGCGCTGTCGGGCGCTTCGCCGTCGGCGAAGGTACTCACCACGAAGAGCGAGCGATGGGAATGGCGTAGGTCGTTTTCACCCAGGTCGTCAAGGGCGCGCACTTGAACGTGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGc                                                                                                                                <  2:881683/150‑1 (MQ=255)
                       cgAGCGCTGTCGGGCGCTTCGCCGTCGCCGAAGGTACTCACCACGAAGAGCGCGCGATGGGCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAg                                                                                                                             <  1:163427/150‑1 (MQ=255)
                       cgAGCGCTGTCGGGCGCTTCGCCGTCGCCGAAGGTACTCACCACGAAGAGCGCGCGATGGGCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAg                                                                                                                             <  1:1114824/150‑1 (MQ=255)
                       cgAGCGCTATCAGGCGCTTCGCCGTCGCCGAAGGTACTCACCCCGAAGAGCGCGCGATGGGCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGNCTGCAGCTGCTCGGCGCTCTGCCAg                                                                                                                             <  2:1000520/150‑1 (MQ=255)
                            gcTGTCGGGCGCTTCGCCGTCGCCGAAGGTACTCACCACGAAGAGCGCGCGATGGGCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCaa                                                                                                                       >  2:1607872/1‑150 (MQ=255)
                                  gggCGCTTCGCCGTCGCCGAAGGTACTCACCACGAAGAGCGCGCGATGGGCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTc                                                                                                                 <  2:542926/150‑1 (MQ=255)
                                  gggCGCTTCGCCGTCGCCGAAGGTACTCACCACGAAGAGCGCGCGATGGGCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCAGGCAgctggcctg                                                                                                                                        >  1:1647526/1‑123 (MQ=255)
                                  gggCGCTTCGCCGTCGCCGAAGGTACTCACCACGAAGAGCGCGCGATGGGCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCAGGCAgctggcctg                                                                                                                                        <  2:1647526/128‑6 (MQ=255)
                                   ggCGCTTCGCCGTCGCCGAAGGTACTCACCACGAAGAGCGCGCGATGGGCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCg                                                                                                                <  2:2027877/150‑1 (MQ=255)
                                   ggCGCTTCGCCGTAGCAGAAGCTACTTACTACGAAGGGCGCGCCAGGGGTCTGGCCAAGGTCGTTTTCTCCCAGCTCGGAAAGGGCGCGCACCTGAACGGGCAGTCCATCTGCGTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCg                                                                                                                <  2:1159649/150‑1 (MQ=255)
                                    taccTTCGGCGACGGCGAAGGTACTCACCAGGAAGAGCGCGCGATGGGCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCgg                                                                                                               >  1:1045939/4‑150 (MQ=255)
                                      gcTTCGCCGTCGCCGAAGGTACTCACCACGAAGAGCGCGCGATTGGCCGGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCa                                                                                                             >  2:696197/1‑150 (MQ=255)
                                      gcTTCGCCGTCGCCGAAGGTACTCACCACGAAGAGCGCGCGATGGGCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGcc                                                                                                             >  2:1059464/1‑149 (MQ=255)
                                      ccTTCGGCGACGGCGAAGGTACTCACCACGAAGAGCGCGCGATGGGCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGcc                                                                                                             >  2:1219359/2‑149 (MQ=255)
                                       cTTCGCCGTCGCCGAAGGTACTCACCACGAAGAGCGCGCGATGGGCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCaa                                                                                                            >  1:1138634/1‑150 (MQ=255)
                                        ttCGCCGTCGCCGAAGGTACTCACCACGAAGAGCGCGCGATGGGCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCaaa                                                                                                           <  2:130212/150‑1 (MQ=255)
                                         tcgccgtcgccgAAGGTACTCACCACGAAGAGCGCGCGATGGGCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAg                                                                                                          <  2:1702286/150‑1 (MQ=255)
                                             cgtcgccgAAGGTACTCACCACGAAGAGCGCGCGATGGGCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCAc                                                                                                      <  2:2066291/150‑1 (MQ=255)
                                              gtcgccgAAGGTACTCACCACGAAGAGCGCGCGATGGGCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCAct                                                                                                     >  1:617321/1‑150 (MQ=255)
                                                      aGGTACTCACCACGAAGAGCGCGCGATGGGCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCACTctggctgg                                                                                             <  2:617321/150‑1 (MQ=255)
                                                        gTACTCACCACGAAGAGCGCGCGATGGGCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCaaa                                                                                                           >  2:1796268/1‑134 (MQ=255)
                                                        gTACTCACCACGAAGAGCGCGCGATGGGCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCaaa                                                                                                           <  1:1796268/134‑1 (MQ=255)
                                                          aCTCACCACGAAGAGCGCGCGATGGGCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCACTCTGGCTGGCaaa                                                                                         <  2:1393355/150‑1 (MQ=255)
                                                            tCACCACGAAGAGCGCGCGATGGGCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCACTCTGGCTGGCAAAGc                                                                                       >  2:2034449/1‑150 (MQ=255)
                                                            tCACCACGAAGAGCGCGCGATGGGCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCACTCTGGCTGGCAAAGc                                                                                       <  2:1138634/150‑1 (MQ=255)
                                                            ccaccacGAAGAGCGCGCGATGGGCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCACTCTGGCTGGCAAAGc                                                                                       <  2:373462/149‑1 (MQ=255)
                                                               ccacGAAGAGCGCGCGATGGGCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCACTCTGGCTGGCAAAGCCGa                                                                                    >  2:1979393/1‑150 (MQ=255)
                                                                 acGAAGAGCGCGCGATGGGCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCg                                                                                                                >  1:2255271/1‑120 (MQ=255)
                                                                 acGAAGAGCGCGCGATGGGCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCg                                                                                                                <  2:2255271/120‑1 (MQ=255)
                                                                  cGAAGAGCGCGCGATGGGCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCACTCTGGCTGGCAAAGccagagt                                                                                 <  2:2063353/150‑6 (MQ=255)
                                                                   gAAGAGCGCGCGATGGGCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCActc                                                                                                    <  2:1088439/130‑1 (MQ=255)
                                                                   gAAGAGCGCGCGATGGGCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCActc                                                                                                    >  1:1088439/1‑130 (MQ=255)
                                                                   gAAGAGCGCGCGATGGGCCTGGCGCAGGCCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCACTCTGGCTGGCAAAGCCGATcag                                                                                >  1:2108399/1‑150 (MQ=255)
                                                                       gagCGCGCGATGGGCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCACTCTGGCTGGCAAAGCCGATcagcca                                                                             <  1:2034449/150‑1 (MQ=255)
                                                                                  gggCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCACTCTGGCTGGCAAAGCCGATCAGCCAGCCGTCGCCAt                                                                  >  2:1838231/1‑150 (MQ=255)
                                                                                    gCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCACTCTGGCTGGCAAAGCCGATCAGCCAGCCGTCGCCATTa                                                                <  1:1607872/150‑1 (MQ=255)
                                                                                     ccTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCACTCtggctggc                                                                                            >  1:102603/1‑121 (MQ=255)
                                                                                     ccTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCACTCtggctggc                                                                                            <  2:102603/121‑1 (MQ=255)
                                                                                      cTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCAc                                                                                                      <  1:763786/110‑1 (MQ=255)
                                                                                      cTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCAc                                                                                                      >  2:763786/1‑110 (MQ=255)
                                                                                       gggggcAGGTCGTTTTCACCCAGCTCGGCATGGGCGGGGACCTGCACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTTCCAGGCCAACTGTTCGGCAAAGCCACTCTGGCTGGCAAAGCCGATCAGCCAGCCGTCGCCATTACCg                                                             <  2:1739424/146‑1 (MQ=255)
                                                                                           gcAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCACTCTGGCTGGCAAAGCCGATCAGCCAGCCGTCGCCATTACCGCTGt                                                         <  1:2218950/150‑1 (MQ=255)
                                                                                             aGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCACTCTGGCTGGCAAAGCCGATCAGCCAGCCGTCGCCATTACCGCTGTTg                                                       <  2:1045939/150‑1 (MQ=255)
                                                                                               gTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCACTCTGGCTGGCAAAGCCGATCAGCCAGCCGTCGCCATTa                                                                <  1:183162/139‑1 (MQ=255)
                                                                                               gTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCACTCTGGCTGGCAAAGCCGATCAGCCAGCCGTCGCCATTa                                                                >  2:183162/1‑139 (MQ=255)
                                                                                               gTCGTTTTCACACAGCTCGACAATTGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCACTCTGGCTGGCAAAGCCGATCAGCCAGCCGTCGCCATTACCGCTGTTGGc                                                     >  2:1419143/1‑150 (MQ=255)
                                                                                               gTCGTGTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCACTCTGGCTGGCAAAGCCGATCAGCCAGCCGTCGCCATTACCGCTGTTGGc                                                     <  1:174549/150‑1 (MQ=255)
                                                                                                    tttCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCACTCTGGCTGGCAAAGCCGATCAGCCAGCCGTCGCCATTACCGCTGTTGGCGACGt                                                >  1:1062176/1‑150 (MQ=255)
                                                                                                     ttCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCACTCTGGCTGGCAAAGCCGATCAGCCAGCCGTCGCCATTACCGCTGTTGGCGACGtt                                               <  1:1838231/150‑1 (MQ=255)
                                                                                                     ttCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGATGCTAGGCGGATACTGGTGCCTACTGTTCGGCAAAGCCACTCTGCCTGGCAAAGCCGATCAGCCAGCCGTCGCCATTACCGCTGTTGGCGACGtt                                               <  1:1621969/150‑1 (MQ=255)
                                                                                                        aCCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCACTCTGGCTGGCAAAGCCGATCAGCCAGCCGTCGCCATTACCGCTGTTgg                                                      <  2:702079/140‑1 (MQ=255)
                                                                                                        aCCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCACTCTGGCTGGCAAAGCCGATCAGCCAGCCGTCGCCATTACCGCTGTTgg                                                      >  1:702079/1‑140 (MQ=255)
                                                                                                         cccAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTc                                                                                                                 <  2:1774990/80‑1 (MQ=255)
                                                                                                         cccAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTc                                                                                                                 >  1:1774990/1‑80 (MQ=255)
                                                                                                          ccAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCACTCTGGCTGGCAAAGCCGATCAGCCAGCCGTCGCCATTACCGCTGTTGGCGACGTTGCCAc                                          >  1:2328646/1‑150 (MQ=255)
                                                                                                          ccAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCACTCTGGCTGGCAAAGCCGATCAGCCAGCCGTCGCCATTACCGCTGTTGGCGACGTTGCCAc                                          >  2:2023915/1‑150 (MQ=255)
                                                                                                            aGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCACTctggct                                                                                               >  1:439610/1‑95 (MQ=255)
                                                                                                            aGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCACTctggct                                                                                               <  2:439610/95‑1 (MQ=255)
                                                                                                            aGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCACTCTGGCTGGCAAAGCCGa                                                                                    <  1:2151237/106‑1 (MQ=255)
                                                                                                            aGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCACTCTGGCTGGCAAAGCCGa                                                                                    >  2:2151237/1‑106 (MQ=255)
                                                                                                             gCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCACTCTGGCTGGCAAAGCCGATCAGCCAGCCGTCGCCATTACCGCTGTTGGCGACGTTGCCACggg                                       >  2:1549714/1‑150 (MQ=255)
                                                                                                             gCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGCTCGGCAAAGCCACTCTGGCTGGCAAAGCCGATCAGCCAGCCGTCGCCATTACCGCTGTTGGCGACGTTGCCACggg                                       <  1:160804/150‑1 (MQ=255)
                                                                                                                     aGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCACTCTGGCTGGCAAAGCCGATCAGCCAGCCGTCGCCATTACCGCTGTTGGCGACGTTGCCACGGGCTgcgcgc                               >  2:2242588/1‑150 (MQ=255)
                                                                                                                         cgcgcACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCACTCTGGCTGGCAAAGCCGATCAGCCAGCCGTCGCCATTACCGCTGTTGGCGACGTTGCCACGGGCTGCGCGCACtt                           <  1:2194849/150‑1 (MQ=255)
                                                                                                                              aCCTGAACGGGCAGGCCAGCTGCCTGAAGCTGCTCGGCGCTCTGCCAGGCAAACTGTTCGGCAAAGCCACTCTGGCTGGCAAAGCCGATCAGCCAGCCGTCGCCATTACCGCTGTTGGCGACGTTGCCACGGGCTGCGCGCACTTGACGt                      >  1:1377842/1‑150 (MQ=255)
                                                                                                                               ccTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCaaa                                                                                                           >  2:967136/1‑64 (MQ=255)
                                                                                                                               ccTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCaaa                                                                                                           <  1:967136/64‑1 (MQ=255)
                                                                                                                                cTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCACTCTGGCTGGCAAAGCCGATCAGCCAGCCGTCGCCATTACCGCTGTTGGCGACGTTGCCACGGGCTGCGCGCACTTGACGttt                    >  2:2362800/1‑150 (MQ=255)
                                                                                                                                  gAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCACTCTGGCTGGCAAAGCCGATCAGCCAGCCGTCGCCATTACCGCTGTTGGCGACGTTGCCACGGGCTGCGCGCACTTGACGTTTCt                  >  2:1432296/1‑150 (MQ=255)
                                                                                                                                  gAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCACTCTGGCTGGCAAAGCCGATCAGCCAGCCGTCGCCATTACCGCTGTTGGCGACGTTGCCACGGGCTGCGCGCACTTGACGTTTCt                  >  2:467323/1‑150 (MQ=255)
                                                                                                                                   aaCGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCACTCTGGCTGGCAAAGCCGATCAGCCAGCCGTCGCCATTACCGCTGTTGGCGACGTTGCCACGGGCTGCGCGCACTTGACGTTTCtt                 <  2:2108399/150‑1 (MQ=255)
                                                                                                                                   aaCGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCAAACTGTTCGGCCAAGCCACTCTGGCTGGCAAAGCCGATCAGCCAGCCGTCGCCATTACCGCTGTTGGCGACGTTGCCACGGGCTGCGCGCACTTGACGTTTCtt                 >  1:925907/1‑150 (MQ=255)
                                                                                                                                          cgagcAGCTGCAGGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCACTCTGGCTGGCAAAGCCGATCAGCCAGCCGTCGCCATTACCGCTGTTGGCGACGTTGCCACGGGCTGCGCGCACTTGACGTTTCTTgcggcgg          >  1:769543/5‑150 (MQ=255)
                                                                                                                                          aGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCACTCTGGCTGGCAAAGCCGATCAGCCAGCCGTCGCCATTACCGCTGTTGGCGACGTTGCCACGGGCTGCGCGCACTTGACGTTTCTTgcggcgg          <  2:1292813/150‑1 (MQ=255)
                                                                                                                                           ggCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCACTCTGGCTGGCAAAGCCGATCAGCCAGCCGTCGCCATTACCGCTGTTGGCGACGTTGCCACGGGCTGCGCGCACTTGACGTTTCTTgcggcggc         >  1:1646840/1‑150 (MQ=255)
                                                                                                                                            gCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCACTCTGGCTGGCAAAGCCGATCAGCCAGCCGTCGCCATTACCGCTGTTGGCGACGTTGCCACGGGCTGCGCGCACTTGACGTTTCTTgcggcggcg        >  2:1606753/1‑150 (MQ=255)
                                                                                                                                            cccAGCTGCCTGCAGCTGCTCGGCGCTCCGCCAGGCCAACTGTTCGGCAAAGCCACTCTGGCTGGCAAAGCCGATCAGCCAGCCGTCGCCATTACCGCTGTTGGCGACGTTGCCACGGGCTGCGCGCACTTGACGTTTCTTgcggcggcg        <  1:1219359/149‑1 (MQ=255)
                                                                                                                                                 ctgcctgcAGCTGCTCGGCACTCTGCCAGGCCAACTGTTCGGCAAAGCCACTCTGGCTGGCAAAGCCGATCAGCCAGCCGTCGCCATTACCGCTGTTGGCGACGTTGCCACGGGCTGCGCGCACTTGACGTTTCTTGCGGCGGCGGTCGa   >  2:216622/1‑150 (MQ=255)
                                                                                                                                                  tgcctgcAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAATCCACTCTGGCTGTCAAAGCCGATCAGCCAGCCGTCGCCATTACCGCTGTTGGCGACGTTGCCACGGGCTGCGTGCACTTGACGTTTCTTGCGGCGGCGGTCGAg  <  2:498354/150‑1 (MQ=255)
                                                                                                                                                  tgcctgcAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCACTCTGGCTGGCAAAGCCGATCAGCCAGCCGTCGCCATTACCGCTGTTGGCGACGGTGCCACGGGCTGCGCGCACTTGACGTTTCTTGCGGCGGCGGTCGAg  <  1:1043838/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
CGAGTACCTTGCGCTCGAAGGCGCGAGCGCTGTCGGGCGCTTCGCCGTCGCCGAAGGTACTCACCACGAAGAGCGCGCGATGGGCCTGGCGCAGGTCGTTTTCACCCAGCTCGGCAAGGGCGCGCACCTGAACGGGCAGGCCAGCTGCCTGCAGCTGCTCGGCGCTCTGCCAGGCCAACTGTTCGGCAAAGCCACTCTGGCTGGCAAAGCCGATCAGCCAGCCGTCGCCATTACCGCTGTTGGCGACGTTGCCACGGGCTGCGCGCACTTGACGTTTCTTGCGGCGGCGGTCGAG  >  NC_002947/996206‑996500

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