Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,704,306 A→T 15.3% intergenic (‑243/+29) 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,3060AT15.3% 302.2 / 41.2 119intergenic (‑243/+29)rnc/lepBRNase III/leader peptidase (signal peptidase I)
Reads supporting (aligned to +/- strand):  ref base A (45/55);  new base T (8/10);  total (53/66)
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.89e-01

Only 100 of 122 total aligned reads displayed.

TGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCA  >  NC_000913/2704207‑2704405
                                                                                                   |                                                                                                   
tgccaACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATa                                                                                                     >  1:1306657‑M2/6‑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)
        ttAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCGTTATGGCGAc                                                                                             <  1:2428268/100‑1 (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)
             atctgggggAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGt                                                                                        <  2:2156647/96‑1 (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:1606236/100‑2 (MQ=38)
                 gggggAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCGTTATGGCGACAACGTGaat                                                                                    <  2:81749/100‑2 (MQ=38)
                 gggggAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGaa                                                                                     <  1:1881186/99‑1 (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)
                                acacAACGACGGGGGATGTCGTTAGCCAAGGGAGATTTATCTCGTAAATAATTCAAGTTGTCGGAGTAACGGCGCGAACGTGAACGAAGATgacgagcaa                                                                     >  1:2531299/1‑92 (MQ=255)
                                 cacaACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAAt                                                                    <  1:1703360/100‑1 (MQ=255)
                                  gcaACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAAt                                                                    <  1:1221523/98‑1 (MQ=255)
                                  acaACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATg                                                                   <  1:95719/100‑1 (MQ=255)
                                  acaACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATg                                                                   <  2:1685594/100‑1 (MQ=255)
                                     acgacgGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGAt                                                                >  1:865970/1‑100 (MQ=255)
                                       gacgGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATgc                                                              <  2:1618009/100‑1 (MQ=255)
                                         ctggggATGTCGTGAGCCACGGGAGATTTATATCATAAATAAATCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATgccg                                                            >  2:1751788/3‑100 (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)
                                              gATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATTGATGCCGCCAAt                                                       <  2:1446978/100‑1 (MQ=255)
                                                tttCGTTAGCCACGGGAGATTTATCTCATAAATAATGCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGc                                                     <  1:874755/98‑1 (MQ=255)
                                                tGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGc                                                     >  2:1289591/1‑100 (MQ=255)
                                                tGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGc                                                     <  1:1775644/100‑1 (MQ=255)
                                                 gTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCg                                                    >  1:2157171/1‑100 (MQ=255)
                                                   cGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGAc                                                  >  1:201498/1‑100 (MQ=255)
                                                      tAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTa                                               <  2:2240961/100‑1 (MQ=255)
                                                        gCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAg                                             >  2:2070594/1‑100 (MQ=255)
                                                        gCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAg                                             >  1:1223486/1‑100 (MQ=255)
                                                          cACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAgcg                                           <  2:2084060/100‑1 (MQ=255)
                                                            cGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAgcgc                                          <  2:1411393/99‑1 (MQ=255)
                                                            cGGGAGATTTATCTCATAAATAATACACGATCTCGCCATAACGCCGACACCGTGAACGAAGATGCCTATTAGTGGCTACCGCCAATGCGAccaaag                                             >  1:2109803/1‑91 (MQ=255)
                                                             gggAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAg                                        <  2:1978548/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                             >  2:1667792/1‑100 (MQ=255)
                                                                        ctcATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGcc                             <  1:94670/100‑1 (MQ=255)
                                                                          cATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCAt                           >  1:224902/1‑100 (MQ=255)
                                                                          cATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCAt                           >  1:5465/1‑100 (MQ=255)
                                                                            taaataaTTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTc                         >  2:275276/1‑100 (MQ=255)
                                                                            taaataaTTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTc                         >  1:518070/1‑100 (MQ=255)
                                                                                taatTCACGTTGTCGCCAGAACGCCGACAACGTGAACGAAGATGGCTATTATTGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCct                     >  1:1025115/1‑100 (MQ=255)
                                                                                 aatTCACGTTGTCGCCATAGCGGCGACAACGTGGACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCctt                    <  1:186532/100‑1 (MQ=255)
                                                                                 aatTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCctt                    >  2:1846066/1‑100 (MQ=255)
                                                                                 aatTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCAGCCAATGCGACTTAAGCGCAGACAATTGGGCCATTCGCct                     >  1:1165969/1‑99 (MQ=255)
                                                                                   ttCACGTTGTCGCCGTTATGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACCTAAGCGCAGACCAGTCGGCCATTCGCcttct                  >  1:2252559/1‑100 (MQ=38)
                                                                                    tCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCcttctt                 >  2:4737/1‑100 (MQ=255)
                                                                                    tCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCcttctt                 <  1:1311316/100‑1 (MQ=255)
                                                                                     cACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTg                >  1:197752/1‑100 (MQ=255)
                                                                                      aCGTTGTCGCCGTTATGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGc               <  1:2158411/100‑1 (MQ=38)
                                                                                      aCGTTGTCGCCGTTATGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCCTCTTGc               >  2:1280917/1‑100 (MQ=38)
                                                                                      aCGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGc               >  1:1292600/1‑100 (MQ=255)
                                                                                      aCGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGc               <  2:197754/100‑1 (MQ=255)
                                                                                      aCGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGc               <  2:1632373/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:1433035/1‑100 (MQ=37)
                                                                                        gTTGTCGCCGTTATGGCGACAACGTGAATGAAGATGGCTATTAATGGATGCCGCCAACGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCtt             >  1:1433151/1‑100 (MQ=38)
                                                                                        gTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCtt             >  2:839659/1‑100 (MQ=255)
                                                                                        gTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCCCCTTCTTGCtt             >  1:1264325/1‑100 (MQ=255)
                                                                                             cGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGa        >  1:286197/1‑100 (MQ=255)
                                                                                             cGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGa        >  2:1738583/1‑100 (MQ=255)
                                                                                              gCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATCCGACTTAAGCGCAGACCAGTCGGCCAATCGCCTTCTTGCTTATCg         >  1:345155/1‑98 (MQ=255)
                                                                                              gCCATAACGGCGACAACGTGAACGAAAATGGCTATTAATGGACGCCGCCAATGCCACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGaa       >  1:962102/1‑100 (MQ=255)
                                                                                                cATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGc     <  1:2247668/100‑1 (MQ=255)
                                                                                                cATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGc     <  2:60558/100‑1 (MQ=255)
                                                                                                cATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGACGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGc     <  1:1812624/100‑1 (MQ=255)
                                                                                                 aTAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCt    <  2:2579987/100‑1 (MQ=255)
                                                                                                 aTAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCt    <  1:1549869/100‑1 (MQ=255)
                                                                                                 aTAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCt    <  2:1880591/100‑1 (MQ=255)
                                                                                                  tAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTc   >  2:2275181/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
TGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCA  >  NC_000913/2704207‑2704405

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: