Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,331,223 (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,331,2151.A97.6% 927.4 / 10.4 248intergenic (‑123/‑136)rplU/ispB50S ribosomal subunit protein L21/octaprenyl diphosphate synthase
Reads supporting (aligned to +/- strand):  ref base . (4/2);  new base A (136/106);  total (140/108)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 291 total aligned reads displayed.

AATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAG‑‑AAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTT  >  NC_000913/3331123‑3331302
                                                                                             |                                                                                          
aaTGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACgaag                                                                                             >  2:1561860/1‑93 (MQ=255)
aaTGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACgaag                                                                                             >  2:4653276/1‑93 (MQ=255)
 aTGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACgaag‑‑a                                                                                          <  1:4964883/93‑1 (MQ=255)
 aTGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACgaag‑‑a                                                                                          >  1:4130627/1‑93 (MQ=255)
  tGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACgaag‑‑aa                                                                                        >  1:3999773/1‑93 (MQ=255)
     ttCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAG‑‑aaaaa                                                                                     >  1:5534983/1‑93 (MQ=255)
      tCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAG‑‑aaaaaa                                                                                    <  1:1703856/93‑1 (MQ=255)
      tCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAG‑‑aaaaaa                                                                                    >  2:4451942/1‑93 (MQ=255)
       cAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAG‑‑aaaaaaa                                                                                  >  2:5217064/1‑93 (MQ=255)
       cAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAG‑‑aaaaaaa                                                                                  <  1:15784/93‑1 (MQ=255)
        agagTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAG‑‑aaaaaaaa                                                                                 >  2:2099953/1‑93 (MQ=255)
         gagTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACTAAG‑‑aaaaaaaaa                                                                                <  2:1359713/93‑2 (MQ=255)
         gagTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAG‑‑aaaaaaaaa                                                                                >  1:2855238/1‑92 (MQ=255)
         gagTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAG‑‑aaaaaaaaa                                                                                >  2:6139519/1‑92 (MQ=255)
         gagTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAG‑‑aaaaaaaaa                                                                                <  2:3787638/93‑2 (MQ=255)
          agTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAG‑‑aaaaaaaaac                                                                               <  1:1156751/93‑3 (MQ=255)
              cgcTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACaa                                                                            >  2:3202507/1‑93 (MQ=255)
                cTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACt                                                                          >  1:1760818/1‑93 (MQ=255)
                 tataAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACtt                                                                         <  1:2805905/93‑1 (MQ=255)
                 tataAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACtt                                                                         >  1:5901238/1‑93 (MQ=255)
                  ataAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTg                                                                        >  1:5731843/1‑93 (MQ=255)
                   taAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTTAATACACGAAGAAAAAAAAACACAACTTGt                                                                       >  2:3938583/1‑93 (MQ=255)
                   taAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGt                                                                       >  1:3223447/1‑93 (MQ=255)
                   taAATATTCACAATAGGGCGCAAATATTACGCAAAACGCACGCCTTTGACAAGGGCTACAGTCAAAACCCGAAAAAAAAAAAAACACAACTTg                                                                        <  2:5844242/93‑1 (MQ=255)
                    aaaTATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTa                                                                      <  2:3387062/93‑1 (MQ=255)
                    aaaTATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAG‑‑AAAAAACAACACAACTTGtg                                                                      >  2:6438313/1‑92 (MQ=255)
                     aaTTTTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTAc                                                                     >  1:4707025/1‑93 (MQ=255)
                       tatTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACgg                                                                   >  2:4048607/1‑93 (MQ=255)
                         ttCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTa                                                                 >  1:2000690/1‑93 (MQ=255)
                         ttCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAAAACAACTTGTACGGTa                                                                 >  2:4500840/1‑93 (MQ=255)
                          tCACAATAGGGCGCGAATATTACGCAGAACCGACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGtca                                                                <  1:6616239/93‑3 (MQ=255)
                           cacaATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAAc                                                               <  1:2396795/93‑1 (MQ=255)
                           cacaATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAAc                                                               >  2:5939085/1‑93 (MQ=255)
                            acaATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACg                                                              >  2:5798069/1‑93 (MQ=255)
                            acaATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACg                                                              >  2:3971808/1‑93 (MQ=255)
                            acaATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACg                                                              >  1:2071471/1‑93 (MQ=255)
                            acaATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACg                                                              >  2:4835483/1‑93 (MQ=255)
                              aaTAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGtt                                                            >  2:5532114/1‑93 (MQ=255)
                               aTAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGttt                                                           >  1:6498501/1‑93 (MQ=255)
                                tAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTa                                                          <  2:5650570/93‑1 (MQ=255)
                                 aGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTAt                                                         >  1:1014299/1‑93 (MQ=255)
                                  gggCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATc                                                        >  2:885403/1‑93 (MQ=255)
                                  gggCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATc                                                        >  1:6031133/1‑93 (MQ=255)
                                  gggCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATc                                                        <  1:5113591/93‑1 (MQ=255)
                                     cgcgAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCtgt                                                     <  2:2894174/93‑1 (MQ=255)
                                      gcgAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACCCGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCtgtg                                                    <  2:2647356/93‑1 (MQ=255)
                                      gcgAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCtgtg                                                    <  1:5229510/93‑1 (MQ=255)
                                        gAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGcc                                                  >  1:5779910/1‑93 (MQ=255)
                                        gAATATTACGCAAAACGCAAGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGcc                                                  >  1:6571161/1‑93 (MQ=255)
                                        aaaTATTACACAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGcc                                                  <  2:3620000/92‑1 (MQ=255)
                                         aaTATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCa                                                 >  2:1463451/1‑93 (MQ=255)
                                          atatTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCAt                                                >  2:6502308/1‑93 (MQ=255)
                                          atatTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCAt                                                <  2:3096895/93‑1 (MQ=255)
                                          atatTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCAt                                                <  2:1807651/93‑1 (MQ=255)
                                           tatTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCAtt                                               >  1:4269462/1‑93 (MQ=255)
                                           aatTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCAtt                                               <  2:3813840/92‑1 (MQ=255)
                                            atTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCAttt                                              <  1:3812022/93‑1 (MQ=255)
                                            atTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCAttt                                              <  2:255790/93‑1 (MQ=255)
                                             ttACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCAtttt                                             <  2:712644/93‑1 (MQ=255)
                                             ttACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCAtttt                                             >  2:4778014/1‑93 (MQ=255)
                                             ttACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCAtttt                                             >  2:5642918/1‑93 (MQ=255)
                                              tACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCAttttt                                            <  1:690051/93‑1 (MQ=255)
                                                 gCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCa                                         <  1:1552692/93‑1 (MQ=255)
                                                  cAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAg                                        >  1:132718/1‑93 (MQ=255)
                                                   aaaaCGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGt                                       >  1:4453673/1‑93 (MQ=255)
                                                   aaaaCGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGt                                       >  2:2340811/1‑93 (MQ=255)
                                                      acgcacgcCTTTGACAAGTGCTACAGTCAATACACGAAG‑‑ACAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACa                                    >  2:1861411/1‑93 (MQ=255)
                                                        gcacgcCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAAt                                  >  2:4362129/1‑93 (MQ=255)
                                                         cacgcCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATc                                 <  2:1712350/93‑1 (MQ=255)
                                                         cacgcCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATc                                 >  2:3737286/1‑93 (MQ=255)
                                                           cgcCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCAc                               >  2:2718859/1‑93 (MQ=255)
                                                             ccTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCAccc                             >  1:300376/1‑93 (MQ=255)
                                                             ccTTTGACAAGTGCTACAGTCAATACACAAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCAccc                             <  1:4923021/93‑1 (MQ=255)
                                                              cTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCt                            >  1:2630502/1‑93 (MQ=255)
                                                               tttGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCta                           <  2:6637290/93‑1 (MQ=255)
                                                               tttGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCta                           >  1:5349418/1‑93 (MQ=255)
                                                               tttGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCta                           <  2:2797261/93‑1 (MQ=255)
                                                                ttGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCtat                          <  2:2050282/93‑1 (MQ=255)
                                                                 tGACAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCtata                         <  1:6449523/93‑1 (MQ=255)
                                                                    cAAGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTc                      >  1:4386377/1‑93 (MQ=255)
                                                                     aaGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTcc                     >  2:5027404/1‑93 (MQ=255)
                                                                      aGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCt                    >  1:6126506/1‑93 (MQ=255)
                                                                      aGTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCt                    >  2:559536/1‑93 (MQ=255)
                                                                       gTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTa                   <  1:3333246/93‑1 (MQ=255)
                                                                       gTGCTACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTa                   >  1:4661369/1‑93 (MQ=255)
                                                                           tACAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCa               <  1:1463286/93‑1 (MQ=255)
                                                                            aCAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCAt              <  1:6144716/93‑1 (MQ=255)
                                                                            aCAGTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCAt              >  2:6350823/1‑93 (MQ=255)
                                                                               gTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATaaa           <  2:123764/93‑1 (MQ=255)
                                                                               gTCAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATaaa           <  1:5545238/93‑1 (MQ=255)
                                                                                 cAATACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAAcc         <  1:6254616/93‑1 (MQ=255)
                                                                                  aaTACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAAccc        <  2:3550991/93‑1 (MQ=255)
                                                                                  aaTACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAAccc        <  1:4790367/93‑1 (MQ=255)
                                                                                  aaTACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAAccc        <  1:244035/93‑1 (MQ=255)
                                                                                   aTACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCt       >  1:789185/1‑93 (MQ=255)
                                                                                   aTACACGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAAACCt       >  1:809277/1‑93 (MQ=255)
                                                                                     acacGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTaa     <  2:2951850/93‑1 (MQ=255)
                                                                                      cacGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAg    >  2:2410945/1‑93 (MQ=255)
                                                                                      cacGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAg    >  1:3502845/1‑93 (MQ=255)
                                                                                      cacGAAGAAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAg    >  1:1098347/1‑93 (MQ=255)
                                                                                             |                                                                                          
AATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAG‑‑AAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTT  >  NC_000913/3331123‑3331302

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