Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,704,304 A→G 17.0% intergenic (‑241/+31) rnc ← / ← lepB RNase III/leader peptidase (signal peptidase I)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,704,3040AG17.0% 263.8 / 43.5 117intergenic (‑241/+31)rnc/lepBRNase III/leader peptidase (signal peptidase I)
Reads supporting (aligned to +/- strand):  ref base A (45/52);  new base G (9/11);  total (54/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

Only 100 of 121 total aligned reads displayed.

GGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCT  >  NC_000913/2704205‑2704403
                                                                                                   |                                                                                                   
ggTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCa                                                                                                     >  2:2540806/1‑100 (MQ=255)
  tgccaACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATa                                                                                                   >  1:1306657‑M2/6‑100 (MQ=255)
   gccaACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATCAATAATTCACGTTGTCGCCATaa                                                                                                  >  2:1772069‑M2/5‑100 (MQ=255)
    ccaACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAAc                                                                                                 <  1:1009244‑M2/97‑1 (MQ=255)
     caACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACg                                                                                                <  2:1478914‑M2/98‑1 (MQ=255)
        ccTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCg                                                                                             <  2:2277924/100‑1 (MQ=255)
         cTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGa                                                                                            >  1:1864878/1‑100 (MQ=255)
          ttAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGAc                                                                                           <  2:884269/100‑1 (MQ=255)
          ttAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGAc                                                                                           >  1:2084046/1‑100 (MQ=255)
             aaCTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAAc                                                                                        <  2:710876/100‑1 (MQ=255)
              aCTTCGGGGGAATATTCTATACTCAACGACGGGGTATCTCGTTAGCCGCGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACg                                                                                       <  1:1747921/100‑1 (MQ=255)
              aCTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCGTTATGGCGACAACg                                                                                       <  2:2149098/100‑1 (MQ=255)
              aCTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCGTTATGGCGACAACg                                                                                       <  2:1611293/100‑1 (MQ=38)
              aCTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCGTTATGGCGACAACg                                                                                       <  1:700133/100‑1 (MQ=38)
              aCTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACg                                                                                       <  2:1013610/100‑1 (MQ=255)
               cTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACg                                                                                       <  1:2311932/99‑1 (MQ=255)
               cTTCGGGGGAATATTCTATACACAACGACACGGGATGTCGTTAGACACGGGAGATTTATCTCATAAAGAATTCACAATATCGCCATAACGGCGACAACGt                                                                                      >  1:273317/1‑100 (MQ=255)
                ttCGGGGGAATATTCTATACACAACGACGGGGGATTTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGt                                                                                      <  1:2225816/99‑1 (MQ=255)
                ttCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATTATTCACGTTGTCGCCGTGATGGCGACAAcct                                                                                      <  2:1747044/99‑3 (MQ=39)
                ttCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCGTTATGGCGACAACGt                                                                                      <  1:542060/99‑1 (MQ=38)
                ttCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCGTTATGGCGACAACGTg                                                                                     <  2:2489563/100‑1 (MQ=38)
                   gggggAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCGTTATGGCGACAACGTGaat                                                                                  <  2:81749/100‑2 (MQ=38)
                   gggggAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCGTTATGGCGACAACGTGaat                                                                                  <  2:1606236/100‑2 (MQ=38)
                   gggggAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGtagat                                                                                  >  1:2156633/1‑96 (MQ=255)
                   gggggAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGaa                                                                                   <  1:1881186/99‑1 (MQ=255)
                   gggggAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAAc                                                                                  >  1:842470/1‑100 (MQ=255)
                   gggggAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAAc                                                                                  <  2:1232788/100‑1 (MQ=255)
                    ggggAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACg                                                                                 <  2:1865241/100‑1 (MQ=255)
                       gAATATTCTATACACAACCACTAGCCATGCCGTCACCCACGGCAGATTTACCTCATAAATAATTCACGTTGTCGCCGTTATGGCGACAACGTg                                                                                     <  1:2158078/93‑1 (MQ=12)
                        aaTATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACg                                                                                 <  2:1280693/96‑1 (MQ=255)
                          tatTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATg                                                                           <  1:1041894/100‑1 (MQ=255)
                            ttCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGc                                                                         >  1:2573901/1‑100 (MQ=255)
                            ttCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGc                                                                         >  1:2310983/1‑100 (MQ=255)
                              cTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGc                                                                         <  1:1781832/98‑1 (MQ=255)
                                  acacAACGACGGGGGATGTCGTTAGCCAAGGGAGATTTATCTCGTAAATAATTCAAGTTGTCGGAGTAACGGCGCGAACGTGAACGAAGATgacgagcaa                                                                   >  1:2531299/1‑92 (MQ=255)
                                   cacaACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAAt                                                                  <  1:1703360/100‑1 (MQ=255)
                                    gcaACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAAt                                                                  <  1:1221523/98‑1 (MQ=255)
                                    acaACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATg                                                                 <  2:1685594/100‑1 (MQ=255)
                                    acaACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATg                                                                 <  1:2292553/100‑1 (MQ=255)
                                       acgacgGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGAt                                                              >  1:865970/1‑100 (MQ=255)
                                         gacgGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATgc                                                            <  2:1618009/100‑1 (MQ=255)
                                           cgGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATgccg                                                          >  2:1401273/1‑100 (MQ=255)
                                           cgGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATgccg                                                          <  2:401939/100‑1 (MQ=255)
                                             ggggATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATgccgcc                                                        <  1:1017743/100‑1 (MQ=255)
                                             ggggATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATgccgcc                                                        <  1:1734501/100‑1 (MQ=255)
                                              gggATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCa                                                       >  2:27390/1‑100 (MQ=255)
                                               ggATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCaa                                                      >  1:2078774/1‑100 (MQ=255)
                                                  tttCGTTAGCCACGGGAGATTTATCTCATAAATAATGCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGc                                                   <  1:874755/98‑1 (MQ=255)
                                                  tGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGc                                                   >  2:1289591/1‑100 (MQ=255)
                                                  tGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGc                                                   <  1:1775644/100‑1 (MQ=255)
                                                     cGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGAc                                                >  1:201498/1‑100 (MQ=255)
                                                        tAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTa                                             <  2:2240961/100‑1 (MQ=255)
                                                          gCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAg                                           >  1:1223486/1‑100 (MQ=255)
                                                          gCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAg                                           >  2:2070594/1‑100 (MQ=255)
                                                            cACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAgcg                                         <  2:2084060/100‑1 (MQ=255)
                                                              cGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAgcgc                                        <  2:1411393/99‑1 (MQ=255)
                                                              cGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCa                                       <  1:2541075/100‑1 (MQ=255)
                                                              cGGGAGATTTATCTCATAAATAATACACGATCTCGCCATAACGCCGACACCGTGAACGAAGATGCCTATTAGTGGCTACCGCCAATGCGAccaaag                                           >  1:2109803/1‑91 (MQ=255)
                                                               gggAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAg                                      <  2:1361521/100‑1 (MQ=255)
                                                               gggAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAg                                      <  2:812141/100‑1 (MQ=255)
                                                                ggAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGa                                     <  1:1451746/100‑1 (MQ=255)
                                                                 gagaTTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGAc                                    <  2:2379329/100‑1 (MQ=255)
                                                                       tATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCg                              >  2:545285/1‑100 (MQ=255)
                                                                          ctcATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGcc                           <  1:1913850/100‑1 (MQ=255)
                                                                          ctcATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGcc                           >  2:1668241/1‑100 (MQ=255)
                                                                          ctcATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGcc                           >  2:1667792/1‑100 (MQ=255)
                                                                          ctcATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGcc                           <  1:94670/100‑1 (MQ=255)
                                                                            cATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCAt                         >  1:224902/1‑100 (MQ=255)
                                                                              taaataaTTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTc                       >  1:518070/1‑100 (MQ=255)
                                                                              taaataaTTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTc                       >  2:275276/1‑100 (MQ=255)
                                                                               aaataaTTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCg                      >  1:1545160/1‑100 (MQ=255)
                                                                                cttcgttCACGTTGTCGCCGTTATGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGc                     >  2:2158425/6‑100 (MQ=37)
                                                                                  taatTCACGTTGTCGCCAGAACGCCGACAACGTGAACGAAGATGGCTATTATTGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCct                   >  1:1025115/1‑100 (MQ=255)
                                                                                   aatTCACGTTGTCGCCATAGCGGCGACAACGTGGACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCctt                  <  1:186532/100‑1 (MQ=255)
                                                                                   aatTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCctt                  >  2:1846066/1‑100 (MQ=255)
                                                                                   aatTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCctt                  >  2:1993008/1‑100 (MQ=255)
                                                                                   aatTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCAGCCAATGCGACTTAAGCGCAGACAATTGGGCCATTCGCct                   >  1:1165969/1‑99 (MQ=255)
                                                                                     ttCACGTTGTCGCCGTTATGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACCTAAGCGCAGACCAGTCGGCCATTCGCcttct                >  1:2252559/1‑100 (MQ=38)
                                                                                      tCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCcttctt               <  1:1311316/100‑1 (MQ=255)
                                                                                      tCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCcttctt               >  2:4737/1‑100 (MQ=255)
                                                                                       cACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTg              >  2:1537320/1‑100 (MQ=255)
                                                                                        aCGTTGTCGCCGTTATGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGc             <  1:2158411/100‑1 (MQ=38)
                                                                                        aCGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGc             >  1:1292600/1‑100 (MQ=255)
                                                                                        aCGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGc             <  2:197754/100‑1 (MQ=255)
                                                                                         cGTTGTCGCCGTTATGGCGACGACGTGAATGAAGATGGCAGTTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCAATCGCCTTCTTGCt            >  1:345088/1‑100 (MQ=37)
                                                                                         cGTTGTCGCCGTTATGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCt            >  2:2192493/1‑100 (MQ=38)
                                                                                         cGTTGTCGCCGTTATGGCGACAAAGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCt            >  1:2189229/1‑100 (MQ=38)
                                                                                         cGTTGTCGCCATAAAGGCGACAACGTGAACGAAGTTGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTg              <  1:2206353/98‑1 (MQ=255)
                                                                                          gTTGTCGCCGTTATGGCGACAACGTGAATGAAGATGGCTATTAATGGATGCCGCCAACGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCtt           >  1:1433151/1‑100 (MQ=38)
                                                                                          gTTGTCGCCGTTATGGCGACAACGTGAATGAAGATGGCTATTAATGGATGCCGCCAACGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCtt           >  1:1433035/1‑100 (MQ=37)
                                                                                          gTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCtt           >  2:839659/1‑100 (MQ=255)
                                                                                          gTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCCCCTTCTTGCtt           >  1:1264325/1‑100 (MQ=255)
                                                                                               cGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGa      >  1:387525/1‑100 (MQ=255)
                                                                                               cGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGa      >  2:1738583/1‑100 (MQ=255)
                                                                                                gCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGaa     >  1:971819/1‑100 (MQ=255)
                                                                                                gCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATCCGACTTAAGCGCAGACCAGTCGGCCAATCGCCTTCTTGCTTATCg       >  1:345155/1‑98 (MQ=255)
                                                                                                gCCATAACGGCGACAACGTGAACGAAAATGGCTATTAATGGACGCCGCCAATGCCACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGaa     >  1:962102/1‑100 (MQ=255)
                                                                                                  cATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGc   <  2:60558/100‑1 (MQ=255)
                                                                                                   aTAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCt  <  2:2579987/100‑1 (MQ=255)
                                                                                                   aTAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCt  <  2:1880591/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCT  >  NC_000913/2704205‑2704403

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