Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,533,536 A→G 100% A83A (GCT→GCC yjhQ ← GNAT family putative N‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,533,5360AG96.8% 502.5 / 2.7 153A83A (GCT→GCCyjhQGNAT family putative N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (2/3);  new base G (85/63);  total (87/66)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.82e-01

Only 100 of 157 total aligned reads displayed.

CCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGCAGC  >  NC_000913/4533439‑4533635
                                                                                                 |                                                                                                   
  gTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGa                                                                                                 >  2:501257/1‑100 (MQ=255)
       gCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGAtg                                                                                            <  2:1860358/100‑1 (MQ=255)
       gCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGAtg                                                                                            <  2:1404857/100‑1 (MQ=255)
        cAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGAtg                                                                                            <  1:823748/99‑1 (MQ=255)
          gCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGc                                                                                         >  1:721531/1‑100 (MQ=255)
           cccATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCa                                                                                        <  2:1309364/100‑1 (MQ=255)
           cccATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCa                                                                                        <  1:742975/100‑1 (MQ=255)
           cccATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCa                                                                                        <  1:3763862/100‑1 (MQ=255)
           cccATCAACCTCAGGTGCTCCATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCa                                                                                        <  2:1498289/100‑1 (MQ=255)
            ccATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGCCCGCCAGGGGGGCGAGGATGTGCAt                                                                                       >  1:1186719/1‑100 (MQ=255)
             cATCAACCTCAGGTGTTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATc                                                                                      >  2:1021190/1‑100 (MQ=255)
             cATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATc                                                                                      <  1:1739765/100‑1 (MQ=255)
             cATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATc                                                                                      >  2:2087261/1‑100 (MQ=255)
             cATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATc                                                                                      <  1:1263924/100‑1 (MQ=255)
               tCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCaa                                                                                    <  2:3402770/100‑1 (MQ=255)
                 aaCCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACg                                                                                  <  1:1835821/100‑1 (MQ=255)
                 aaCCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACg                                                                                  >  2:2547295/1‑100 (MQ=255)
                  aCCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACgg                                                                                 >  2:2729713/1‑100 (MQ=255)
                   ccTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGCTGTGCATCAACGGt                                                                                >  2:2660577/1‑100 (MQ=255)
                   ccTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGt                                                                                >  2:298280/1‑100 (MQ=255)
                    cTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTg                                                                               >  2:2046616/1‑100 (MQ=255)
                    cTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTg                                                                               <  2:1169946/100‑1 (MQ=255)
                     tCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTg                                                                               >  1:704317/1‑99 (MQ=255)
                     tCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGa                                                                              >  2:872354/1‑100 (MQ=255)
                      cAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGaa                                                                             >  2:1093088/1‑100 (MQ=255)
                       aGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAAt                                                                            <  1:3487335/100‑1 (MQ=255)
                       aGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAAt                                                                            >  1:3125366/1‑100 (MQ=255)
                       aGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAAt                                                                            <  1:947318/100‑1 (MQ=255)
                       aGGTGCTCGATGCCAGTGCGTAGCAACCGTCCTCCCACGCCCATCCCTTGTTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAAt                                                                            <  2:772417/100‑1 (MQ=255)
                        ggTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATc                                                                           <  1:1251378/100‑1 (MQ=255)
                          tGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCa                                                                         >  2:906931/1‑100 (MQ=255)
                          tGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCa                                                                         <  1:1464294/100‑1 (MQ=255)
                           gCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCa                                                                         >  1:1535605/1‑99 (MQ=255)
                           gCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCAt                                                                        <  2:3439122/100‑1 (MQ=255)
                           gCTCGATGCCAATGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCAt                                                                        >  1:1995032/1‑100 (MQ=255)
                            cTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCAtt                                                                       >  2:1110775/1‑100 (MQ=255)
                              cGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTtc                                                                     >  2:1033394/1‑100 (MQ=255)
                              cGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTtc                                                                     >  2:3725182/1‑100 (MQ=255)
                                 tGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTcc                                                                  >  2:541075/1‑100 (MQ=255)
                                 tGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTcc                                                                  <  1:1918845/100‑1 (MQ=255)
                                  gCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCt                                                                 >  2:1714621/1‑100 (MQ=255)
                                   ccAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGACGATGTGCATCAACGGTGAATCCATTTCTCCt                                                                 >  1:2372081/1‑99 (MQ=255)
                                    cAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCttt                                                               <  1:3799115/100‑1 (MQ=255)
                                     aGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTg                                                              <  1:264157/100‑1 (MQ=255)
                                          gTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGt                                                         >  2:3303483/1‑100 (MQ=255)
                                           tATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGtt                                                        >  1:661257/1‑100 (MQ=255)
                                            aTCAACCGTCCTCCCACGCGCATCCCTTGGTATTACCGTCTAACCGCCCGGGGGGCAATGACGTACCTCCACGGTAACTCCATTTCTCCTTTTaaa                                                           >  1:3178703/1‑95 (MQ=255)
                                             tCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGc                                                      <  1:650838/100‑1 (MQ=255)
                                              cAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGcc                                                     <  2:2410521/100‑1 (MQ=255)
                                              cAACCGTCCTCCCACGCCCATCCCTTGGTATTCCCGTATGACCGCCCGGGGGGCCAGCCTGTGCCTCAACGGTGAATCCCCTCATCCTTTGAAGGTTGcc                                                     >  1:912874/1‑100 (MQ=255)
                                                  cGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCggg                                                 >  2:1635859/1‑100 (MQ=255)
                                                    tcctccCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGta                                               <  2:721535/100‑1 (MQ=255)
                                                    tcctccCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGta                                               <  1:1269991/100‑1 (MQ=255)
                                                     cctccCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGCTGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGta                                               >  2:1774198/1‑99 (MQ=255)
                                                     cctccCACGCCCATCCCTTGGTATTCCGGTATGACCACCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGt                                                <  1:98185/98‑1 (MQ=255)
                                                       tccCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGtaaat                                            >  1:3787694/1‑100 (MQ=255)
                                                         ccACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGtaaataa                                          >  1:1085209/1‑100 (MQ=255)
                                                          cACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGtaaataa                                          >  2:3575372/1‑99 (MQ=255)
                                                           aCGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTaaataaaa                                        >  1:1385295/1‑100 (MQ=255)
                                                            cGCCCATCCCTTGGTATTCCGGTATGACCTCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTaaataaaat                                       <  1:2643796/100‑1 (MQ=255)
                                                              cccATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAatat                                     >  1:3735966/1‑100 (MQ=255)
                                                                cATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGc                                   <  1:2438932/100‑1 (MQ=255)
                                                                 aTCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGcc                                  <  1:3387521/100‑1 (MQ=255)
                                                                    ccTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGa                               >  2:11301/1‑100 (MQ=255)
                                                                    ccTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGa                               >  2:296383/1‑100 (MQ=255)
                                                                    ccTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGa                               >  2:2523565/1‑100 (MQ=255)
                                                                       tGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGACGGTTGCCCGGGTAAATAAAATATGCCCGa                               >  1:3239167/1‑97 (MQ=255)
                                                                        ggTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGc                           <  1:1318152/100‑1 (MQ=255)
                                                                        ggTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGc                           >  2:1017886/1‑100 (MQ=255)
                                                                            ttCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTg                       <  1:3277825/100‑1 (MQ=255)
                                                                              ccGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTACATAAAATATGCCCGACGGCTTTGcc                     >  1:826575/1‑100 (MQ=255)
                                                                              ccGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGcc                     <  1:3570560/100‑1 (MQ=255)
                                                                               cGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCt                    <  2:1134243/100‑1 (MQ=255)
                                                                               cGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGATGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCt                    >  1:2331057/1‑100 (MQ=255)
                                                                                ggTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCtt                   <  1:285913/100‑1 (MQ=255)
                                                                                ggTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCtt                   >  1:1545148/1‑100 (MQ=255)
                                                                                ggTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCtt                   <  2:2474669/100‑1 (MQ=255)
                                                                                 gTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTc                  >  2:1872320/1‑100 (MQ=255)
                                                                                 gTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTc                  <  1:765832/100‑1 (MQ=255)
                                                                                 gTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTc                  <  1:356370/100‑1 (MQ=255)
                                                                                 gTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTc                  <  2:1599708/100‑1 (MQ=255)
                                                                                 gTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTc                  <  1:100958/100‑1 (MQ=255)
                                                                                    tGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATa               >  2:837482/1‑100 (MQ=255)
                                                                                    tGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATa               <  1:1237249/100‑1 (MQ=255)
                                                                                    tGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATa               >  1:1405028/1‑100 (MQ=255)
                                                                                    tGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATa               >  2:1231746/1‑100 (MQ=255)
                                                                                     gACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATAc              <  2:2893310/100‑1 (MQ=255)
                                                                                      aCCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATCTGCCCGACGGCTTTGCCTTCATAcc             >  1:2593165/1‑100 (MQ=255)
                                                                                      aCCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATAcc             >  1:1991651/1‑100 (MQ=255)
                                                                                       ccGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATTCCCGACGGCTTTGCCTTCATACCg            >  1:509642/1‑100 (MQ=255)
                                                                                       ccGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCg            >  1:2564678/1‑100 (MQ=255)
                                                                                       cagCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCg            >  2:3522007/3‑100 (MQ=255)
                                                                                        cGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCgg           <  2:2255656/100‑1 (MQ=255)
                                                                                         gCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCggg          >  1:3617641/1‑100 (MQ=255)
                                                                                          ccAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCACACCGGGc         <  2:1988105/100‑1 (MQ=255)
                                                                                           cAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCt        <  2:3252309/100‑1 (MQ=255)
                                                                                            aGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCt        <  1:1171047/99‑1 (MQ=255)
                                                                                             tggggGCGAGGATGTGCAGCAACGGTTAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTag      <  1:1470314/99‑1 (MQ=255)
                                                                                             ggggggCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTag      >  2:1735519/1‑100 (MQ=255)
                                                                                                 aGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTagcagc  <  2:3279540/100‑1 (MQ=255)
                                                                                                 |                                                                                                   
CCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGCAGC  >  NC_000913/4533439‑4533635

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: