Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,333,582:1 (A)8→9 100% intergenic (‑131/‑128) rplU ← / → ispB 50S ribosomal subunit protein L21/octaprenyl diphosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,333,5741.A87.5% 303.0 / 39.6 104intergenic (‑123/‑136)rplU/ispB50S ribosomal subunit protein L21/octaprenyl diphosphate synthase
Reads supporting (aligned to +/- strand):  ref base . (6/7);  new base A (48/43);  total (54/50)
Fisher's exact test for biased strand distribution p-value = 7.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 125 total aligned reads displayed.

CCTTTTCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAG‑‑AAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTG  >  NC_000913/3333475‑3333668
                                                                                                      |                                                                                                
ccTTTTCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACgaag                                                                                                   >  2:2520232/1‑100 (MQ=255)
ccTTTTCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACgaag                                                                                                   <  1:183548/100‑1 (MQ=255)
ccTTTTCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACgaag                                                                                                   >  1:3624291/1‑100 (MQ=255)
  ttttCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACgaag‑‑aa                                                                                               <  1:3861843/100‑1 (MQ=255)
    ttCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTTCTACAGTCAATACACGAAG‑‑aaaa                                                                                             >  2:2040889/1‑100 (MQ=255)
    ttCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAG‑‑aaaa                                                                                             >  2:1390375/1‑100 (MQ=255)
       aaTGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAG‑‑aaaaaaa                                                                                         <  1:1738330/100‑1 (MQ=255)
        aTGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAG‑‑aaaaaaaa                                                                                        >  2:3866461/1‑100 (MQ=255)
        aTGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAG‑‑aaaaaaaa                                                                                        >  1:3927991/1‑100 (MQ=255)
          gATTCAGAGTGCGCTATAAATATTCACATTAGGGCGTGAATATTACGAAAAACGAACGCCTTGGACAAGTGCTACAGTCAATACACGAAG‑‑aaaaaaaaac                                                                                      <  1:3324515/100‑3 (MQ=255)
          gATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAG‑‑aaaaaaaaac                                                                                      <  1:822143/100‑3 (MQ=255)
          gATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCAAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAG‑‑aaaaaaaaac                                                                                      <  2:3558697/100‑3 (MQ=255)
           aTTCAGAGTGCGCTATAAATATTCACAATAGGGTGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAG‑‑aaaaaaaaac                                                                                      <  1:3022738/99‑3 (MQ=255)
            ttCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACgaag‑‑aagatcggaaga                                                                                     >  1:2634133/1‑90 (MQ=255)
             tCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAG‑‑aaaaaaaaacaca                                                                                   <  1:1983472/100‑6 (MQ=255)
             tCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAG‑‑aaaaaaaaacaca                                                                                   <  1:1955768/100‑6 (MQ=255)
              cAGAGTGCGCTATAAATATTCACAATAGTGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAG‑‑AAAAAAAAaacaca                                                                                  >  2:1193116/1‑98 (MQ=255)
              cAGAGTGCGCTATAAATATTCACAATAGGGCGCGACTATTACGCAAAACGCACGCCTTTGAAAAGTGATACAGTCAATACACGAAG‑‑aaaaaaaaacaca                                                                                   >  1:3485167/1‑94 (MQ=255)
              cAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTAACAAATTCTACAGTCAATACACGAAG‑‑aaaaaaaaacacca                                                                                  >  1:1325928/1‑94 (MQ=255)
               agagTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAAc                                                                                  <  1:2443196/100‑1 (MQ=255)
               agagTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAAc                                                                                  <  2:1362454/100‑1 (MQ=255)
                  gTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTg                                                                               <  1:958776/100‑1 (MQ=255)
                   tGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTAATACAGACAATACACGAAGAAAAAAAAACACAACTTg                                                                               >  1:1076262/1‑99 (MQ=255)
                    tcttccgatctatTCACAATAGGGCGCGAATATTACGAAAAACACACGCCTTTGACACGTGCTACAGTCAATACACGAAAAAAAAAAAACACAACTTGTa                                                                             <  2:3752428/90‑1 (MQ=255)
                    gcgcTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTa                                                                             >  2:571543/1‑100 (MQ=255)
                     cgcTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTAc                                                                            >  2:1979423/1‑100 (MQ=255)
                      gcTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTAAATAAAAGAAGAAAAAAAAACACAACTTGTACg                                                                           <  1:677041/100‑1 (MQ=255)
                       cTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTAc                                                                            >  1:2481552/1‑98 (MQ=255)
                       cTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACgg                                                                          >  2:3283456/1‑100 (MQ=255)
                        tataAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACgg                                                                          <  1:1957598/99‑1 (MQ=255)
                         ataAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAAGACAACTTGTACGATa                                                                        >  1:3870953/1‑100 (MQ=255)
                           aaaTATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACAAAACTTGTACgtgaac                                                                      >  2:384909/1‑95 (MQ=255)
                             atatTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACACCTTGTACGGTAACGt                                                                    >  2:3589655/1‑100 (MQ=255)
                              tatggaaaaTAGGGCGCTAATATTACGTAAAACGTACGCTTTTGATAAGTGCTATAGTTTATAAAAAAAAAAAAAAAAACACAACTTGTACGGTAACGt                                                                    <  2:3149295/92‑1 (MQ=255)
                                  cacaATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACAAAAA‑‑AAAAAAAACACAACTTGTACGGTAACGTTTa                                                                 >  1:1049192/1‑98 (MQ=255)
                                  cacaATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGtaaat                                                                >  1:2585363/1‑95 (MQ=255)
                                  cacaATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATc                                                               <  1:2750607/100‑1 (MQ=255)
                                  cacaATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACAAGAAGAAAAAAAAAACACAACTTGTACGGTAACGTTGAt                                                                >  2:691617/1‑100 (MQ=255)
                                    caATAGGGCGCGAATATTACGCAAAACGCACGCCTTTTACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCt                                                              >  1:266155/1‑99 (MQ=255)
                                        aGGGCGCTAATATTACGTAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGcc                                                         <  1:938525/100‑1 (MQ=255)
                                        aGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGcc                                                         >  2:2668572/1‑100 (MQ=255)
                                         gggCGCGAATATTACGCAAAACGCACGCCTTTGACAATTGCTACAGTTAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGcc                                                         <  1:3619037/99‑1 (MQ=255)
                                         gggCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCa                                                        >  1:3389100/1‑100 (MQ=255)
                                            cgcgAATATTACGCAAAACGCACTCCTTTGACAATTTCTACAGTCAATACACGAAAA‑‑AAAAAAAACACAACTTATACGGAAACGTTTATCTGTGCCAtt                                                      >  1:132303/1‑99 (MQ=255)
                                            cgcgAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCAttt                                                     >  2:3920959/1‑100 (MQ=255)
                                              cgAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACATTTATCTGTGCCAttttt                                                   >  2:1728108/1‑100 (MQ=255)
                                               gAATATTACGCAAAACGCTCGCCTTTGACAAGTGCTACAGTCTTTACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCAtttttn                                                  <  1:567308/100‑2 (MQ=255)
                                               gAATATTACGCAAAACGCACGACTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCAttttt                                                   <  2:1860313/99‑1 (MQ=255)
                                                aataatACGTAAAACACACTCATTTGACAAGTACTACAGTCAATACACGAAAA‑‑AAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGACATTTTTTc                                                 <  2:3489429/100‑1 (MQ=255)
                                                aaTATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTc                                                 <  1:3843610/100‑1 (MQ=255)
                                                aaTATTACGCAAAACGCAAGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTc                                                 <  1:2933635/100‑1 (MQ=255)
                                                aaTATTACGAAAAATGCACGCTTTTGTCAAGTGTTACAGTTAATACATGAAGAAAAAAAAACACAACTGGTACGGTAACGTTTATCTGTGCCATTTTTTc                                                 <  1:1580577/100‑1 (MQ=255)
                                                 atatTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCa                                                <  2:3624252/100‑1 (MQ=255)
                                                 atatTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCa                                                <  2:2785233/100‑1 (MQ=255)
                                                 atatTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCa                                                <  2:3389062/100‑1 (MQ=255)
                                                  tatTACGTAAAACGCACGCCTTTGACAATTGTTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAg                                               <  1:3587483/100‑1 (MQ=255)
                                                  tatTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTTAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAg                                               <  1:3123476/100‑1 (MQ=255)
                                                  tatTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGTCAtttttt                                                  >  2:1741085/1‑97 (MQ=255)
                                                  tatTAAGCAAAACGCACGCCTTTTACAAGTGTTATTGTTAATAAATGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAg                                               <  1:861663/100‑1 (MQ=255)
                                                   atTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAG‑‑AAAAAACAACACAACTTGTACGTTAACTTTTATCTGTGCCAtttttt                                                  >  2:43700/1‑96 (MQ=255)
                                                   atTACGCAAAACGCACGCCTTTGACAAGTGCGACAGTAAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGt                                              <  2:1695546/100‑1 (MQ=255)
                                                       cGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTCTTTGTGCAATTTTTTCAGTACaa                                          >  2:1696296/1‑100 (MQ=255)
                                                        gCAAATCGCACGCCTTTGACAAGTGTTACTGTTAATACATGTAG‑‑AAAAAATAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAAt                                         <  1:499765/100‑1 (MQ=255)
                                                        gCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAAt                                         >  2:1885208/1‑100 (MQ=255)
                                                         cAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTTTACGATAACGTTTATCTTTGCCATTTTTTCAGTACAAt                                         >  1:3251482/1‑99 (MQ=255)
                                                         cAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTAAACAGAAACAAATAACTGTGCCATTTTTTCAAAACaa                                          >  1:397698/1‑98 (MQ=255)
                                                         cAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAAAACAACTTTTACCGTAACGCTTATCTGTACCATTTTTTTAGttctact                                        >  2:2123319/1‑94 (MQ=255)
                                                          aaaaCGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACAGAAACATATATCTGTGCCATTTTTTCATTACAATCa                                       >  1:1028110/1‑100 (MQ=255)
                                                           aaaCGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTAAGTAAAATCAc                                      >  2:258046/1‑100 (MQ=255)
                                                           aaaCGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGATAACGTTTATCTGTGCCATTTTTTGAGTACAATCa                                       >  2:1156505/1‑99 (MQ=255)
                                                                cacgcCTTTGACAAGTGCTACAGTCAATACCCGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCAAAGCCATTTTTTAATTACAATCAACCt                                   >  1:868145/1‑98 (MQ=255)
                                                                 acgcCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGCAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCtat                                 >  2:1681134/1‑99 (MQ=255)
                                                                  ggcCTTTGACAAATGCTACAGGCAATACATGAAGAAAAAAAAACACAACTTCTACAGAAAAAATAATATGTGCCATTTTTTCAGTAAAATCAAACTatat                               >  2:1269042/2‑100 (MQ=255)
                                                                    ccTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACTGTAACGTTTATCTATACCAAATTTTAACTACAATAACAATATAtt                              >  1:1712007/1‑99 (MQ=255)
                                                                      tttGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTcc                            <  1:2529473/99‑1 (MQ=255)
                                                                       ttGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTa                          <  2:2869868/100‑1 (MQ=255)
                                                                         gACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGTTAACGTTAAACAATACCAATTTTTCAAAACAATCACCCTATATTCATaa                         >  1:1810925/1‑99 (MQ=255)
                                                                         gACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTaan                        <  1:1050644/100‑2 (MQ=255)
                                                                         gACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAAc                        <  2:1247302/100‑1 (MQ=255)
                                                                         gACAAGTGCTACAGTCAATACACGAAG‑‑AAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAAcc                       <  1:177307/100‑1 (MQ=255)
                                                                          aCAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCATTTTAATTACCCTATATTCCTAAcc                       >  1:1387217/1‑100 (MQ=255)
                                                                          aCAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAAcc                       <  1:1455431/100‑1 (MQ=255)
                                                                          aCAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAAcc                       <  2:2720752/100‑1 (MQ=255)
                                                                           cAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACCGTAACGTTTATCAGTGCCATTTTTTCATTATAATCAACCTATAtttata                          >  2:107047/1‑92 (MQ=255)
                                                                            aaagacTACAGTTAATAAACAAAAA‑‑AAAAAAAACACAACTTGTACTGTAACGTTTATCTGTGCTATTTTTTCAGTACAATCACCCTATATTCCTAACCa                      <  2:3501372/94‑1 (MQ=255)
                                                                              gTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTAAAATAAAATCACCCTATATTCCTAACCATaa                   >  1:1690814/1‑100 (MQ=255)
                                                                              gTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTACCATTTTTTCAGTACAATCACCCTATATTCCTAATCATaa                   >  2:413591/1‑100 (MQ=255)
                                                                                gCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTTCTAACCATAAAc                 >  1:3133362/1‑100 (MQ=255)
                                                                                gCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAAc                 <  1:3814559/100‑1 (MQ=255)
                                                                                  tACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAAccc               <  2:724479/100‑1 (MQ=255)
                                                                                  tACAGTCAATACACGAAG‑‑AAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCt              >  2:387133/1‑100 (MQ=255)
                                                                                        cAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGtt         >  1:1331842/1‑100 (MQ=255)
                                                                                         aaTAAACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGt          <  2:2928401/98‑1 (MQ=255)
                                                                                           tACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGcc      <  2:2204357/100‑1 (MQ=255)
                                                                                           tACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATAACCCTATATTCCTAACCATAAACCCTAAGTTGcc      >  1:1525330/1‑100 (MQ=255)
                                                                                             cacGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCtt    >  2:749156/1‑100 (MQ=255)
                                                                                             cacGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCAATTTTTTCAGTACAATCACTCTATATTCATAACTATAAACCCTAAGTTGCCtt    >  2:870588/1‑100 (MQ=255)
                                                                                               cGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCtt    >  2:845311/1‑98 (MQ=255)
                                                                                               cGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTg  <  1:3292909/100‑1 (MQ=255)
                                                                                               cGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAATACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGACtt    >  2:215611/1‑98 (MQ=255)
                                                                                                      |                                                                                                
CCTTTTCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAG‑‑AAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTG  >  NC_000913/3333475‑3333668

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