Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,344,389 C→G 7.8% Q261H (CAG→CAC)  PP_2941 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,344,3890CG 7.8% 280.5 / 13.7 90Q261H (CAG→CAC) PP_2941hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (45/38);  new base G (3/4);  total (48/42)
Fisher's exact test for biased strand distribution p-value = 7.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

ATTCTCGAAGGGAAGGTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGTTGCATCCAC  >  NC_002947/3344243‑3344536
                                                                                                                                                  |                                                                                                                                                   
aTTCTCGAAGGGAAGGTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGa                                                                                                                                                  >  1:1866247/1‑150 (MQ=255)
 ttCTCGAAGGGAAGGTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGaa                                                                                                                                                 >  2:571005/1‑150 (MQ=255)
     cGAAGGGAAGGTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCAGTGAAGCgg                                                                                                                                             <  1:675820/150‑1 (MQ=255)
      gaagggaaggTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCAGTGAAGCGGc                                                                                                                                            <  2:2086785/150‑1 (MQ=255)
        agggaaggTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGctc                                                                                                                                          >  1:582517/1‑150 (MQ=255)
           gaaggTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGctc                                                                                                                                          <  2:1730035/147‑1 (MQ=255)
           gaaggTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGctc                                                                                                                                          >  1:1730035/1‑147 (MQ=255)
            aaggTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTAta                                                                                                                                      <  2:582517/150‑2 (MQ=255)
                    ggCTGAACAGTACGTGGAAAGCCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTAtcagctcagc                                                                                                                              <  2:1199879/150‑1 (MQ=255)
                     gCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAgcg                                                                                                                             >  2:1580907/1‑150 (MQ=255)
                                  tGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGaa                                                                                                                <  2:1866247/150‑1 (MQ=255)
                                        gTCAGAAAAGCGGAGCTTGTGCTGCCGAGGGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTgg                                                                                                          <  1:233/150‑1 (MQ=255)
                                           aGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCAGTGAAGCGGctgc                                                                                                                                         <  1:340853/116‑3 (MQ=255)
                                           aGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCAGTGAAGCGGctgc                                                                                                                                         >  2:340853/1‑114 (MQ=255)
                                           aGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCAGTGAAGCGGctgc                                                                                                                                         <  1:340392/116‑3 (MQ=255)
                                           aGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCAGTGAAGCGGctgc                                                                                                                                         >  2:340392/1‑114 (MQ=255)
                                           aGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCaa                                                                                                       >  2:471890/1‑150 (MQ=255)
                                            gAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAAt                                                                                                      >  1:665742/1‑150 (MQ=255)
                                                  cGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCtgt                                                                                                <  1:507766/150‑1 (MQ=255)
                                                    gAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGc                                                                                              <  1:1534067/150‑1 (MQ=255)
                                                     aGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCt                                                                                             >  2:1038503/1‑150 (MQ=255)
                                                     aGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCt                                                                                             >  2:1038272/1‑150 (MQ=255)
                                                     aGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAAACGCTTGATTTCGAACCTTCGCAATGGCCGTGcc                                                                                             >  2:295117/1‑149 (MQ=255)
                                                            gctgcCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGa                                                                                      <  1:1038503/150‑1 (MQ=255)
                                                            gctgcCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGa                                                                                      <  1:1038272/150‑1 (MQ=255)
                                                               gcCATGTGAGTCCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCCACCACCTGCGAATCTGTGCAGCCGCTTCACTGATGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATAc                                                                                   <  1:1293957/150‑1 (MQ=255)
                                                                 cAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCa                                                                                 <  1:772246/150‑1 (MQ=255)
                                                                    gTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTc                                                                              >  2:1216401/1‑150 (MQ=255)
                                                                       aGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGaa                                                                           <  2:987030/150‑1 (MQ=255)
                                                                        gTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAg                                                                          >  2:494941/1‑150 (MQ=255)
                                                                         ttCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAgg                                                                         <  1:1216401/150‑1 (MQ=255)
                                                                           cGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGcc                                                                       <  1:494941/150‑1 (MQ=255)
                                                                               gAGGTCTGATCGGACCCACCGCAACCATGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGc                                                                                                                   >  2:1828765/1‑102 (MQ=255)
                                                                               gAGGTCTGATCGGACCCACCGCAACCATGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGc                                                                                                                   <  1:1828765/102‑1 (MQ=255)
                                                                               gAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGc                                                                                                                   <  2:176027/102‑1 (MQ=255)
                                                                               gAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGc                                                                                                                   >  1:176027/1‑102 (MQ=255)
                                                                                 ggTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCggct                                                                 >  1:2291943/1‑150 (MQ=255)
                                                                                   tCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCggctg                                                                >  1:1499517/1‑149 (MQ=255)
                                                                                   tCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCggctg                                                                <  2:1499517/149‑1 (MQ=255)
                                                                                    cTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCggctggc                                                              <  1:950900/150‑1 (MQ=255)
                                                                                     tGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATACCCGCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCggctggct                                                             >  1:1103151/1‑150 (MQ=255)
                                                                                        tCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGc                                                          >  1:164370/1‑150 (MQ=255)
                                                                                         cGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCtgtg                                                                                               <  1:2198527/112‑1 (MQ=255)
                                                                                         cGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCtgtg                                                                                               >  2:2198527/1‑112 (MQ=255)
                                                                                              ccACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTc                                                                     >  2:308743/1‑133 (MQ=255)
                                                                                              ccACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTc                                                                     <  1:308743/133‑1 (MQ=255)
                                                                                              ccACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGt                                                    >  2:1535492/1‑150 (MQ=255)
                                                                                              ccACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGt                                                    >  2:2240819/1‑150 (MQ=255)
                                                                                                aCCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTg                                                  >  1:1166761/1‑150 (MQ=255)
                                                                                                aCCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTAGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGccc                                                        >  2:1289226/1‑144 (MQ=255)
                                                                                                aCCGCAACCACGGGCAGGTCGACCACCTGAGAGCCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTAGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGccc                                                        <  1:1289226/144‑1 (MQ=255)
                                                                                                 ccGCAACCACGGGCGGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCt                                                                      <  2:2012322/129‑1 (MQ=255)
                                                                                                 ccGCAACCACGGGCGGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCt                                                                      >  1:2012322/1‑129 (MQ=255)
                                                                                                 ccGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGc                                                                                                                   >  1:191667/1‑84 (MQ=255)
                                                                                                 ccGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGc                                                                                                                   <  2:191667/84‑1 (MQ=255)
                                                                                                 ccGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACAt                                                                                                             <  2:1395626/90‑1 (MQ=255)
                                                                                                 ccGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACAt                                                                                                             >  1:1395626/1‑90 (MQ=255)
                                                                                                       ccACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCa                                           >  2:538663/1‑150 (MQ=255)
                                                                                                        cACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCAt                                          >  2:1541665/1‑150 (MQ=255)
                                                                                                               aGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCg                                   >  2:595340/1‑150 (MQ=255)
                                                                                                                 gTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCCGATCTATGAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTGGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGc                                 <  1:294634/150‑1 (MQ=255)
                                                                                                                    gACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTg                              <  1:471890/150‑1 (MQ=255)
                                                                                                                        accTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGa                                                                                      <  2:1380471/90‑1 (MQ=255)
                                                                                                                        accTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGa                                                                                      >  1:1380471/1‑90 (MQ=255)
                                                                                                                         ccTGCGAATCTGTGCAGCAGCTTCACTGAAGCGGCTCTATCAGCTAAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTAGCAACCCCAGGTTGACCGCCATCTGCTCGGCTTGTTCGc                         >  2:482195/1‑150 (MQ=255)
                                                                                                                           tGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGccc                       <  2:1103151/150‑1 (MQ=255)
                                                                                                                           tGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGccc                       <  2:164370/150‑1 (MQ=255)
                                                                                                                           tGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGccc                       >  1:1525068/1‑150 (MQ=255)
                                                                                                                            gCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGATGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCg                      <  1:1580907/150‑1 (MQ=255)
                                                                                                                             cGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGcc                                                         >  2:312824/1‑114 (MQ=255)
                                                                                                                             cGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGcc                                                         <  1:312824/114‑1 (MQ=255)
                                                                                                                             cGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCCCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCCCACTCGATACCACTCGAAGGCCTCGGCTGGCCCTCCCAGGCTGAGCGCCAGCTGCTCGGTCTGTTGGcaccg                     >  2:410854/1‑146 (MQ=255)
                                                                                                                              gactCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGGTTGCTCGCCCggg                    >  1:151784/4‑149 (MQ=255)
                                                                                                                              gAACCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGc                    <  1:1186938/150‑1 (MQ=255)
                                                                                                                               aaTCTGTGCAGCCGCTTCACTGAAGCGGCTATATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCg                   >  1:1471552/1‑150 (MQ=255)
                                                                                                                                aTCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGACGTCTGTCGATACCACTCGAAGGAGTCGGCACCGGCGCGTAGCTTTATCGCCATCTGCTCGGCTTGTTCGCCCGGCGt                  >  1:2271591/1‑150 (MQ=255)
                                                                                                                                    gtgCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAg                                                                          >  2:2002100/1‑90 (MQ=255)
                                                                                                                                    gtgCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAg                                                                          <  1:2002100/90‑1 (MQ=255)
                                                                                                                                       gaGCCGCTTCAGTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGt           >  2:190169/2‑150 (MQ=255)
                                                                                                                                       cAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCATCGGCTGGCTCGCCAATGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGt           <  2:1834610/150‑1 (MQ=255)
                                                                                                                                        aGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCggg            >  2:1500531/1‑148 (MQ=255)
                                                                                                                                        aGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCggg            <  1:1500531/148‑1 (MQ=255)
                                                                                                                                          ccGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGTTGc        >  1:1046149/1‑150 (MQ=255)
                                                                                                                                            gCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTTCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGTTGCAt      >  1:1041075/1‑150 (MQ=255)
                                                                                                                                            gCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGTTGCAt      >  2:1002096/1‑150 (MQ=255)
                                                                                                                                             cctCACTGACGCCGCTCTATCAGCTCAGCGCCGTGATTGCGAACCTTGGCATTGGCCGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGTTGCATc     <  2:151784/148‑1 (MQ=255)
                                                                                                                                             cTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGt                             >  2:2209708/1‑126 (MQ=255)
                                                                                                                                             cTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGt                             <  1:2209708/126‑1 (MQ=255)
                                                                                                                                               tCACTGAAGCAACTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGTGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGTTGCATccc   >  1:1144485/1‑149 (MQ=255)
                                                                                                                                                cACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGTTGCATCCAc  <  2:1115411/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
ATTCTCGAAGGGAAGGTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGTTGCATCCAC  >  NC_002947/3344243‑3344536

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