Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,344,389 C→G 6.6% 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 6.6% 238.4 / 7.5 76Q261H (CAG→CACPP_2941hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (37/34);  new base G (2/3);  total (39/37)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.49e-01

CCATTCTCGAAGGGAAGGTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGTTGCATCCACTC  >  NC_002947/3344241‑3344538
                                                                                                                                                    |                                                                                                                                                     
ccATTCTCGAAGGGAAGGTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACt                                                                                                                                                      <  2:1036640/150‑1 (MQ=255)
  aTTCTCGAAGGGAAGGTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGa                                                                                                                                                    >  2:1306637/1‑150 (MQ=255)
       cGAAGGGAAGGTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCAGTGAAGCgg                                                                                                                                               <  1:1098076/150‑1 (MQ=255)
        gaagggaaggTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTGGATGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGc                                                                                                                                              <  1:1277609/150‑1 (MQ=255)
              aaggTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTAtc                                                                                                                                        <  2:1206175/150‑1 (MQ=255)
              aaggTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTAtc                                                                                                                                        >  1:1606061/1‑150 (MQ=255)
              aaggTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTAtc                                                                                                                                        >  2:1813888/1‑150 (MQ=255)
                    aTGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCgg                                                                                                                                               <  1:396639/137‑1 (MQ=255)
                    aTGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCgg                                                                                                                                               >  2:396639/1‑137 (MQ=255)
                           aaCAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCtt                                                                                                                           >  2:676814/1‑150 (MQ=255)
                                 aCGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGc                                                                                                                     <  1:1606017/150‑1 (MQ=255)
                                         aGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTATAGTAGCTCAGCGCCTTGATTGCGAACATTg                                                                                                             <  1:51419/150‑1 (MQ=255)
                                          gTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGGCCCCGATGCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTgg                                                                                                            >  2:51419/1‑150 (MQ=255)
                                           tCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGc                                                                                                           >  2:1135341/1‑150 (MQ=255)
                                              gAAAAGTCGAGCTGGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCCCAGCGCCGTGAGTACGAACATTTGCaaa                                                                                                        >  2:802228/1‑149 (MQ=255)
                                              gAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAAt                                                                                                        >  1:415586/1‑150 (MQ=255)
                                                  aGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCt                                                                                                    <  1:1813888/150‑1 (MQ=255)
                                                   gCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCtg                                                                                                   >  2:378285/1‑150 (MQ=255)
                                                   gCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCtg                                                                                                   >  2:906036/1‑150 (MQ=255)
                                                         cTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCa                                                                                             >  1:880678/1‑150 (MQ=255)
                                                         cTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCCCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCa                                                                                             >  2:853900/1‑150 (MQ=255)
                                                          ttGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCAc                                                                                            <  1:1135341/150‑1 (MQ=255)
                                                               ctgcCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGTTTCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGAt                                                                                       <  1:595126/150‑1 (MQ=255)
                                                               ctgcCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGCGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGAt                                                                                       <  1:583544/150‑1 (MQ=255)
                                                                 gcCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATAc                                                                                     <  1:802228/150‑1 (MQ=255)
                                                                 gcCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATAc                                                                                     <  1:378285/150‑1 (MQ=255)
                                                                  ccAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCTACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATAcc                                                                                    >  2:1347110/1‑150 (MQ=255)
                                                                  ccAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATAcc                                                                                    >  1:410610/1‑150 (MQ=255)
                                                                    aGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTTGACCACCTGCGAATCTGTGCATCCGATTCACTGAAGGGCCTCTATAAGCTCAGCTTCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATAGCAc                                                                                  >  1:133007/1‑150 (MQ=255)
                                                                    aGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCAc                                                                                  >  1:743223/1‑150 (MQ=255)
                                                                      gTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTc                                                                                >  1:948642/1‑150 (MQ=255)
                                                                         aGTTCGCCGAGGTCTGATCGGACCTACCGAAGCGACGGAGGAGTCGAAGGATGCCGAAACGATGCGGCCGCTTCACTGAAGCGGGTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGaa                                                                             <  2:415586/150‑1 (MQ=255)
                                                                           ttCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAgg                                                                           >  2:211310/1‑150 (MQ=255)
                                                                                cGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCg                                                                      <  2:1606061/150‑1 (MQ=255)
                                                                                 gAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCgg                                                                     >  1:213398/1‑150 (MQ=255)
                                                                                             gACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACAtt                                                                                                              <  1:627806/97‑1 (MQ=255)
                                                                                             gACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACAtt                                                                                                              >  2:627806/1‑97 (MQ=255)
                                                                                              aCCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCgg                                                                                                                                               <  1:1078997/63‑1 (MQ=255)
                                                                                              aCCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCgg                                                                                                                                               >  2:1078997/1‑63 (MQ=255)
                                                                                               cccACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAgg                                                       <  1:1163797/150‑1 (MQ=255)
                                                                                                ccACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGt                                                      >  1:1615738/1‑150 (MQ=255)
                                                                                                  aCCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTg                                                    >  1:891650/1‑150 (MQ=255)
                                                                                                     gCAACCACGGGCAGGTCTACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAgc                                                 <  1:1347110/150‑1 (MQ=255)
                                                                                                             gggCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTg                                         >  1:235683/1‑150 (MQ=255)
                                                                                                              ggCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGc                                        <  2:54293/150‑1 (MQ=255)
                                                                                                                 aGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTAGAAGACCTCGGCTGACTCGCCCAGGTTGAACGGCATCTTCTCg                                     >  2:1545606/1‑150 (MQ=255)
                                                                                                                  ggTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCgg                                    <  1:211310/150‑1 (MQ=255)
                                                                                                                   gTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAgg                                                                           >  2:915006/1‑110 (MQ=255)
                                                                                                                   gTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAgg                                                                           <  1:915006/110‑1 (MQ=255)
                                                                                                                      gACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCgg                                                                     <  1:1690837/113‑1 (MQ=255)
                                                                                                                      gACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCgg                                                                     >  2:1690837/1‑113 (MQ=255)
                                                                                                                         caccTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGt                               <  1:1741521/148‑1 (MQ=255)
                                                                                                                         caccTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGt                               >  2:1741521/1‑148 (MQ=255)
                                                                                                                           ccTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCt                                  >  2:1870703/1‑143 (MQ=255)
                                                                                                                           ccTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGc                           >  2:1352010/1‑150 (MQ=255)
                                                                                                                           ccTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGc                           <  2:743223/150‑1 (MQ=255)
                                                                                                                           ccTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCGTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCt                                  <  1:1870703/143‑1 (MQ=255)
                                                                                                                              gCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCg                        >  1:1428798/1‑150 (MQ=255)
                                                                                                                              gCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCg                        >  2:13959/1‑150 (MQ=255)
                                                                                                                                      gtgCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTcgcg                <  1:1306637/150‑1 (MQ=255)
                                                                                                                                        agaGCCGCTTCAGTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTc                                      >  1:782777/3‑126 (MQ=255)
                                                                                                                                        agaGCCGCTTCAGTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTc                                      <  2:782777/124‑1 (MQ=255)
                                                                                                                                          aGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTc                                                                                              <  1:231694/68‑1 (MQ=255)
                                                                                                                                          aGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTc                                                                                              >  2:231694/1‑68 (MQ=255)
                                                                                                                                          aGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGtt            <  2:357306/150‑1 (MQ=255)
                                                                                                                                           gCCGCTTCAGTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTcgc                 <  2:657067/144‑1 (MQ=255)
                                                                                                                                           gCCGCTTCAGTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTcgc                 >  1:657067/1‑144 (MQ=255)
                                                                                                                                           gCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAgg                                                       >  2:810950/1‑106 (MQ=255)
                                                                                                                                           gCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAgg                                                       <  1:810950/106‑1 (MQ=255)
                                                                                                                                              gCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGTTGCAt        <  1:676814/150‑1 (MQ=255)
                                                                                                                                               cTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGTTGCATc       <  1:13959/150‑1 (MQ=255)
                                                                                                                                               cTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGCGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGcc                                              >  1:1916329/1‑111 (MQ=255)
                                                                                                                                               cTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGCGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGcc                                              <  2:1916329/111‑1 (MQ=255)
                                                                                                                                                  cACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGTTGCATCCAc    <  2:1428798/150‑1 (MQ=255)
                                                                                                                                                   aCTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCggctg                                                                  <  1:1441496/87‑1 (MQ=255)
                                                                                                                                                   aCTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCggctg                                                                  >  2:1441496/1‑87 (MQ=255)
                                                                                                                                                    cTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTAACAGGTTGCCTCCACTc  >  2:1847305/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCATTCTCGAAGGGAAGGTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGTTGCATCCACTC  >  NC_002947/3344241‑3344538

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