Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,008,491 2 bp→CG intergenic (+58/‑174) 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,4910GC99.5% 1600.4 / NA 435intergenic (+58/‑175)b4555/yigMhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base C (194/238);  total (196/239)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009134,008,4920CG100.0% 1612.1 / NA 428intergenic (+59/‑174)b4555/yigMhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (195/233);  total (195/233)

Only 100 of 462 total aligned reads displayed.

GGCTACAGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGG‑‑AAAGCTGCATTATGGGCAGTGGCCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTA  >  NC_000913/4008403‑4008580
                                                                                                   ||                                                                                          
    cagccTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCgcacca                                                                                       <  2:6322331/92‑1 (MQ=255)
     agccTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAg                                                                                      <  1:9290859/93‑1 (MQ=255)
      gccTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGc                                                                                     >  1:3561769/3‑96 (MQ=255)
      gccTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGCAAGGCTGCATTATGGGCAGTGCGCACCAGc                                                                                     >  2:3006521/3‑96 (MQ=255)
       ccTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCg                                                                                    >  1:5733034/2‑96 (MQ=255)
       ccTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCg                                                                                    >  1:4988502/2‑96 (MQ=255)
        cTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGt                                                                                   >  2:9355629/1‑96 (MQ=255)
        cTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGCGCGCACCAGCGt                                                                                   >  1:11378010/1‑96 (MQ=255)
        cTGAAGGAATGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGt                                                                                   >  1:6490533/1‑96 (MQ=255)
         tGAAGGTTTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCCGGC‑‑AAAGGCTGCATTATGGGCAGTGAGCACCAGCGTc                                                                                  >  1:10000040/1‑96 (MQ=255)
         tGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTc                                                                                  >  1:10854513/1‑96 (MQ=255)
             ggATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGaa                                                                              >  1:6141413/1‑96 (MQ=255)
                          tttggtgACGGGATGAACGCCGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCTCCAGCGTCTGAAAGACCTTCATccc                                                                 <  2:3506327/92‑1 (MQ=255)
                          tttggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATccc                                                                 <  2:4703266/92‑1 (MQ=255)
                           ttggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCg                                                                <  1:14114330/93‑1 (MQ=255)
                           ttggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCGGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCg                                                                >  2:667754/4‑96 (MQ=255)
                           ctggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCg                                                                >  1:8869207/4‑96 (MQ=255)
                            tggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATcccta                                                               <  1:2979521/94‑3 (MQ=255)
                            tggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGa                                                               >  1:13187070/3‑96 (MQ=255)
                              gtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGc                                                             <  1:7038455/96‑1 (MQ=255)
                               tgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCt                                                           <  2:4543996/96‑1 (MQ=255)
                               tgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGTG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCATCGTCTGCAAGACCTTCATCCCGAGCt                                                           >  1:6961764/1‑96 (MQ=255)
                                gACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTg                                                         >  1:12337876/1‑96 (MQ=255)
                                  cGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTgga                                                       >  1:2747296/1‑94 (MQ=255)
                                   gggATGAACGAGGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTggaa                                                      >  2:14237029/1‑93 (MQ=255)
                                   gggATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTggaa                                                      <  1:9130859/96‑4 (MQ=255)
                                   gggATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTggaa                                                      <  2:10508543/96‑4 (MQ=255)
                                   gggATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTggaa                                                      >  2:2391646/1‑93 (MQ=255)
                                        gAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGAt                                                  >  2:9542865/1‑96 (MQ=255)
                                          acgacgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTg                                                <  1:13336172/96‑1 (MQ=255)
                                           cgacgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTgg                                               >  2:3437713/1‑96 (MQ=255)
                                           cgacgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTga                                               <  1:13792932/96‑2 (MQ=255)
                                            gccgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGt                                              <  2:13414290/94‑1 (MQ=255)
                                            gacgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGt                                              <  1:524372/96‑1 (MQ=255)
                                            gacgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGt                                              >  2:13565387/1‑96 (MQ=255)
                                             acgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTa                                             >  1:6395801/1‑96 (MQ=255)
                                             acgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTa                                             >  2:198522/1‑96 (MQ=255)
                                               gCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACt                                           <  2:4288595/96‑1 (MQ=255)
                                               gCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACt                                           <  1:1999619/96‑1 (MQ=255)
                                               gCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCACCCCGGGCTGGAAGTGATTGGTACt                                           <  1:4551446/96‑1 (MQ=255)
                                                ccGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGATGTGATTGGTACTa                                          >  2:4844077/1‑96 (MQ=255)
                                                 cGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAGGTGATTGGTACTaa                                         >  1:5445168/1‑96 (MQ=255)
                                                 cGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTaa                                         <  2:5827692/96‑1 (MQ=255)
                                                 cGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCAATATGGGAAGTGCGCACAAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGtaataa                                         >  1:98990/1‑96 (MQ=255)
                                                 cGAAATGCTGTCGAGGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTaa                                         >  1:5114308/1‑96 (MQ=255)
                                                   aaaTGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATg                                       >  2:13125727/1‑96 (MQ=255)
                                                   aaaTGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATg                                       >  1:4869010/1‑96 (MQ=255)
                                                   aaaTGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATg                                       <  1:11681880/96‑1 (MQ=255)
                                                     aTGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGcc                                     <  2:11831566/96‑1 (MQ=255)
                                                     aTGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGcc                                     <  2:1212676/96‑1 (MQ=255)
                                                     aTGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGcc                                     <  1:7967824/96‑1 (MQ=255)
                                                      tGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcg                                    <  1:10492689/96‑1 (MQ=255)
                                                      tGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcg                                    <  2:8611170/96‑1 (MQ=255)
                                                       ggtGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcga                                   >  1:7521313/3‑96 (MQ=255)
                                                       gCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcga                                   <  2:784704/96‑1 (MQ=255)
                                                        cTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgac                                  <  1:1902855/96‑1 (MQ=255)
                                                        cTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgac                                  >  1:6755947/1‑96 (MQ=255)
                                                        cTGTCGATGCAGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgac                                  <  1:2393129/96‑1 (MQ=255)
                                                         tGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGGGATTGGGACTAATGCcgacg                                 >  2:4132733/1‑96 (MQ=255)
                                                          gTCGATGGGGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgacga                                >  1:4671930/1‑96 (MQ=255)
                                                           tCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCgacgac                               <  2:10162357/96‑1 (MQ=255)
                                                           tCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCgacgac                               >  1:4409288/1‑96 (MQ=255)
                                                             gATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGAcgc                             >  2:9173504/1‑96 (MQ=255)
                                                              aTGGCGGGG‑‑AAAGGCTGCATTATTGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGGGCTGGAAGTGATTGGTACTAATGCCGACGAcgcg                            >  2:12291817/1‑96 (MQ=255)
                                                              aTGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGAcgcg                            >  2:8359228/1‑96 (MQ=255)
                                                              aTGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAAACCTTCATCCCGAGCTGGAAGGGATTGGTACTAATGGCGGCGAcgcg                            >  1:3960574/1‑96 (MQ=255)
                                                              aTGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGCCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGAcgcg                            >  2:7091006/1‑96 (MQ=255)
                                                              aTGGCGGGG‑‑AAAGTCTGCATTATGGGCAGTGCGCCCCAGCGTCTGAACGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGAcgcg                            <  2:12519757/96‑1 (MQ=255)
                                                                     gggg‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGcc                       >  2:7465530/2‑96 (MQ=255)
                                                                     gggg‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGcc                       >  1:3224311/2‑96 (MQ=255)
                                                                       gg‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCTTCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGATGCCGc                     >  2:12783958/1‑96 (MQ=255)
                                                                       gg‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGc                     >  2:3275642/1‑96 (MQ=255)
                                                                            aaaggCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCAt                   >  1:3852446/5‑96 (MQ=255)
                                                                            aaaggCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCAt                   >  2:10913418/5‑96 (MQ=255)
                                                                             aaggCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGCTTGGTACTAATGCCGACGACGCGGTGCCGCAtt                  <  2:14066289/93‑1 (MQ=255)
                                                                             aaggCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCAtt                  >  2:2782277/4‑96 (MQ=255)
                                                                              aggCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTa                 <  2:5166531/94‑1 (MQ=255)
                                                                              aggCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTa                 >  2:9738272/3‑96 (MQ=255)
                                                                                gCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATc               <  2:6141413/96‑1 (MQ=255)
                                                                                 tcTGCATTATGGGCAGTGCGCACCAGCGACTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATc               <  1:12248005/95‑1 (MQ=255)
                                                                                  cTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCt              >  1:9627120/1‑96 (MQ=255)
                                                                                  cTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCt              >  1:7229577/1‑96 (MQ=255)
                                                                                   tGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTg             <  2:11461145/96‑1 (MQ=255)
                                                                                   tGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTg             <  2:6507775/96‑1 (MQ=255)
                                                                                    gCATTATGGGCAGTGCGCCCCAGCGTCTGAAAGACCTTCATCCCGTGCGGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGc            <  1:4245751/96‑1 (MQ=255)
                                                                                    gCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGc            <  2:9200765/96‑1 (MQ=255)
                                                                                     cATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCg           <  1:12681184/96‑1 (MQ=255)
                                                                                      aTTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGt          <  2:5005219/96‑1 (MQ=255)
                                                                                      aTTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGt          <  2:5482252/96‑1 (MQ=255)
                                                                                       ttATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTa         <  1:11990528/96‑1 (MQ=255)
                                                                                         aTGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATTTGCGTaaa       <  2:12971223/96‑1 (MQ=255)
                                                                                         aTGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTaaa       >  2:11164189/1‑96 (MQ=255)
                                                                                         aTGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTaaa       <  1:1425508/96‑1 (MQ=255)
                                                                                          tGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAAc      >  1:8222077/1‑96 (MQ=255)
                                                                                           gggCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGCCGACGCGGTGCGGCATTATGTGCGTaagtg     >  1:10237775/1‑93 (MQ=255)
                                                                                           gggCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAAct     <  2:3907276/96‑1 (MQ=255)
                                                                                           gggCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAAct     >  2:1678771/1‑96 (MQ=255)
                                                                                            ggCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAActc    >  2:682961/1‑96 (MQ=255)
                                                                                            ggCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAActc    <  2:9621549/96‑1 (MQ=255)
                                                                                             gCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAActct   >  2:760669/1‑96 (MQ=255)
                                                                                                   ||                                                                                          
GGCTACAGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGG‑‑AAAGCTGCATTATGGGCAGTGGCCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTA  >  NC_000913/4008403‑4008580

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