Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,363,476 T→A 5.3% L173Q (CTG→CAG)  adhP → alcohol dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,363,4760TA 5.3% 236.2 / 2.7 75L173Q (CTG→CAG) adhPalcohol dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (37/34);  new base A (2/2);  total (39/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CCGACTACTCGGTTAAAGTGCCCGACGGCCTCGACTCCGCCGCCGCCAGCAGCATCACCTGCGCCGGCGTCACGACCTACAAGGCAGTGAAAATCTCCAACGTCCGCCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCAAAGTCATCGCCATCGACGTCAACGAAGAGCAACTGCGCTTCGCCAGCGAGATGGGCGCGGATCTGGTCGTCAACCCGCTCAACGAAGATGCCGCGAAGGTCATC  >  NC_002947/4363332‑4363624
                                                                                                                                                |                                                                                                                                                    
ccGACTACTCGGTTAAAGTGCCCGACGGCCTCGACTCCGCCGCCGCCAGCAGCATCACCTGCGCCGGCGGCACGACCTACAAGGCAGTGAAAATCTCCAACGTCCGCCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCa                                                                                                                                                 >  1:798776/1‑150 (MQ=255)
  gACTACTCGGTTAAAGTGCCCGACGGCCTCGACTCCGCCGCCGCCAGCAGCATCACCTGCGCCGGCGTCACGACCTACAAGGCAGTGAAAATCTCCAACGTCCGCCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAAc                                                                                                                                               >  2:1906112/1‑150 (MQ=255)
  gACTACTCGGTTAAAGTGCCCGACGGCCTCGACTCCGCCGCCGCCAGCAGCATCACCTGCGCCGGCGTCACGACCTACAAGGCAGTGAAAATCTCCAACGTCCGCCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAAc                                                                                                                                               <  2:476627/150‑1 (MQ=255)
     tactCGGTTAAAGTGCCCGACGGCCTCGACTCCGCCGCCGCCAGCAGCATCACCTGCGCCGGCGTCACGACCTACAAGGCAGTGAAAATCTCCAACGTCCGCCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTg                                                                                                                                            <  2:1381363/150‑1 (MQ=255)
        tCGGTTAAAGTGCCCGACGGCCTCGACTCCGCCGCCGCCAGCAGCATCACCTGCGCCGGCGTCACGACCTACAAGGCAGTGAAAATCTCCAACGTCCGCCCCGGCCAATGGTTCGCCATCTACGGGCTCGGCCGCCTGGGCAACCTGGcc                                                                                                                                         >  1:2004591/1‑150 (MQ=255)
        tCGGTTAAAGTGCCCGACGGCCTCGACTCCGCCGCCGCCAGCAGCAGCACCTGCGCCGGCGTCACGACCTACAAGGCAGTGAAAATCTCCAACGTCCGCCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGcc                                                                                                                                         >  1:1537774/1‑150 (MQ=255)
           gTTAAAGTGCCCGACGGCCTCGACTCCGCCGCCGCCAGCAGCATCACCTGCGCCGGCGTCACGACCTACAAGGCAGTGAAAATCTCCAACGTCCGCCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTg                                                                                                                                      >  1:303798/1‑150 (MQ=255)
            ttAAAGTGCCCGACGGCCTCGACTCCGCCGCCGCCAGCAGCATCACCTGCGCCGGCGTCACGACCTACAAGGCAGTGAAAATCTCCAACGTCCGCCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGc                                                                                                                                     <  1:1604845/150‑1 (MQ=255)
                    cccGACGGCCTCGACTCCGCCGCCGCCAGCAGCATCACCTGCGCCGGCGTCACGACCTACAAGGCAGTGAAAATCTCCAACGTCCGCCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGcc                                                                                                                             <  2:1255054/150‑1 (MQ=255)
                     ccGACGGCCTCGACTCCGCCGCCGCCAGCAGCATCACCTGCGCCGGCGCCACGACCTACAAGGCAGTGAAAATCTCCAACGTCCGCCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCa                                                                                                                            <  2:1565285/150‑1 (MQ=255)
                          ggCCTCGCCTCCGCCGCCGCCAGCAGCACCACCTGCGCCGGCGTCACCACCTACGAGGCAGTAAAAATGTGCGACATGCTGTTCCGCCAACAGATCGCCATCTACGGGCTCGGCGGCCTCGGCAACCTGGCCCTGCAATAGGCCAAGaag                                                                                                                       <  2:1951474/150‑2 (MQ=255)
                           gCCTCGACTCCGCCGCCGCCAGCAGCATCACCTGCGCCGGCGTCACGACCTACAAGGCAGTGAAAATCTCCAACGTCCGCCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAAtg                                                                                                                      >  1:626481/1‑150 (MQ=255)
                              tCGACTCCGCCGCCGCCAGCAGCATCACCTGCGCCGGCGTCACGACCTACAAGGCAGTGAAAATCTCCAACGTCCGCCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATgtgt                                                                                                                   >  2:1178412/1‑150 (MQ=255)
                               cGACTCCGCCGCCGCCAGCAGCATCACCTGCGCCGGCGTCACGACCTACAAGGCAGTGAAAATCTCCAACGTCCGCCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGtt                                                                                                                  <  1:593264/150‑1 (MQ=255)
                                  cTCCGCCGCCGCCAGCAGCATCACCTGCGCCGGCGTCACGACCTACAAGGCAGTGAAAATCTCCAACGTCCGCCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCaa                                                                                                               >  1:934222/1‑150 (MQ=255)
                                  cTCCGCCGCCGCCAGCAGCATCACCTGCGCCGGCGTCACGACCTACAACGCAGTGAAAATCTCCAACGTCCGCCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCaa                                                                                                               <  2:1566646/150‑1 (MQ=255)
                                   tCCGCCGCCGCCAGCAGCATCACCTGCGCCGGCGTCACGACCTACAAGGCAGTGAAAATCTCCAACGTCCGCCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAAc                                                                                                              >  2:317669/1‑150 (MQ=255)
                                       ccgccgccAGCAGCATCACCTGCGCCGGCGTCACGACCTACAAGGCAGTGAAAATCTCCAACGTCCGCCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCa                                                                                                          >  2:1578779/1‑150 (MQ=255)
                                         gcctccAGAAGCATCACCTGCGCCGGCGCCAAGACCTACAAGGCAGTGAAAATCTCCAACGTCCGCCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCaaa                                                                                                        <  2:413228/150‑1 (MQ=255)
                                         gccgccAGCAGCATCACCTGCGCCGGCGTCACGACCTACAAGGCAGTGAAAATCTCCAACGTCCGCCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCaaa                                                                                                        >  1:189653/1‑150 (MQ=255)
                                         gccgccAGCAGCATCACCTGCGCCGGCGTCACGACCTACAAGGCAGTGAAAATCTCCAACGTCCGCCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCaaa                                                                                                        <  2:1701411/150‑1 (MQ=255)
                                             ccAGCAGCATCACCTGCGCCGGCGTCACGACCTACAAGGCAGTGAAAATCTCCAACGTCCGCCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCAAAGTCa                                                                                                    >  1:1751965/1‑150 (MQ=255)
                                              cagcagCATCACCTGCGCCGGCGTCACGACCTACAAGGCAGTGAAAATCTCCAACGTCCGCCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCAAAGTCAt                                                                                                   >  1:826408/1‑150 (MQ=255)
                                               agcagcATCACCTGCGCCGGCGTCACGACCTACAAGGCAGTGAAAATCTCCAACGTCCGCCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCAAAGTCATc                                                                                                  <  2:189653/150‑1 (MQ=255)
                                                        aCCTGCGCCGGCGTCACGACCTACAAGGCAGTGAAAATCTCCAACGTCCGCCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCAAAGTCATCGCCATCGAc                                                                                         >  2:50043/1‑150 (MQ=255)
                                                        aCCTGCGCCGGCGTCACGACCTACAAGGCAGTGAAAATCTCCAACGTCCGCCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCAAAGTCATCGCCATCGAc                                                                                         >  1:1726531/1‑150 (MQ=255)
                                                         ccTGCGCCGGCGTCACGACCTACAAGGCAGTGAAAATCTCCAACGTCCGCCCCGGCCAATGGATCGCCATCTACGGACTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCAAAGTCATCGCCATCGACg                                                                                        >  1:1536630/1‑150 (MQ=255)
                                                          cTGCGCCGGCGTCACGACCTACAAGGCAGTGAAAATCTCCAACGTCCGCCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCAAAGTCATCGCCATCGACGt                                                                                       <  2:826408/150‑1 (MQ=255)
                                                          cTGCGCCGGCGTCACGACCTACAAGGCAGTGAAAATCTCCAACGTCCGCCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCAAAGTCATCGCCATCGACGt                                                                                       <  1:1178412/150‑1 (MQ=255)
                                                               ccGGCGTCACGACCTACAAGGCAGTGAAAATCTCCAACGTCCGCCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCAAAGTCATCGCCATCGACGTCAACg                                                                                  >  1:1327820/1‑150 (MQ=255)
                                                                   cGTCACGACCTACAAGGCAGTGAAAATCTCCAACGTCCGCCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCAAAGTCATCGCCATCGACGTCAACGAaga                                                                              <  2:1537774/150‑1 (MQ=255)
                                                                    gTCACGACCTACAAGGCAGTGAAAATCTCCAACGTCCACCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCAAAGTCATCGCCATCGACGTCAACGAagag                                                                             >  1:1891441/1‑150 (MQ=255)
                                                                       aCGACCTACAAGGCAGTGAAAATCTCCAACGTCCGCCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCAAAGTCATCGCCATCGACg                                                                                        >  1:777645/1‑136 (MQ=255)
                                                                       aCGACCTACAAGGCAGTGAAAATCTCCAACGTCCGCCCCGGCCAATGGATCGCAATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCAAAGTCATCGCCATCGACg                                                                                        <  2:777645/136‑1 (MQ=255)
                                                                         gACCTACAAGGCAGTGAAAATCTCCAACGTCCGCCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCAAAGTCATCGCCATCGACGTCAACGAAGAGCAACt                                                                        <  1:317669/150‑1 (MQ=255)
                                                                             tACAAGGCAGTGAAAATCTCCAACGTCCGCCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCAAAGTCATCGCCATCGACGTCAACGAAGATCAACTgcgc                                                                    <  2:355123/150‑1 (MQ=255)
                                                                                 aGGCAGTGAAAATCTCCAACGTCCGCCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAAc                                                                                                              >  2:997360/1‑104 (MQ=255)
                                                                                 aGGCAGTGAAAATCTCCAACGTCCGCCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAAc                                                                                                              <  1:997360/104‑1 (MQ=255)
                                                                                  ggCAGTGAAAATCTCCAACGTCCGCCCCGGCCAATGGATCGCCATCTACGGGCTCGGAGAACAGGGCAACCTGGCCCTTCAAGATGCCAAGAATGTGTTCAACGCCAAAGTCATCGCCATCGACGTCAACGAAGAGCAACTGCGCTTCGc                                                               >  1:871109/1‑150 (MQ=255)
                                                                                  ggCAGTGAAAATCTCCAACGTCCGCCCCGGCCAATGGATCGCCATCTACGGGCTAGGAAGCCAGGGCAACCTATGCGTGCAATATGCCAAGAATGTGTTCAACGCCAAAGTCATCGCCATCGACGTCAACGAAGAACAACTGCGCTTCGc                                                               >  1:871054/1‑150 (MQ=255)
                                                                                         aaaaTCTCCAACGTCCGCCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCt                                                                                                                                       <  1:1523545/71‑1 (MQ=255)
                                                                                         aaaaTCTCCAACGTCCGCCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCt                                                                                                                                       >  2:1523545/1‑71 (MQ=255)
                                                                                            aTCTCCAACGTCCGCCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCAAAGTCATCGCCATCGACGTCAACGAAGAGCAACTGCGCTTCGCCAGCGAGATg                                                     >  2:1458532/1‑150 (MQ=255)
                                                                                             tctcCAACGTCCGCCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCAAAGTCATCGCCATCGACGTCAACGAAGAGCAACTGCGCTTCGCCAGCGAGATgg                                                    >  1:899185/1‑150 (MQ=255)
                                                                                              ctcCAACGTCCGCCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCAAAGTCATCGCCATCGACGTCAACGAAGAGCAACTGCGCTTCGCCAGCGAGATggg                                                   >  2:1018386/1‑150 (MQ=255)
                                                                                               tcCAACGTCCGCCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCACTATGCCAAGAATGTGTTCAACGCCAAAGTCATCGCCATCGACGTCAACGAAGAGCAACTGCGCTTCGCCAGCGAGATgggg                                                  >  1:333247/1‑149 (MQ=255)
                                                                                                 cAACGTCCGCCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCAAAGTCATCGCCATCGACGTCAACGAAGAGCAACTGCGCTTCGCCAGCGAGATGGgcgc                                                >  2:432525/1‑150 (MQ=255)
                                                                                                 cAACGTCCGCCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCAAAGTCATCGCCATCGACGTCAACGAAGAGCAACTGCGCTTCGCCAGCGAGATGGgcgc                                                <  2:955669/150‑1 (MQ=255)
                                                                                                   aCGTCCGCCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCAAAGTCATCGCCATCGACGTCAACGAAGAGCAACTGCGCTTCGCCAGCgaga                                                       <  2:220135/141‑1 (MQ=255)
                                                                                                   aCGTCCGCCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCAAAGTCATCGCCATCGACGTCAACGAAGAGCAACTGCGCTTCGCCAGCgaga                                                       >  1:220135/1‑141 (MQ=255)
                                                                                                       ccGCCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCAAAGTCATCGCCATCGACGTCAACGTAAAGCAACTGCGCTTCGCCAGCGAGATGGGCGCGGATCt                                          >  1:2239505/1‑150 (MQ=255)
                                                                                                         gCCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCAAAGTCATCGCCATCGACGTCAACGAAGAGCAACTGCGCTTCGCCAGCGAGATGGGCGCGGATCTgg                                        <  1:50043/150‑1 (MQ=255)
                                                                                                           cccGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCAAAGTCATCGCCATCGACGTCAACGAAGAGCAACTGCGCTTCGCCAGCGAGATGGGCGCGGATCTGgtc                                      <  1:1458532/150‑1 (MQ=255)
                                                                                                           cccGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCAGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCAAAGTCATCGCCATCGACGTCAACGAAGAGCAACTGCGCTTCGCCAGCGAGATGGGCGCGGATCTGgtc                                      <  2:871109/150‑1 (MQ=255)
                                                                                                           cccGGCCAATGCATCGCCATCTACGGGCTCGGCGGCCAGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCAAAGTCATCGCCATCGACGTCAACGAAGAGCAACTGCGCTTCGCCAGCGAGATGGGCGCGGATCTGgtc                                      <  2:871054/150‑1 (MQ=255)
                                                                                                             cGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCAAAGTCATCGCCATCGACGTCAACGAAGAGCAACTGCGCTTCGCCAGCGAGATGGGCGCGGATCTGgtcgt                                    <  2:697579/150‑1 (MQ=255)
                                                                                                                   aTGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCAAAGTCATCGCCATCGACGTCAACGAAGAGCAACTGCGCTTCGCCAGCGAGATGGGCGCGGATCTGGTCGTCAAccc                              >  1:71803/1‑150 (MQ=255)
                                                                                                                     ggATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCAAAGTCATCGCCATCGACGTCAACGAAGAGCAACTGCGCTTCGCCAGCGAGATGGGCGCGGATCTGGTCGTCAACCCga                            <  2:626481/150‑2 (MQ=255)
                                                                                                                      gATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCAAAGTCATCGCCATCGACGTCAACGAAGAGCAACTGCGCTTCGCCAGCGAGATGGGCGCGGATCTGGTCGTCAACCCGCt                           >  2:1923188/1‑150 (MQ=255)
                                                                                                                      gATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCAAAGTCATCGCCATCGACGTCAACGAAGAGCAACTGCGCTTCGCCAGCGAGATGGGCGCGGATCTGGTCGTCAACCCGCt                           >  2:1615210/1‑150 (MQ=255)
                                                                                                                        tCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCAAAGTCATCACCATCGACGTCAACGAAGAGCAACTGCGCTTCGCCAGCGAGATGGGCGCGGATCTGGTCGTCACCCCGCTCa                         >  2:1986644/1‑150 (MQ=255)
                                                                                                                           ccATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGaa                                                                                                                        <  1:826362/52‑1 (MQ=255)
                                                                                                                           ccATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGaa                                                                                                                        >  2:826362/1‑52 (MQ=255)
                                                                                                                           ccATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCAAAGTCATCGCCATCGACGTCAACGAAGAGCAACTGCGCTTCGCCa                                                             >  2:1206524/1‑111 (MQ=255)
                                                                                                                           ccATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCAAAGTCATCGCCATCGACGTCAACGAAGAGCAACTGCGCTTCGCCa                                                             <  1:1206524/111‑1 (MQ=255)
                                                                                                                            cATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCAAAGTCATCGCCATCGACGTCAACGAAGAGCAACTGCGCTTCGCCAGCGAGATGGGCGCGGATCTGGTCGTCAACCCGCTCAACGa                     >  1:1701949/1‑150 (MQ=255)
                                                                                                                             aTCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCAAAGTCATCGCCATCGACGTCAACGAAGAGCAACTGCGCTTCGCCAGCGAGATGGGCGCGGATCTGGTCGTCAACCCGCTCAACGaa                    <  1:1615210/150‑1 (MQ=255)
                                                                                                                             aTCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCAAAGTCATCGCCATCGACGTCAACGAAGAGCAACTGCGCTTCGCCAGCGAGATGGGCGCGGATCTGGTCGTCAACCCGCTCAACGaa                    <  1:888568/150‑1 (MQ=255)
                                                                                                                               cTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCAAAGTCATCGCCATCGACGTCAACGAAGAGCAACTGCGCTTCGCCAGCGAGATGGGCGCGGATCTGGTCGTCAACCCGCTCAACGAAGa                  <  2:2004591/150‑1 (MQ=255)
                                                                                                                                  cGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTc                                                                                                                 <  1:1170966/52‑1 (MQ=255)
                                                                                                                                  cGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTc                                                                                                                 >  2:1170966/1‑52 (MQ=255)
                                                                                                                                     gCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCAAAGTCATCGCCATCGACGTCAACGAAGAGCAACTGCGCTTCGCCAGCGAGATGGGCGCGGATCTGGTCGTCAACCCGCTCAACGAAGATGCcgc            >  1:932464/1‑150 (MQ=255)
                                                                                                                                        cggcggCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCAAAGTCATCGCCATCGACGTCAACGAAGAGCAACTGCGCTTCGCCAGCGAGATGGGCGCGGATCTGGTCGTCAACCCGCTCAACGAAGATGCCGCGaa         <  1:1923188/150‑1 (MQ=255)
                                                                                                                                         ggcggcCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCAAAGTCATCGCCATCGACGTCAACGAAGAGCAACTGCGCTTCGCCAGCGAGATGGGCGCGGATCTGGTCGTCAACCCGCTCAACGAAGATGCCGCGAAg        <  2:899185/150‑1 (MQ=255)
                                                                                                                                               cTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCAAAGTCATCGCCATCGACGTCAACGAAGAGCAACTGCGCTTCGCCAGCGAGATGGGCGCGGATCTGGTCGTCAACCCGCTCAACGAAGATGCCGCGAAGGTCATc  <  2:303798/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
CCGACTACTCGGTTAAAGTGCCCGACGGCCTCGACTCCGCCGCCGCCAGCAGCATCACCTGCGCCGGCGTCACGACCTACAAGGCAGTGAAAATCTCCAACGTCCGCCCCGGCCAATGGATCGCCATCTACGGGCTCGGCGGCCTGGGCAACCTGGCCCTGCAATATGCCAAGAATGTGTTCAACGCCAAAGTCATCGCCATCGACGTCAACGAAGAGCAACTGCGCTTCGCCAGCGAGATGGGCGCGGATCTGGTCGTCAACCCGCTCAACGAAGATGCCGCGAAGGTCATC  >  NC_002947/4363332‑4363624

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: ‑