Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,038,117 C→G 16.8% H126Q (CAC→CAG cbdA → cytochrome bd‑II oxidase, subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,038,1170CG16.8% 243.9 / 35.8 91H126Q (CAC→CAGcbdAcytochrome bd‑II oxidase, subunit I
Reads supporting (aligned to +/- strand):  ref base C (44/31);  new base G (10/5);  total (54/37)
Fisher's exact test for biased strand distribution p-value = 7.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTT  >  NC_000913/1038019‑1038216
                                                                                                  |                                                                                                   
gcACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCAcc                                                                                                    >  2:2428867/1‑100 (MQ=255)
gcACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCAcc                                                                                                    >  1:658438/1‑100 (MQ=255)
gcACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCAcc                                                                                                    >  2:2000311/1‑100 (MQ=255)
   ccGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACctgc                                                                                                 <  1:679331/100‑1 (MQ=255)
      cTGGCGATTGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTgg                                                                                              <  1:993018/100‑1 (MQ=255)
      cTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTgg                                                                                              <  1:2049902/100‑1 (MQ=255)
      cTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTgg                                                                                              <  2:733299/100‑1 (MQ=255)
       tGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTgg                                                                                              <  2:2115521/99‑1 (MQ=255)
       tGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGt                                                                                             <  1:1929036/100‑1 (MQ=255)
       tGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTGCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGt                                                                                             <  2:1137816/100‑1 (MQ=255)
        ggCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTg                                                                                            >  1:1270438/1‑100 (MQ=255)
        ggCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTg                                                                                            >  1:1615008/1‑100 (MQ=255)
        ggCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTg                                                                                            >  1:187527/1‑100 (MQ=255)
        ggCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTg                                                                                            <  1:2220448/100‑1 (MQ=255)
         gCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCAGGTGCTGGtat                                                                                           <  1:1227637/100‑3 (MQ=255)
          cGATGGAAGCATTAATGGCCTTCTTCCTCGAAGCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCAGGTGCTGGtatt                                                                                          <  2:1425634/100‑4 (MQ=37)
           gATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCAGGTGCTGGtattt                                                                                         <  2:1473085/100‑5 (MQ=37)
           gATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCAGGTGCTGGtattt                                                                                         <  1:1706894/100‑5 (MQ=37)
            aTGGAAGCATTAATGGCTTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTGGGCTGGCAACGGCTGAATAAATACCAGCAGGTGCTGGtattta                                                                                        <  2:178616/100‑6 (MQ=37)
             tGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGtn                                                                                       <  2:1562911/100‑2 (MQ=255)
                ccgatctAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCt                                                                                    <  2:535744/94‑1 (MQ=255)
                aaGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCt                                                                                    >  1:1057998/1‑100 (MQ=255)
                aaGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCt                                                                                    >  1:2119117/1‑100 (MQ=255)
                    aTTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCtggtg                                                                                <  1:2484184/100‑1 (MQ=255)
                    aTTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGATGAATAAATACCAGCACCTGCTGGTGACGTGGCtggtg                                                                                >  2:729236/1‑100 (MQ=255)
                     ttAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCtggtgg                                                                               <  2:1036593/100‑1 (MQ=255)
                      tAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGc                                                                              <  1:638201/100‑1 (MQ=255)
                       aaTGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCtagatcgg                                                                             >  1:535751/1‑93 (MQ=255)
                       aaTGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCg                                                                             >  2:2089850/1‑100 (MQ=255)
                        aTGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGt                                                                            <  1:2232308/100‑1 (MQ=255)
                          ggCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTc                                                                          <  2:2119103/100‑1 (MQ=255)
                           gCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCg                                                                         >  2:775006/1‑100 (MQ=255)
                            ccTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCgg                                                                        >  1:411671/1‑100 (MQ=255)
                               tcttcCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTc                                                                     <  1:744372/100‑1 (MQ=255)
                                cttcCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCa                                                                    >  2:701035/1‑100 (MQ=255)
                                 ttcCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGACTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCaa                                                                   <  1:2009625/100‑1 (MQ=255)
                                  tcCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCaaa                                                                  >  2:338264/1‑100 (MQ=255)
                                   ccTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAAt                                                                 <  1:59660/100‑1 (MQ=255)
                                   ccTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAAt                                                                 >  2:2336533/1‑100 (MQ=255)
                                   ccTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGCCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGtt                                                                      >  1:2254567/1‑95 (MQ=255)
                                    cTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAAtc                                                                >  1:2491056/1‑100 (MQ=255)
                                    cTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAAtc                                                                >  2:1893675/1‑100 (MQ=255)
                                    cTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAAtc                                                                >  1:1039232/1‑100 (MQ=255)
                                    cTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAAtc                                                                >  1:98138/1‑100 (MQ=255)
                                     tCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAAtc                                                                <  2:200578/99‑1 (MQ=255)
                                          tCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTg                                                          <  1:124200/100‑1 (MQ=255)
                                           ccACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTg                                                          <  2:1624056/99‑1 (MQ=255)
                                            cACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTgn                                                         <  2:571060/99‑2 (MQ=255)
                                             aCCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGt                                                       >  2:384977/1‑100 (MQ=255)
                                              ccTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGtt                                                      >  1:809410/1‑100 (MQ=255)
                                              ccTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTATCTGCGtt                                                      >  2:2420950/1‑100 (MQ=255)
                                                tttGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTtgt                                                    >  1:1295839/1‑100 (MQ=255)
                                                   gTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAATTACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCGCTCTGCGTTGTggn                                                 <  2:1396208/100‑2 (MQ=255)
                                                   gTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTgg                                                  >  2:1729668/1‑99 (MQ=255)
                                                           nttcttcttcGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAAt                                         >  1:2158079/2‑100 (MQ=255)
                                                           gTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAAt                                         <  1:627478/100‑1 (MQ=255)
                                                            ttcttcttcGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATg                                        <  2:230568/100‑1 (MQ=255)
                                                             tgttcttcGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGc                                       >  1:1733924/3‑100 (MQ=255)
                                                               ttcttcGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCa                                     >  2:695928/1‑100 (MQ=255)
                                                                  ttcGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTCGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTATCTGCGTTGTGGATATTGAATGACAACg                                  >  1:1221452/1‑100 (MQ=255)
                                                                   tcGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACgg                                 >  1:182227/1‑100 (MQ=255)
                                                                     ggCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGCTATTGCATgaataaggtt                               >  1:852267/1‑91 (MQ=255)
                                                                     ggCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACgg                                 >  2:2419842/1‑98 (MQ=255)
                                                                      gCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTg                              <  1:2451895/100‑1 (MQ=255)
                                                                        tGGCAACGGCTCAATACATACCAGCACCTGTTGGTAACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGa                            <  1:1566991/100‑1 (MQ=255)
                                                                         ggCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGa                            >  1:2166846/1‑99 (MQ=255)
                                                                           cAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGc                         >  1:984261/1‑100 (MQ=255)
                                                                           cAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGc                         <  2:2075993/100‑1 (MQ=255)
                                                                           cAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGc                         <  1:1416917/100‑1 (MQ=255)
                                                                            aaCGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCa                        <  1:1129530/100‑1 (MQ=255)
                                                                              cGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAAt                      >  2:150743/1‑100 (MQ=255)
                                                                                 cTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATAcc                   >  1:1664300/1‑100 (MQ=255)
                                                                                 cTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATAcc                   >  2:384954/1‑100 (MQ=255)
                                                                                   gccgcgtcaCCAGCAGGAGCTGGTGACGTGGCTGATGGGGTTCGGTTCAAATCTCTCGGCATTATAGATATGGAATGGCAACGGTGGGAAGCAATAc                    >  2:1763144/9‑97 (MQ=255)
                                                                                   gccacgtcaCCAGCAGGTGCTGGTGACGTGGGTGGTGGCGTCCGGTTCAAATCTCGCTGCGTTGTGTATTTTGAATGCCGACGGTTGGATGCAATAccc                  >  2:852707/9‑99 (MQ=255)
                                                                                   gAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCg                 >  1:2464117/1‑100 (MQ=255)
                                                                                    aataaataCCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGa                >  2:2388886/1‑100 (MQ=255)
                                                                                    aataaataCCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGa                >  2:453290/1‑100 (MQ=255)
                                                                                     cacgtcaCCAGCAGGTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGAc               >  1:418693/7‑100 (MQ=255)
                                                                                     ataaataCCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCACTACCCGa                >  2:1169433/1‑99 (MQ=255)
                                                                                       cgtcaCCAGCAGGTGCTTGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCaa                       >  2:969232/5‑90 (MQ=255)
                                                                                       cgtcaCCAGCAGGTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCg             >  2:152469/5‑100 (MQ=255)
                                                                                        gtcaCCAGCAGGTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCgg            >  1:311655/4‑100 (MQ=255)
                                                                                        gtcaCCAGCAGGTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCgg            >  1:95032/4‑100 (MQ=255)
                                                                                        gtcaCCAGCAGGTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCgg            >  1:2002524/4‑100 (MQ=255)
                                                                                        aataCCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCgg            >  1:2530305/1‑100 (MQ=255)
                                                                                        aataCCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACACGACCgg            >  2:2060054/1‑100 (MQ=255)
                                                                                         tcaCCAGCAGGTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGt           >  2:431702/3‑100 (MQ=255)
                                                                                         ataCCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGt           >  2:2297852/1‑100 (MQ=255)
                                                                                          caCCAGCAGGTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTg          >  2:1389797/2‑100 (MQ=255)
                                                                                            ccAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGGTGGATGCAATACCCGACCGGTgcg        >  2:599645/1‑100 (MQ=255)
                                                                                             cAGCATCTGCTGGTGGCGTGGCTGGTGGCGTTCGGTTCAAATCTCGCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTgcg        <  1:632993/99‑1 (MQ=255)
                                                                                             cAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTgcgc       <  2:654951/100‑1 (MQ=255)
                                                                                              aGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCa      >  2:609157/1‑100 (MQ=255)
                                                                                                  ccTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCAtttt  <  1:2011383/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
GCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTT  >  NC_000913/1038019‑1038216

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: