Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,008,473 +G intergenic (+40/‑193) b4555 → / → yigM hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,008,4721.G98.9% 1404.2 / 3.9 371intergenic (+39/‑194)b4555/yigMhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (2/2);  new base G (185/182);  total (187/184)
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

Only 100 of 379 total aligned reads displayed.

AACTGGCTACAGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGAAAGCTGCATTATGGGCAGTGGCCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGC  >  NC_000913/4008399‑4008561
                                                                                   |                                                                                           
actgggCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAg                                                                            <  2:6198175/92‑1 (MQ=255)
 ctgggCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGt                                                                           >  1:7196977/4‑96 (MQ=255)
  tgggCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTg                                                                          >  2:14020536/3‑96 (MQ=255)
   gggCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTgc                                                                         <  2:12041227/95‑2 (MQ=255)
   gggCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTgc                                                                         <  2:3477950/95‑2 (MQ=255)
    ggCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTgcg                                                                        >  1:7076684/1‑94 (MQ=255)
    ggCTACAGCGTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTgcg                                                                        <  2:13437271/96‑3 (MQ=255)
     gCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTgcgc                                                                       >  2:7685925/1‑93 (MQ=255)
      cTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTgcgca                                                                      >  2:198723/1‑92 (MQ=255)
       tACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTAt                                                                                  >  2:14250077/1‑83 (MQ=255)
         cagccTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCgcacca                                                                   <  1:14422890/92‑1 (MQ=255)
         cagccTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCgcacca                                                                   <  2:3572644/92‑1 (MQ=255)
         aagccTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCgcacca                                                                   >  1:11336053/5‑96 (MQ=255)
          agccTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAg                                                                  >  1:8056151/4‑96 (MQ=255)
           gccTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGc                                                                 >  2:8000064/3‑96 (MQ=255)
           gccTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGc                                                                 >  2:7169485/3‑96 (MQ=255)
           gccTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGc                                                                 <  2:9198338/94‑1 (MQ=255)
            ccTGAAGGATTGTATTGCGTTT‑‑GCTGACGGGATGAACGACGCCGAAATGCTGTCGATTGCGGGGAAAGCCTGCATTATGGGCAGTGCGCACCAGCg                                                                <  2:10184571/95‑1 (MQ=255)
            ccTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCg                                                                >  1:5733034/2‑96 (MQ=255)
            ccTGAAGGATTGTACTGCGTTTGGGGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCAGTATGGGCAGTGCGCACCAGCg                                                                >  2:4580742/2‑96 (MQ=255)
             cTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGt                                                               >  1:4717965/1‑96 (MQ=255)
             cTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGt                                                               >  2:9355629/1‑96 (MQ=255)
               tGAAGGTTTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCCGGCAAAGGCTGCATTATGGGCAGTGAGCACCAGCGTc                                                              >  1:10000040/1‑96 (MQ=255)
               tGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTc                                                              >  1:10854513/1‑96 (MQ=255)
               tGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTc                                                              >  2:12779655/1‑96 (MQ=255)
               tGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGGGCGCACCAGCGTc                                                              >  1:330447/1‑96 (MQ=255)
                 aaGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTg                                                            >  1:6964135/1‑96 (MQ=255)
                   ggATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGaa                                                          >  1:6141413/1‑96 (MQ=255)
                               tttt‑‑GTGACGGGATGAACGACGCCTAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATccc                                             <  1:13245395/94‑1 (MQ=255)
                                 tttggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGTGCAGTGAGCACCAGTGTCTGAAAGACCGTCATccc                                             >  1:5092274/5‑96 (MQ=255)
                                 tttggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATccc                                             <  2:9934253/92‑1 (MQ=255)
                                  ttggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCGGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCg                                            >  2:667754/4‑96 (MQ=255)
                                  ctggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCg                                            >  1:8869207/4‑96 (MQ=255)
                                   tggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGa                                           >  1:13187070/3‑96 (MQ=255)
                                     gtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGc                                         >  2:5019823/1‑96 (MQ=255)
                                       tgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCt                                        <  2:8539751/96‑1 (MQ=255)
                                       tgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCt                                        >  1:12160187/1‑96 (MQ=255)
                                        gACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTg                                      >  2:11592110/1‑96 (MQ=255)
                                        gACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTg                                      >  2:11053208/1‑96 (MQ=255)
                                        gACGCGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTg                                      <  1:2012447/96‑1 (MQ=255)
                                         aCGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTgg                                     >  1:12434257/1‑95 (MQ=255)
                                          cGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTgga                                    >  1:8685149/1‑94 (MQ=255)
                                           gggATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTggaa                                   <  1:9130859/96‑4 (MQ=255)
                                           gggATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTggaa                                   >  2:13587071/1‑93 (MQ=255)
                                           gggATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGAGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTggaa                                   >  2:6383718/1‑93 (MQ=255)
                                           gggATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGGGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTggaa                                   >  2:5675779/1‑93 (MQ=255)
                                              aTGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTg                                 >  2:1559585/1‑96 (MQ=255)
                                               tGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGa                                <  2:7039376/96‑1 (MQ=255)
                                                gAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGGCCTTCATCCCGAGCTGGAAGTGAt                               <  1:3650023/96‑1 (MQ=255)
                                                 aaCGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGAtt                              >  2:9511239/1‑96 (MQ=255)
                                                  acgacgCCGACATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTg                             <  1:10446427/96‑1 (MQ=255)
                                                  acgacgCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTg                             <  1:1342954/96‑1 (MQ=255)
                                                   cgacgCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTgg                            >  2:3437713/1‑96 (MQ=255)
                                                   cgacgCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTga                            <  1:13792932/96‑2 (MQ=255)
                                                    gccgCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGt                           <  2:13414290/94‑1 (MQ=255)
                                                    gacgCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGt                           >  1:3907276/1‑96 (MQ=255)
                                                     acgCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTa                          >  1:5201485/1‑96 (MQ=255)
                                                     acgCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTa                          >  1:6395801/1‑96 (MQ=255)
                                                       gCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACt                        >  2:327963/1‑96 (MQ=255)
                                                       gCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACt                        <  2:4288595/96‑1 (MQ=255)
                                                       gCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGGGATTGGTACt                        >  1:2883237/1‑96 (MQ=255)
                                                       gCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCACCCCGGGCTGGAAGTGATTGGTACt                        <  1:4551446/96‑1 (MQ=255)
                                                        ccGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTa                       <  2:583212/96‑1 (MQ=255)
                                                        ccGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTa                       <  1:9500971/96‑1 (MQ=255)
                                                         cGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAGGTGATTGGTACTaa                      >  1:5445168/1‑96 (MQ=255)
                                                          gAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAAt                     <  2:313102/96‑1 (MQ=255)
                                                          gAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGGCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAAt                     >  2:1055763/1‑96 (MQ=255)
                                                           aaaTGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCCCCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATg                    <  2:2593835/96‑1 (MQ=255)
                                                           aaaTGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATg                    >  2:9593158/1‑96 (MQ=255)
                                                           aaaTGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATg                    >  2:13125727/1‑96 (MQ=255)
                                                           aaaTGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATg                    <  1:11681880/96‑1 (MQ=255)
                                                           aaaTGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATg                    <  1:2478460/96‑1 (MQ=255)
                                                           aaaTGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATg                    >  1:4869010/1‑96 (MQ=255)
                                                             aTGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGcc                  <  1:3112074/96‑1 (MQ=255)
                                                             aTGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGcc                  <  2:11831566/96‑1 (MQ=255)
                                                               ggtGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcga                >  1:7521313/3‑96 (MQ=255)
                                                               gCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcga                <  2:11280995/96‑1 (MQ=255)
                                                                cTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGTCCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgac               <  1:14065070/96‑1 (MQ=255)
                                                                cTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgac               >  1:6755947/1‑96 (MQ=255)
                                                                cTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgac               <  1:9543973/96‑1 (MQ=255)
                                                                cTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgac               <  1:1902855/96‑1 (MQ=255)
                                                                cTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgac               >  2:4433339/1‑96 (MQ=255)
                                                                 tGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGGCCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgacn              <  1:4189712/96‑2 (MQ=255)
                                                                 tGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgacg              <  1:12779655/96‑1 (MQ=255)
                                                                 tGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgacg              <  2:7485192/96‑1 (MQ=255)
                                                                  gTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGCCCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgacga             <  2:2379647/96‑1 (MQ=255)
                                                                   tCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCgacgac            >  1:4409288/1‑96 (MQ=255)
                                                                    cGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGacgacg           <  2:4841377/96‑1 (MQ=255)
                                                                     gATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGAcgc          >  2:6557579/1‑96 (MQ=255)
                                                                     gATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGAcgc          >  2:4046019/1‑96 (MQ=255)
                                                                     gATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGAcgc          >  2:13001925/1‑96 (MQ=255)
                                                                     gATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGAcgc          <  1:11590987/96‑1 (MQ=255)
                                                                      aTGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGAcgcg         <  2:5226627/96‑1 (MQ=255)
                                                                      aTGGCGGGGAAAGTCTGCATTATGGGCAGTGCGCCCCAGCGTCTGAACGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGAcgcg         <  2:12519757/96‑1 (MQ=255)
                                                                         ggcggggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGt       >  1:9971213/5‑96 (MQ=255)
                                                                         ggcggggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGt       <  2:12289099/92‑1 (MQ=255)
                                                                            ggggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGcc    >  2:12841156/2‑96 (MQ=255)
                                                                            ggggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGcc    <  2:881401/95‑1 (MQ=255)
                                                                             gggAAAGACTGCATTATGGAGAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGATCTGGAAGTGATTGGTACTAATGCCGACGACTCGGTGCCg   >  1:6505205/1‑96 (MQ=255)
                                                                              ggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGc  <  1:14237029/96‑1 (MQ=255)
                                                                                   |                                                                                           
AACTGGCTACAGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGAAAGCTGCATTATGGGCAGTGGCCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGC  >  NC_000913/4008399‑4008561

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: