Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,019,857 C→T intergenic (+56/‑28) b3836 → / → b3838 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,019,8570CT99.4% 648.8 / NA 181intergenic (+56/‑28)b3836/b3838hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (108/71);  total (109/72)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 188 total aligned reads displayed.

AAGACGCGAAGCGCCACGATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGCCGCCTTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGG  >  NC_000913/4019767‑4019942
                                                                                          |                                                                                     
aagaCGCGAAGCGCCACGATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCtcgtcg                                                                                     >  2:1345716/1‑93 (MQ=255)
 agaCGCGAAGCGCCACGATAAAGAGCAGGTGTAATCCGTGTTTGATATCGCTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCtcgtcgt                                                                                    <  1:3294866/93‑2 (MQ=255)
   gcgcgAAGCGCCACGATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCtcgtcgttc                                                                                  <  1:171796/92‑4 (MQ=255)
         aGCGCCACGATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGc                                                                          >  1:2484788/1‑93 (MQ=255)
          gcgcCACGATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGcc                                                                         >  1:2823776/1‑93 (MQ=255)
            gcCACGATCAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGc                                                                       >  2:4563537/1‑93 (MQ=255)
            gcCACGATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGTCCTCGTCG‑‑TTCTGGcgccgc                                                                       >  1:3466660/1‑93 (MQ=255)
            gcCACGATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGc                                                                       >  2:4563538/1‑93 (MQ=255)
             ccACGATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCa                                                                      >  1:3003427/1‑93 (MQ=255)
              cACGATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCaa                                                                     <  1:375810/93‑1 (MQ=255)
              cACGATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCaa                                                                     <  1:800520/93‑1 (MQ=255)
              cACGATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCaa                                                                     >  2:2617117/1‑93 (MQ=255)
               aCGATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAAc                                                                    >  1:2546817/1‑93 (MQ=255)
               aCGATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAAc                                                                    <  2:3179537/93‑1 (MQ=255)
                 gATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGa                                                                  <  2:2746535/93‑1 (MQ=255)
                 gATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGa                                                                  >  2:4281093/1‑93 (MQ=255)
                 gATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGa                                                                  >  2:2935824/1‑93 (MQ=255)
                  aTAAAGAGCAGGTGTATTCCGTGTTTGATATCGGTTTTAGCTAACTGCTATTGGTGTTCATCGTCTGCGTCGTCG‑‑TTCTGGGGCCGCAACGAc                                                                 <  2:821913/93‑1 (MQ=255)
                    aaaGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTg                                                               >  1:3996292/1‑93 (MQ=255)
                     aaGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGc                                                              >  1:635342/1‑93 (MQ=255)
                      agagCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTGGTCG‑‑TTCTGGGGCCGCAACGACTGcc                                                             >  2:3795562/1‑93 (MQ=255)
                       gagCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCt                                                            <  2:2075046/93‑1 (MQ=255)
                           aGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTgg                                                        >  2:1909834/1‑93 (MQ=255)
                            ggTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTTTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGc                                                       <  1:164808/93‑1 (MQ=255)
                              tgtAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCgg                                                     >  2:116904/1‑93 (MQ=255)
                              tgtAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCgg                                                     >  1:3860048/1‑93 (MQ=255)
                                tAATCTGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTa                                                   <  1:16101/93‑1 (MQ=255)
                                tAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTa                                                   >  2:3849200/1‑93 (MQ=255)
                                tAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTa                                                   >  1:1115017/1‑93 (MQ=255)
                                 aaTCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGtca                                                  >  2:2107178/1‑91 (MQ=255)
                                 aaTCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTaa                                                  >  2:310454/1‑93 (MQ=255)
                                 aaTCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTaa                                                  >  1:584270/1‑93 (MQ=255)
                                 aaTCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTaa                                                  <  1:2589241/93‑1 (MQ=255)
                                 aaTCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTaa                                                  >  1:1340830/1‑93 (MQ=255)
                                  aTCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTATGGCGGTaaa                                                 >  2:2493355/1‑93 (MQ=255)
                                   tCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTaaaa                                                >  1:3716171/1‑93 (MQ=255)
                                   tCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTaaaa                                                >  2:3257889/1‑93 (MQ=255)
                                   tCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTaaaa                                                >  1:3936660/1‑93 (MQ=255)
                                    ccGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTaaaaa                                               <  2:531222/93‑1 (MQ=255)
                                     cGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAAc                                              <  1:22436/93‑1 (MQ=255)
                                       tgtTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACgg                                            >  1:120249/1‑93 (MQ=255)
                                        gtTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTTATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACgga                                           >  2:1856164/1‑92 (MQ=255)
                                        gtTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGt                                           >  2:4051428/1‑93 (MQ=255)
                                         tttGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTa                                          >  1:2934545/1‑93 (MQ=255)
                                         tttGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTa                                          >  2:4479703/1‑93 (MQ=255)
                                         tttGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGGAAAAACGGTa                                          >  2:194513/1‑93 (MQ=255)
                                          ttGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAg                                         <  1:1508986/93‑1 (MQ=255)
                                           tGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGc                                        >  1:1535952/1‑93 (MQ=255)
                                            gATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCg                                       >  1:709983/1‑93 (MQ=255)
                                            gATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCg                                       >  1:1411182/1‑93 (MQ=255)
                                             atatCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCgg                                      <  1:775316/93‑1 (MQ=255)
                                             atatCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCgg                                      >  2:575507/1‑93 (MQ=255)
                                              tatCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCggg                                     >  2:740322/1‑93 (MQ=255)
                                              tatCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGAggg                                     >  2:1507088/1‑93 (MQ=255)
                                                tCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCG‑‑GCGTTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCt                                   >  1:4570590/1‑93 (MQ=255)
                                                  ggTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTgg                                 >  1:542154/1‑93 (MQ=255)
                                                   gTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGa                                <  2:2272141/93‑1 (MQ=255)
                                                    ttttAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGAt                               >  1:3082827/1‑93 (MQ=255)
                                                     tttAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGAtt                              <  1:1801467/93‑1 (MQ=255)
                                                      ttAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTc                             >  1:4245752/1‑93 (MQ=255)
                                                       tAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTcg                            <  1:2630523/93‑1 (MQ=255)
                                                        aGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTcgc                           <  2:4596527/93‑1 (MQ=255)
                                                        aGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTcgc                           >  1:4430773/1‑93 (MQ=255)
                                                         gCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTcgcg                          >  2:2164619/1‑93 (MQ=255)
                                                         gCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTcgcg                          >  1:4364093/1‑93 (MQ=255)
                                                         gCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTcgcg                          >  1:4400933/1‑93 (MQ=255)
                                                          cGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTcgcgc                         <  2:358477/93‑1 (MQ=255)
                                                           gAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTcgcgcg                        >  2:1863471/1‑93 (MQ=255)
                                                           gAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTcgcgcg                        >  2:1030554/1‑93 (MQ=255)
                                                            aaCTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCgcgt                       >  1:188473/1‑93 (MQ=255)
                                                             aCTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCgcgtt                      >  2:1378666/1‑93 (MQ=255)
                                                             aCTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCgcgtt                      >  2:2434776/1‑93 (MQ=255)
                                                                  tATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCgcgttgcgtt                 >  2:2683804/1‑93 (MQ=255)
                                                                   aTTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTc                <  1:3658354/93‑1 (MQ=255)
                                                                   aTTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTc                >  2:580563/1‑93 (MQ=255)
                                                                    ttGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCa               <  2:1139800/93‑1 (MQ=255)
                                                                     tGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCAc              <  2:3768118/93‑1 (MQ=255)
                                                                      ggTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACt             >  2:1559617/1‑93 (MQ=255)
                                                                      ggTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACt             <  2:2563243/93‑1 (MQ=255)
                                                                      ggTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACt             <  1:4476928/93‑1 (MQ=255)
                                                                      ggTGTGCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACt             <  2:1826852/93‑1 (MQ=255)
                                                                       gtgtTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTg            <  1:316559/93‑1 (MQ=255)
                                                                        tgtTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTgg           >  1:1323598/1‑93 (MQ=255)
                                                                         gtTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGc          >  2:1516521/1‑93 (MQ=255)
                                                                         gtTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGc          >  1:812450/1‑93 (MQ=255)
                                                                          ttCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCg         <  2:2264298/93‑1 (MQ=255)
                                                                          ttCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCg         >  1:3889641/1‑93 (MQ=255)
                                                                           tcatcaTCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGa        <  2:3753908/93‑1 (MQ=255)
                                                                           tcatcaTCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGa        <  1:2437061/93‑1 (MQ=255)
                                                                           tcatcaGCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTAGATTCGCGCGTTGCGTTCACTGGCGa        <  1:1260328/93‑1 (MQ=255)
                                                                             atcatcGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACa      >  2:3551388/1‑93 (MQ=255)
                                                                             atcatcGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACa      <  1:1000836/93‑1 (MQ=255)
                                                                              tcatcGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACaa     >  1:1046716/1‑93 (MQ=255)
                                                                              tcatcGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACaa     >  2:2018624/1‑93 (MQ=255)
                                                                              tcatcGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACaa     <  1:2201066/93‑1 (MQ=255)
                                                                              tcatcGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCCGGCTGGATTCGCGCGTTGCGTTCACTGGCGACaa     >  1:2022595/1‑93 (MQ=255)
                                                                               catcGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAAc    >  1:461963/1‑93 (MQ=255)
                                                                                atcGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACg   <  2:4439070/93‑1 (MQ=255)
                                                                                atcGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACg   >  1:478812/1‑93 (MQ=255)
                                                                                 tcGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCCGGCTGGATTTGCGCGTTGCGCTCACTGGCGACaaacg  >  2:94829/1‑90 (MQ=255)
                                                                                          |                                                                                     
AAGACGCGAAGCGCCACGATAAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGCCGCCTTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGG  >  NC_000913/4019767‑4019942

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: