Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,788,127 A→G 7.9% Q136R (CAG→CGG)  ahpF → alkyl hydroperoxide reductase subunit F

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,788,1270AG 7.9% 209.6 / 7.0 76Q136R (CAG→CGG) ahpFalkyl hydroperoxide reductase subunit F
Reads supporting (aligned to +/- strand):  ref base A (37/33);  new base G (3/3);  total (40/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 = 3.03e-01

CATCGTTGGTGCTGGCACTGCTGCAAGTGGGTGGCCACCCGTCAAAGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTACCTCAACGGTGAGGTGTTTGGTCAGGGCCGTATGG  >  NC_002947/2787980‑2788276
                                                                                                                                                    |                                                                                                                                                     
caTCGTTGGTGCTGGCACTGCTGCAAGTGGGTGGCCACCCGTCAAAGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAgg                                                                                                                                                     >  2:1582318/1‑150 (MQ=255)
 aTCGTTGGTGCTGGCACTGCTGCAAGTGGGTGGCCACCCGTCAAAGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACATAATCTTCGCTGTCTTGCCAGAACTGCCCGGACGTCGTCCAGGc                                                                                                                                                    >  1:1114365/1‑150 (MQ=255)
            tGGCACTGCTGCAAGTGGGTGGCCACCCGTCAAAGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGa                                                                                                                                         >  1:739630/1‑150 (MQ=255)
              gCACTGCTGCAAGTGGGTGGCCACCCGTCAAAGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATg                                                                                                                                       >  2:1617005/1‑150 (MQ=255)
                  tgctgcAAGTGGGTGGCCACCCGTCAAAGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCgg                                                                                                                                   >  1:65563/1‑150 (MQ=255)
                    ctgcAAGTGGGTGGCCACCCGTCAAAGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTa                                                                                                                                 >  2:102311/1‑150 (MQ=255)
                               tggCCACCCGTCAAAGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCaa                                                                                                                      >  2:1877221/1‑150 (MQ=255)
                                gtccgcccGTCAAAGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAAc                                                                                                                     <  2:731737/145‑1 (MQ=255)
                                ggCCACCCGTCAAAGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAAc                                                                                                                     <  2:603958/150‑1 (MQ=255)
                                  ccACCCGTCAAAGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGt                                                                                                                   >  1:1303799/1‑150 (MQ=255)
                                    aCCCGTCAAAGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTgc                                                                                                                 <  2:1074338/150‑1 (MQ=255)
                                    aCCCGTCAAAGGCCAGTGCCGAAGTATTGTGGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTgc                                                                                                                 >  1:1243460/1‑150 (MQ=255)
                                     cccGTCAAAGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTgcg                                                                                                                >  1:264585/1‑150 (MQ=255)
                                      ccGTCAAAGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGACGGTACTGATGGCTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTgcgc                                                                                                               >  2:2433153/1‑150 (MQ=255)
                                       cGTCAAAGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGcc                                                                                                              >  2:1647464/1‑150 (MQ=255)
                                        gTCAAAGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCa                                                                                                             >  1:1629456/1‑150 (MQ=255)
                                         tCAAAGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCAt                                                                                                            <  1:473169/150‑1 (MQ=255)
                                            aaGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGtt                                                                                                         >  2:1370472/1‑150 (MQ=255)
                                            aaGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGtt                                                                                                         >  2:757798/1‑150 (MQ=255)
                                             aGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTg                                                                                                        <  1:106691/150‑1 (MQ=255)
                                             aGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTg                                                                                                        >  2:1588307/1‑150 (MQ=255)
                                               gCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCACGCTGTCGTGCCAGAACTGGCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGGACGCAAACCCAACGTGCGTCACGTTGcc                                                                                                      >  2:731779/1‑150 (MQ=255)
                                                        gAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGt                                                                                             <  1:1037371/150‑1 (MQ=255)
                                                             gATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCt                                                                                        >  1:56607/1‑150 (MQ=255)
                                                              aTCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTAc                                                                                                                                <  2:1597980/109‑1 (MQ=255)
                                                              aTCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTAc                                                                                                                                >  1:1597980/1‑109 (MQ=255)
                                                               tCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTGAAACCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGt                                                                                      >  1:1307982/1‑150 (MQ=255)
                                                                 gAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGtta                                                                                    >  1:1273691/1‑149 (MQ=255)
                                                                 gAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTc                                                                                    >  1:2351943/1‑150 (MQ=255)
                                                                    cAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTTCTCAACCCGAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAg                                                                                 <  1:1617005/150‑1 (MQ=255)
                                                                      gATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGa                                                                               <  2:2295614/150‑1 (MQ=255)
                                                                              cccTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGa                                                                               >  2:154968/1‑142 (MQ=255)
                                                                              cccTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGa                                                                               <  1:154968/142‑1 (MQ=255)
                                                                               ccTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCTAGGGCGAAGTGGa                                                                      >  1:2162091/1‑150 (MQ=255)
                                                                                  ggAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCACCGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATc                                                                   <  2:1054894/150‑1 (MQ=255)
                                                                                  ggAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATc                                                                   <  1:1877221/150‑1 (MQ=255)
                                                                                  ggAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATc                                                                   <  1:835917/150‑1 (MQ=255)
                                                                                       gTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCa                                                              >  2:2430092/1‑150 (MQ=255)
                                                                                          aaTTCAACTTCGAAACCTACTTCTCGCTGTCGTGCAAGTGCTGCCCGGACGTCGTTCAGCCAGTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGa                                                           <  2:1145734/150‑1 (MQ=255)
                                                                                          aaTTCAACTTCGAAACCTACTTCTCGCTGTCATGCCAGAACTGCCCGCTATCGAGGTAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGa                                                           <  2:1243460/150‑1 (MQ=255)
                                                                                              cAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATcct                                                       <  1:1402788/150‑3 (MQ=255)
                                                                                               aaCTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATg                                                      >  1:590788/1‑150 (MQ=255)
                                                                                               aaCTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATg                                                      >  1:2384683/1‑150 (MQ=255)
                                                                                                   tCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCgg                                                  <  2:739630/150‑1 (MQ=255)
                                                                                                          cTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAg                                           <  1:1647464/150‑1 (MQ=255)
                                                                                                          cTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAg                                           <  2:1114365/150‑1 (MQ=255)
                                                                                                            aCTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTa                                         <  1:757798/150‑1 (MQ=255)
                                                                                                              ttCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCaa                                                                                                                      <  2:2116817/71‑1 (MQ=255)
                                                                                                              ttCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCaa                                                                                                                      >  1:2116817/1‑71 (MQ=255)
                                                                                                                 tcGCTGTCATGCCCGAACTGCACGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTAc                                    <  2:1307982/150‑1 (MQ=255)
                                                                                                                   gCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTACCt                                  <  2:373326/150‑1 (MQ=255)
                                                                                                                   gCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTACCt                                  >  2:1784289/1‑150 (MQ=255)
                                                                                                                            ccAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTTCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTACCTCAACGGTGa                         >  2:1355427/1‑150 (MQ=255)
                                                                                                                            ccAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTACCTCAACGGTGa                         <  1:1588307/150‑1 (MQ=255)
                                                                                                                            ccAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTACCTCAACGGTGa                         <  1:102311/150‑1 (MQ=255)
                                                                                                                            ccAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTACCTCAACGGTGa                         >  1:350642/1‑150 (MQ=255)
                                                                                                                             cAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCATGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTACCTCAACGGTGAg                        <  1:426163/150‑1 (MQ=255)
                                                                                                                             cAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTACCTCAACGGTGAg                        <  1:1582318/150‑1 (MQ=255)
                                                                                                                                  gtGCCTGGACGACGTCCGGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTACCTCAACGGTGAGGTGtt                   >  2:793143/2‑150 (MQ=255)
                                                                                                                                  gtGCCTGGACGACGTCCGGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTACCTCAACGGTGAGGTGtt                   >  2:1824737/2‑150 (MQ=255)
                                                                                                                                   tGCCTGGACGACGTCCGGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTACCTCAACGGTGAGGTGttt                  <  1:1376254/150‑1 (MQ=255)
                                                                                                                                   tGCCTGGACGACGTCCGGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTACCTCAACGGTGAGGTGttt                  <  1:793143/150‑1 (MQ=255)
                                                                                                                                   tGCCTGGACGACGTCCGGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTACCTCAACGGTGAGGTGttt                  <  1:184209/150‑1 (MQ=255)
                                                                                                                                   tGCCCGGACGTCGTCCAGGCACTGAACCTTATGCCGGTACTCCACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTACCTCAACGGTGAGGTGttt                  <  1:731779/150‑1 (MQ=255)
                                                                                                                                   tGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGGAAGATGATGGCGGTACCGAGTATCTACCTCAACGGTGAGGTGttt                  <  2:264585/150‑1 (MQ=255)
                                                                                                                                   tGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTACCTCAACGGTGAGGTGttt                  <  2:2249950/150‑1 (MQ=255)
                                                                                                                                    gCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTACCTCAACGGTGAGGTGTTTg                 <  2:56607/150‑1 (MQ=255)
                                                                                                                                    gCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTACCTCAACGGTGAGGTGTTTg                 <  1:1784289/150‑1 (MQ=255)
                                                                                                                                     cctggACGACGTCCGGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTACCTCAACGGTGAGGTGTTTgg                >  1:457522/4‑150 (MQ=255)
                                                                                                                                     cccGGACGTCGTCCAGGCACTGACCCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTACCTCAACGGTGAGGTGTTTgg                >  1:918729/1‑150 (MQ=255)
                                                                                                                                     cccGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTACCTCAACGGTGAGGTGttt                  <  1:2072814/148‑1 (MQ=255)
                                                                                                                                     cccGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTACCTCAACGGTGAGGTGttt                  >  2:2072814/1‑148 (MQ=255)
                                                                                                                                      ccGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTACCTCAACGGTGAGGTGTTTGGt               >  2:2135060/1‑150 (MQ=255)
                                                                                                                                           cgtcgtCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTACCTCAACGGTGAGGTGTTTGGTCAggg          >  1:1098775/1‑150 (MQ=255)
                                                                                                                                               gtcCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTACCTCAACGGTGAGGTGTTTGGTCAGGTCCGt      >  1:335462/1‑150 (MQ=255)
                                                                                                                                                   cAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTACCTCAACGGTGAGGTGTTTGGTCAGGGCCGTATg   >  1:2330026/1‑150 (MQ=255)
                                                                                                                                                    aGGCACTGAACCTGATGGCGGTACTCAGCCCCAACGTGCGCCATGTTGGCATTGACTGTTCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGTGTATCTACCTCAACGGTGAGGTGTTTGGTCAGGGCCGTATgg  <  2:65563/150‑1 (MQ=255)
                                                                                                                                                    aGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAg                                                                          <  1:1330838/78‑1 (MQ=255)
                                                                                                                                                    aGGCACTGAACCTGATGGCGGTACGCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAg                                                                          >  2:1330838/1‑78 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CATCGTTGGTGCTGGCACTGCTGCAAGTGGGTGGCCACCCGTCAAAGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTACCTCAACGGTGAGGTGTTTGGTCAGGGCCGTATGG  >  NC_002947/2787980‑2788276

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: