Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,704,308 C→T 16.0% intergenic (‑245/+27) 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,3080CT16.0% 287.9 / 39.0 119intergenic (‑245/+27)rnc/lepBRNase III/leader peptidase (signal peptidase I)
Reads supporting (aligned to +/- strand):  ref base C (46/53);  new base T (8/11);  total (54/65)
Fisher's exact test for biased strand distribution p-value = 8.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 120 total aligned reads displayed.

CCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCA  >  NC_000913/2704209‑2704405
                                                                                                   |                                                                                                 
ccaACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAAc                                                                                                   <  1:1009244‑M2/97‑1 (MQ=255)
 caACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACg                                                                                                  <  2:1478914‑M2/98‑1 (MQ=255)
  aACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACgg                                                                                                 <  1:1311192‑M2/99‑1 (MQ=255)
     cTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGa                                                                                              <  2:763262/100‑1 (MQ=255)
      ttAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCGTTATGGCGAc                                                                                             <  1:2428268/100‑1 (MQ=255)
      ttAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGAc                                                                                             <  2:884269/100‑1 (MQ=255)
      ttAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGAc                                                                                             >  1:2084046/1‑100 (MQ=255)
          aCTTCGGGGGAATATTCTATACTCAACGACGGGGTATCTCGTTAGCCGCGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACg                                                                                         <  1:1747921/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:81749/100‑2 (MQ=38)
               gggggAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGtagat                                                                                    >  1:2156633/1‑96 (MQ=255)
               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)
                          cTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGc                                                                           <  1:1781832/98‑1 (MQ=255)
                                gcaACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAAt                                                                    <  1:1221523/98‑1 (MQ=255)
                                acaACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATg                                                                   <  1:2292553/100‑1 (MQ=255)
                                acaACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATg                                                                   <  1:95719/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:401939/100‑1 (MQ=255)
                                       cgGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATgccg                                                            >  2:1401273/1‑100 (MQ=255)
                                         ggggATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATgccgcc                                                          <  1:1017743/100‑1 (MQ=255)
                                         ggggATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATgccgcc                                                          <  1:1734501/100‑1 (MQ=255)
                                          gggATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCa                                                         >  2:27390/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)
                                                 cGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGAc                                                  <  2:375502/100‑1 (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: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                             >  2:1668241/1‑100 (MQ=255)
                                                                      ctcATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGcc                             <  1:1913850/100‑1 (MQ=255)
                                                                      ctcATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGcc                             <  1:94670/100‑1 (MQ=255)
                                                                        cATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCAt                           >  1:5465/1‑100 (MQ=255)
                                                                          taaataaTTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTc                         >  1:518070/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                 >  2:4737/1‑100 (MQ=255)
                                                                                   cACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCcttctt                 <  2:2435412/99‑1 (MQ=255)
                                                                                   cACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTg                >  1:197752/1‑100 (MQ=255)
                                                                                   cACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTg                >  2:1537320/1‑100 (MQ=255)
                                                                                    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:387525/1‑100 (MQ=255)
                                                                                           cGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGa        >  2:1738583/1‑100 (MQ=255)
                                                                                           cGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGa        >  1:286197/1‑100 (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)
                                                                                                tAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTc   >  2:2275181/1‑100 (MQ=255)
                                                                                                 aaCGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCCCCTTCTTGCTTATCGAAGCTCa  >  1:2194380/1‑100 (MQ=255)
                                                                                                   |                                                                                                 
CCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCA  >  NC_000913/2704209‑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: