Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,333,582 (A)8→9 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.A97.8% 295.9 / ‑0.5 89intergenic (‑123/‑136)rplU/ispB50S ribosomal subunit protein L21/octaprenyl diphosphate synthase
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base A (46/41);  total (47/42)
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

CCTTTTCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAG‑‑AAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTG  >  NC_000913/3333475‑3333668
                                                                                                    |                                                                                               
ccTTTTCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACgaag                                                                                                  >  1:503322/1‑100 (MQ=255)
ccTTTTCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACgaag                                                                                                  <  1:657981/100‑1 (MQ=255)
 cTTTTCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACgaag‑‑a                                                                                               <  2:2709648/100‑1 (MQ=255)
     tCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAG‑‑aaaaa                                                                                           >  1:1839702/1‑100 (MQ=255)
      cAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAG‑‑aaaaaa                                                                                          >  2:894075/1‑100 (MQ=255)
        aTGATTCAGAGTGCGCTATAAATATTAACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAG‑‑aaaaaaaa                                                                                        <  2:365884/100‑1 (MQ=255)
          gATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAG‑‑aaaaaaaaac                                                                                      >  1:2368765/1‑98 (MQ=255)
           aTTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAG‑‑aaaaaaaaaca                                                                                     >  2:202106/1‑97 (MQ=255)
           aTTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAG‑‑aaaaaaaaaca                                                                                     >  2:85953/1‑97 (MQ=255)
             tCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAG‑‑aaaaaaaaacaca                                                                                   >  2:533661/1‑95 (MQ=255)
             tCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAG‑‑AAAAAAAACACaa                                                                                   >  1:616199/1‑100 (MQ=255)
                gagTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACt                                                                                 <  2:1079516/100‑1 (MQ=255)
                    gcgcTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTa                                                                             >  1:2604521/1‑100 (MQ=255)
                      gcTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTAc                                                                            <  2:2957715/99‑1 (MQ=255)
                      gcTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACg                                                                           >  1:289332/1‑100 (MQ=255)
                      gatataAATATTCACAATAGGGCGCGAATATTACCCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACg                                                                           <  2:2271716/98‑1 (MQ=255)
                        tataAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTATAc                                                                            >  2:1875575/1‑97 (MQ=255)
                          taAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTaa                                                                       >  2:2415309/1‑100 (MQ=255)
                          taAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAAACGAAG‑‑AAAAAAAACACAACTTGTACGGTAAc                                                                      <  1:1272484/100‑1 (MQ=255)
                            aaTATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACg                                                                     <  1:280932/100‑1 (MQ=255)
                                  ctaaaTAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATc                                                               <  2:2241143/97‑1 (MQ=255)
                                  cacaATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGttt                                                                  >  2:1438631/1‑97 (MQ=255)
                                  cacaATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTAt                                                                >  1:735100/1‑99 (MQ=255)
                                  cacaATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACAAAAAAAAAAAAACACAACTTGTACGGTAACGTTTATc                                                               <  2:434924/100‑1 (MQ=255)
                                   acaATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTAtgt                                                              <  1:2095093/100‑3 (MQ=255)
                                   acaATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATc                                                               <  2:719486/99‑1 (MQ=255)
                                   acaATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCt                                                              >  2:3003256/1‑100 (MQ=255)
                                   acaAAAGGGCGCGAATATTACGCAAAAAGCACGCCTTTGACAAGTGCTACAGTCAATAAAAAAAAAAAAAAAAACACAACTTGTACGGTAACGTTTa                                                                 <  2:1032724/97‑1 (MQ=255)
                                    caATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTaacag                                                             >  1:2241166/1‑96 (MQ=255)
                                    caATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTAt                                                                >  2:2587323/1‑97 (MQ=255)
                                    caATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCtg                                                             <  1:582667/100‑1 (MQ=255)
                                    caATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTAATCtg                                                             >  2:2873016/1‑100 (MQ=255)
                                      aTAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCtgtg                                                           >  1:1416666/1‑100 (MQ=255)
                                         gggCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGcc                                                         >  2:1479737/1‑100 (MQ=255)
                                            cgcgAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCAttt                                                     >  1:2312667/1‑100 (MQ=255)
                                            cgcgAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCAttt                                                     >  2:1369251/1‑100 (MQ=255)
                                            cgcgAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCCGTGCCAttt                                                     >  1:2470108/1‑100 (MQ=255)
                                             gcgAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCAtttt                                                    <  1:2545118/100‑1 (MQ=255)
                                             gcgAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTCTACGGTAACGTTTATCTGTGCCAtttt                                                    >  1:2016065/1‑100 (MQ=255)
                                              cgAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCAttttt                                                   >  2:1307545/1‑100 (MQ=255)
                                                aaTATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCAttttttn                                                 <  1:2798029/100‑2 (MQ=255)
                                                aaTATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTc                                                 <  1:1687073/100‑1 (MQ=255)
                                                aaTATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACAAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTc                                                 <  2:28295/100‑1 (MQ=255)
                                                 atatTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCa                                                <  2:1548944/100‑1 (MQ=255)
                                                 atatTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTATACCATTTTTTCa                                                >  2:351089/1‑100 (MQ=255)
                                                 atatTACGAAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCa                                                <  1:120489/100‑1 (MQ=255)
                                                  tatTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAg                                               <  1:2315433/100‑1 (MQ=255)
                                                  tatTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGCAACGTTTATCTGTGCCATTTTTTCAg                                               <  1:1881739/100‑1 (MQ=255)
                                                   atTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTAACGTAAAGTTTCTCCTTGCCATTTTTTCAGt                                              >  1:3030208/1‑100 (MQ=255)
                                                    ttACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGt                                              >  2:1700319/1‑99 (MQ=255)
                                                     tACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTAc                                            <  2:626618/100‑1 (MQ=255)
                                                     tACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTAc                                            <  1:2821977/100‑1 (MQ=255)
                                                      aCGCAAAACGCATGCCTTTGTTAAGTGCTTCAGTCAATACGCGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAg                                               <  2:958721/96‑1 (MQ=255)
                                                       cGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACaa                                          >  2:3071054/1‑100 (MQ=255)
                                                         cAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAAt                                         >  1:659792/1‑99 (MQ=255)
                                                         cAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAAt                                         >  1:2568688/1‑99 (MQ=255)
                                                           aaaCGCACGCCTTTGACAAGTGCTACATTTAATAAAAAAAAAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATc                                        <  2:1877162/98‑1 (MQ=255)
                                                             acgcacgcCTTTGACAAGTGCTAAAGTCAAAACACGAAAAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCAccc                                    <  1:1283900/100‑1 (MQ=255)
                                                               gcacgcCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCta                                  <  2:77350/100‑1 (MQ=255)
                                                                cacgcCTTTGACAAGTGCTACAGTCAGTACACGAAGAAAAAAAAACACAACTTGTACGGCAACGCATAAATGTGCCATTTTTTCAGTACAATcccaat                                   >  1:2360216/1‑93 (MQ=255)
                                                                    ccTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTATGCCATTTTTTCAGTACAATCACCCTATATTc                             >  2:869957/1‑100 (MQ=255)
                                                                     cTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTcc                            >  2:1464390/1‑100 (MQ=255)
                                                                     cTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTcc                            <  1:202107/100‑1 (MQ=255)
                                                                     cTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTcc                            <  1:2415332/100‑1 (MQ=255)
                                                                      tttGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCt                           <  1:2606292/100‑1 (MQ=255)
                                                                      tttGACAAGTGCTACAGTCAATAAACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTcc                            <  2:3043365/99‑1 (MQ=255)
                                                                       ttGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCt                           <  1:1479751/100‑1 (MQ=255)
                                                                       ttGAAAAGTGCTACAGTAAATACACAAAGAAAAAAAAACACAACTTATACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTcc                            <  2:3073185/98‑1 (MQ=255)
                                                                        tGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTaa                         <  2:1427489/100‑1 (MQ=255)
                                                                        tGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTa                          >  2:1077326/1‑99 (MQ=255)
                                                                         gACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAAc                        <  2:1839682/100‑1 (MQ=255)
                                                                         gACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAAc                        >  1:2584698/1‑100 (MQ=255)
                                                                         gACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTATGCCATTTTTTCAGTACAATCACCCTATATTCCTAAc                        >  2:1335904/1‑100 (MQ=255)
                                                                         gACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGCCACGTTTATCTGTGCCATTTTTTCAGTACACTCTCCCTATATTCCTaa                         >  1:787338/1‑99 (MQ=255)
                                                                           cAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTAAGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCa                      >  2:2358409/1‑100 (MQ=255)
                                                                           cAAGTGCTACAGTCAATACACGAAAAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCa                      <  2:2157207/100‑1 (MQ=255)
                                                                            acttGCTACATTCAATACACAATGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAAcc                       <  2:1065149/95‑1 (MQ=255)
                                                                              gTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATaa                   >  2:2091634/1‑100 (MQ=255)
                                                                              gTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATaa                   >  2:2183240/1‑100 (MQ=255)
                                                                                 cTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTACCCATAAAc                 >  1:490215/1‑99 (MQ=255)
                                                                                 cTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAAcc                >  2:1201674/1‑100 (MQ=255)
                                                                                  tACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAATACAATCACCCTATATTCCTAACCATAAAccc               >  1:3078815/1‑100 (MQ=255)
                                                                                  tACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATATTTTCAGTACAATCACCCTATATTCCTAACCATAAAccc               >  2:1718215/1‑100 (MQ=255)
                                                                                   aCAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCt              <  2:2368742/100‑1 (MQ=255)
                                                                                    cAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTa             <  2:2365386/100‑1 (MQ=255)
                                                                                    cAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTa             >  1:2229887/1‑100 (MQ=255)
                                                                                      gTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTaa            >  1:1413862/1‑99 (MQ=255)
                                                                                      gTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAg           <  1:894084/100‑1 (MQ=255)
                                                                                        cAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGtt         >  2:329853/1‑100 (MQ=255)
                                                                                        cAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGtt         >  1:2470089/1‑100 (MQ=255)
                                                                                         aaTACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGt          <  2:3071583/98‑1 (MQ=255)
                                                                                          aTACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGc       <  1:1235067/100‑1 (MQ=255)
                                                                                           tACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGcc      >  2:17474/1‑100 (MQ=255)
                                                                                           tACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGcc      <  1:1473163/100‑1 (MQ=255)
                                                                                             cacGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCtt    >  1:2913010/1‑100 (MQ=255)
                                                                                               cGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTg  >  2:1845026/1‑100 (MQ=255)
                                                                                               cGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTg  >  2:2992383/1‑100 (MQ=255)
                                                                                               cGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTg  <  1:2515082/100‑1 (MQ=255)
                                                                                               cGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTg  >  2:879510/1‑100 (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 < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: