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.A89.2% 315.1 / 31.6 102intergenic (‑123/‑136)rplU/ispB50S ribosomal subunit protein L21/octaprenyl diphosphate synthase
Reads supporting (aligned to +/- strand):  ref base . (4/7);  new base A (42/49);  total (46/56)
Fisher's exact test for biased strand distribution p-value = 7.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 118 total aligned reads displayed.

CTTTTCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAG‑AAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTGTTC  >  NC_000913/3333476‑3333671
                                                                                                     |                                                                                                  
 ttttCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACg‑aag‑aa                                                                                                  >  2:1523539/1‑100 (MQ=255)
  tttCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAG‑aaa                                                                                                 >  1:3385839/1‑100 (MQ=255)
   ttCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAG‑aaaa                                                                                                >  2:89749/1‑100 (MQ=255)
    tCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAG‑aaaaa                                                                                               <  2:412372/100‑1 (MQ=255)
     cAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAG‑aaaaaa                                                                                              <  2:2367349/100‑1 (MQ=255)
     cAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAG‑aaaaaa                                                                                              >  1:1954043/1‑100 (MQ=255)
     cAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAG‑aaaaaa                                                                                              >  1:2619299/1‑100 (MQ=255)
      aaTGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAG‑aaaaaaa                                                                                             >  2:3067917/1‑100 (MQ=255)
      aaTGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACA‑AAG‑aaaaaaa                                                                                             <  1:955656/100‑1 (MQ=255)
      aaTGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCGTTGACAAGTGCTACAGTCAATAC‑ACG‑AAG‑aaaaaaa                                                                                             >  2:2663144/1‑100 (MQ=255)
       aTGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAG‑aaaaaaaa                                                                                            <  2:2504022/100‑1 (MQ=255)
        tGATTCAAAGTGCGCTATAAATATTCACAATAGGGCGCAAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACA‑AAG‑aaaaaaaa                                                                                            <  1:3800263/99‑1 (MQ=255)
          aTTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAG‑aaaaaaaaaca                                                                                         >  1:372149/1‑97 (MQ=255)
          aTTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAG‑AAAAAAAa‑aac                                                                                        <  2:278822/100‑1 (MQ=255)
             cAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACaa                                                                                      >  1:2117593/1‑100 (MQ=255)
                 gTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACtt                                                                                   >  1:2697955/1‑99 (MQ=255)
                    cgcTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTAc                                                                               <  2:3352905/100‑1 (MQ=255)
                      cTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTa                                                                                >  1:2057087/1‑97 (MQ=255)
                           aaTATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACg                                                                        <  2:1563196/100‑1 (MQ=255)
                           aaTATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACg                                                                        <  2:944217/100‑1 (MQ=255)
                             tatTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGtt                                                                      <  2:228015/100‑1 (MQ=255)
                              atTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGttt                                                                     >  2:3593000/1‑100 (MQ=255)
                              atTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGttt                                                                     >  1:2839358/1‑100 (MQ=255)
                               ttCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑AACAACTTGTACGTTAACGt                                                                       >  2:3611101/1‑97 (MQ=255)
                               ntCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTa                                                                    >  1:720731/2‑100 (MQ=255)
                                tCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGttt                                                                     >  1:2221504/1‑98 (MQ=255)
                                tCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAG‑AAAAAAAA‑CACAACTTGTACGGTAACGTTTATc                                                                  <  2:3147344/100‑1 (MQ=255)
                                 cacaATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGttt                                                                     >  1:3498122/1‑97 (MQ=255)
                                 cacaATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATc                                                                  >  2:3051499/1‑100 (MQ=255)
                                 cacaATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATc                                                                  >  2:3133257/1‑100 (MQ=255)
                                  acaATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTTACATTTATc                                                                  >  2:765469/1‑99 (MQ=255)
                                    aaTAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCt                                                                 <  1:631970/98‑1 (MQ=255)
                                    aaTAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAG‑AAAAAAAA‑CACAACTTGTACGGTAACGTTTATCtgtg                                                              >  2:1681452/1‑100 (MQ=255)
                                     aTAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGCCAACGTTTATc                                                                  >  1:1698613/1‑96 (MQ=255)
                                         ggCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGGAACGTTTATCTGTGCCAt                                                          >  2:3570960/1‑100 (MQ=255)
                                          gcgcGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCtgt                                                               >  1:3323826/1‑94 (MQ=255)
                                          gcgcGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCAtt                                                         >  2:3468975/1‑100 (MQ=255)
                                           cgcgAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTTCCAttt                                                        >  2:334716/1‑100 (MQ=255)
                                           cgcgAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCCGTGCCAttt                                                        >  2:406203/1‑100 (MQ=255)
                                            gcgAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCAtttt                                                       >  1:1847207/1‑100 (MQ=255)
                                              gAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCAtttttt                                                     >  2:2232353/1‑100 (MQ=255)
                                              gAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCAtttttt                                                     <  2:1602013/100‑1 (MQ=255)
                                               aaTATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTc                                                    <  1:3484451/100‑1 (MQ=255)
                                                ccgatctgCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAG‑AAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAg                                                  <  2:949522/93‑1 (MQ=255)
                                                atatTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTc                                                    <  1:2571882/99‑1 (MQ=255)
                                                atatTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCa                                                   <  2:2478356/100‑1 (MQ=255)
                                                atatTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTCATCTGTGCCATTTTTTCa                                                   >  2:1904030/1‑100 (MQ=255)
                                                atatTACGCAAAAAGCACGCCTTTGACAAGTGCTACTGTCAATAC‑AAA‑AAAAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTc                                                    <  1:858742/99‑1 (MQ=255)
                                                 tatTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCa                                                   >  2:564877/1‑99 (MQ=255)
                                                 tatTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAg                                                  <  2:1098608/100‑1 (MQ=255)
                                                 tatTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAAAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAg                                                  <  1:1051148/100‑1 (MQ=255)
                                                  atTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ATG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGt                                                 <  2:2255032/100‑1 (MQ=255)
                                                  atTACGCAAAACGAAAGACTTTGACAAGTGCTACAGTCAATAA‑AAAAAAA‑AAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAg                                                  <  1:336640/99‑1 (MQ=255)
                                                   ttACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTa                                                <  1:1983409/100‑1 (MQ=255)
                                                    tACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGt                                                 >  1:850178/1‑98 (MQ=255)
                                                       gCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAG‑AAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAgagatcgg                                           >  1:949517/1‑93 (MQ=255)
                                                         aaaaCGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCa                                          >  2:3777921/1‑100 (MQ=255)
                                                          aaaCGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCAc                                         >  2:1618534/1‑100 (MQ=255)
                                                           aaCGCACGCCTTTGACAAGTGCTACAGTCAATAC‑AAG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCAcc                                        <  1:1149494/100‑1 (MQ=255)
                                                             cgcacgcCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTAAGTACAATCAcac                                       >  1:1588144/1‑97 (MQ=255)
                                                             aaaacaaatttGACAAGTGTTACAGTCAAAAAAAAA‑AAA‑AAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCAccc                                       <  1:1667168/91‑1 (MQ=255)
                                                              gcacgcCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCta                                     <  2:1461039/100‑1 (MQ=255)
                                                               cacgcCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAG‑AAAAAAAACCACAACTTGTACGGTAACGTTTATCTGTGACATTTTTCCAGTACACTCACCCta                                     >  1:1023686/1‑99 (MQ=255)
                                                                aggcCTTAGACAAGTGCTACCGGCAAAAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCta                                     <  1:3062587/96‑1 (MQ=255)
                                                                   ccTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTCCGGTAACGTTTATCTGTGCCATTTTTTCCGTACAATCACCCTATAtt                                 >  2:1044475/1‑99 (MQ=255)
                                                                   ccTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTc                                >  2:1629813/1‑100 (MQ=255)
                                                                    cTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTcc                               <  2:1199732/100‑1 (MQ=255)
                                                                    cTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTcc                               <  2:1065144/100‑1 (MQ=255)
                                                                    cTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTcc                               <  2:3473697/100‑1 (MQ=255)
                                                                     tttGACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCt                              <  2:3083192/100‑1 (MQ=255)
                                                                     tttCACAAGTGTGACAGTCAATAA‑ACA‑AAAAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTcc                               <  1:1757641/99‑1 (MQ=255)
                                                                      ttGACAAGTGCTACAGTCAATAC‑ACA‑AAAAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTa                             <  1:2573253/100‑1 (MQ=255)
                                                                      ttGACAAGTGCTACAGTAAAAAC‑AAA‑AAAAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTcc                               <  1:1935111/98‑1 (MQ=255)
                                                                       tGACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTaa                            <  2:3241796/100‑1 (MQ=255)
                                                                        gACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTTAGTATATTCACCCTATATTCCTa                             >  1:265329/1‑98 (MQ=255)
                                                                        gACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAAc                           <  2:481580/100‑1 (MQ=255)
                                                                        gACAAGTGCTACAGTCAATAC‑AAG‑AAGAAAAAAAAA‑CACAAATTGTACGGTAACGTTTATATGTGCCATTTTTTAAGTACAATCACCATATATTCCTAAc                           >  1:3796244/1‑100 (MQ=255)
                                                                         aCAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAAcc                          <  2:2138849/100‑1 (MQ=255)
                                                                          cAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCa                         >  1:1033386/1‑100 (MQ=255)
                                                                          cAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCa                         >  1:1033229/1‑100 (MQ=255)
                                                                          cAAGTGCTACAGTCAATAA‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCa                         <  1:1023958/100‑1 (MQ=255)
                                                                          cAAGGGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCa                         <  2:1800579/100‑1 (MQ=255)
                                                                           aaGTGCTACAGTCAATAC‑AAG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCAt                        <  1:2232343/100‑1 (MQ=255)
                                                                            aGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATa                       >  2:2252966/1‑100 (MQ=255)
                                                                             gTGCTAAAGTCAATTC‑ATG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGGCATTTGTTCCGTACAAGCACTCCACATTCCTAACCATaa                      <  1:1651403/100‑1 (MQ=255)
                                                                                 tACAGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAAccc                  <  2:2619314/100‑1 (MQ=255)
                                                                                    aGTCAATAC‑ACG‑AGGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATATGTGCCATTTTTTTAGTACAATCACCCGATATTCCGAAACATAAACCCTaa               >  1:2885436/1‑100 (MQ=255)
                                                                                    aGTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTa                <  1:328775/99‑1 (MQ=255)
                                                                                     gTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTTCTGTAACGTTTATCAGTGCCAATTTTTCAGCACAAACCCCCAATATTCCTTATCCTAAAccc                  >  2:563722/1‑96 (MQ=255)
                                                                                     gTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTTTTTTTTTTTTCCCTATATTCCTAACCATAAACCCt                 >  1:408863/1‑97 (MQ=255)
                                                                                     gTCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAg              >  2:2715166/1‑100 (MQ=255)
                                                                                      tCAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGt             <  2:1343239/100‑1 (MQ=255)
                                                                                       cAATAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGtt            >  1:893432/1‑100 (MQ=255)
                                                                                        aaTAC‑ACG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTg           <  2:1507989/100‑1 (MQ=255)
                                                                                          tAC‑ACG‑AAGAAAAAAAAA‑CACAACTCGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGcc         <  2:1062892/100‑1 (MQ=255)
                                                                                           ac‑acG‑AAG‑AAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCtt       <  1:1681445/100‑1 (MQ=255)
                                                                                               cG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTg     <  2:446829/100‑1 (MQ=255)
                                                                                               cG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTg     <  2:3595711/100‑1 (MQ=255)
                                                                                               cG‑AAGAAAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTg     <  2:3228765/100‑1 (MQ=255)
                                                                                                  aag‑aaAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTGTTc  <  2:74776/100‑1 (MQ=255)
                                                                                                     |                                                                                                  
CTTTTCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAG‑AAAAAAAA‑CACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTGTTC  >  NC_000913/3333476‑3333671

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: ‑