Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,866,940 C→A 5.3% K204N (AAG→AAT) ‡ rpoS ← RNA polymerase, sigma S (sigma 38) factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,866,9400CA 5.3% 551.8 / 16.2 151K204N (AAG→AAT) ‡rpoSRNA polymerase, sigma S (sigma 38) factor
Reads supporting (aligned to +/- strand):  ref base C (74/69);  new base A (4/4);  total (78/73)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.71e-01

Only 100 of 153 total aligned reads displayed.

TCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGC  >  NC_000913/2866794‑2867083
                                                                                                                                                  |                                                                                                                                               
      tGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCtt                                                                                                                                               <  4:575178/143‑1 (MQ=255)
      tGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCtt                                                                                                                                               >  3:575178/1‑143 (MQ=255)
      tGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCtt                                                                                                                                               <  3:226200/143‑1 (MQ=255)
                         cGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTc                                                                                                                       <  2:285112/148‑1 (MQ=255)
                          gTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTAt                                                                                                                                             >  4:766777/1‑125 (MQ=255)
                          gTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTAt                                                                                                                                             <  3:766777/125‑1 (MQ=255)
                           ttCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCtt                                                                                                                        >  2:232942/1‑145 (MQ=255)
                           ttCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCtt                                                                                                                        <  1:232942/145‑1 (MQ=255)
                            tctcTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGc                                                                                                                    >  2:881079/1‑148 (MQ=255)
                               ctTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACt                                                                                                                 <  1:881079/148‑1 (MQ=255)
                                 ttttCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCt                                                                                                                                                <  4:892053/115‑1 (MQ=255)
                                 ttttCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCt                                                                                                                                                >  3:892053/1‑115 (MQ=255)
                                 ttttCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATctct                                                                                                                         <  1:164782/138‑1 (MQ=255)
                                           ccAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTg                                                                                                               >  3:1091583/1‑138 (MQ=255)
                                           ccAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTg                                                                                                               <  4:1091583/138‑1 (MQ=255)
                                            cAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATctct                                                                                                                         <  2:840174/127‑1 (MQ=255)
                                            cAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATctct                                                                                                                         >  1:840174/1‑127 (MQ=255)
                                            cAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTcc                                                                                                    <  1:333683/148‑1 (MQ=255)
                                             aGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTgg                                                                                                              >  1:1027216/1‑137 (MQ=255)
                                             aGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTgg                                                                                                              <  2:1027216/137‑1 (MQ=255)
                                              ggATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTg                                                                                                                                 >  3:871181/1‑117 (MQ=255)
                                                aTGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTc                                                                                                                       <  4:397273/125‑1 (MQ=255)
                                                aTGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTc                                                                                                                       >  3:397273/1‑125 (MQ=255)
                                                aTGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCt                                                                                                >  3:102835/1‑148 (MQ=255)
                                                aTGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCt                                                                                                >  2:849873/1‑148 (MQ=255)
                                                   tCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCAGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGATCATGGt                                                                                                      >  2:644252/1‑139 (MQ=255)
                                                   tCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGAATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGt                                                                                                      <  1:644252/139‑1 (MQ=255)
                                                         aaCGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATg                                                                                            >  2:401249/1‑143 (MQ=255)
                                                          aCGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTgg                                                                                                              >  3:69151/1‑124 (MQ=255)
                                                          aCGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTgg                                                                                                              <  4:69151/124‑1 (MQ=255)
                                                          aCGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATg                                                                                            <  3:254424/142‑1 (MQ=255)
                                                           cGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTAt                                                                                             >  2:665649/1‑140 (MQ=255)
                                                           cGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTAt                                                                                             <  1:665649/140‑1 (MQ=255)
                                                             cTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGATAAGCCAGTTga                                                                                                                                     <  1:405320/98‑2 (MQ=255)
                                                              tttttCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCa                                                                                  <  4:102835/148‑1 (MQ=255)
                                                               ttttCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATggg                                                                                          >  1:204328/1‑139 (MQ=255)
                                                               ttttCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATggg                                                                                          <  2:204328/139‑1 (MQ=255)
                                                               ttttCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATggg                                                                                          <  2:347206/139‑1 (MQ=255)
                                                               ttttCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCAc                                                                                 >  3:1027862/1‑148 (MQ=255)
                                                                   cGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGc                                                                             <  4:1027862/148‑1 (MQ=255)
                                                                          accaccCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACg                                                                                <  1:1035875/138‑1 (MQ=255)
                                                                          accaccCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACg                                                                                >  2:1035875/1‑138 (MQ=255)
                                                                            caccCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCAt                                                                                                         <  2:233414/111‑1 (MQ=255)
                                                                            caccCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCAt                                                                                                         >  1:233414/1‑111 (MQ=255)
                                                                              cccAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTAt                                                                                                                                             <  1:1007810/73‑1 (MQ=255)
                                                                              cccAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTAt                                                                                                                                             >  2:1007810/1‑73 (MQ=255)
                                                                              cccAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGATAAGCCAGTTgat                                                                                                                                    <  1:130070/82‑3 (MQ=255)
                                                                                       gtgtCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCa                                                         <  4:230861/148‑1 (MQ=255)
                                                                                        tgtCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGt                                                                          <  4:751741/130‑1 (MQ=255)
                                                                                        tgtCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTaaa                                                               <  1:480341/141‑1 (MQ=255)
                                                                                        tgtCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTaaa                                                               >  2:480341/1‑141 (MQ=255)
                                                                                             aCCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCtt                                                                                                                        <  4:181020/79‑1 (MQ=255)
                                                                                             aCCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCtt                                                                                                                        >  3:181020/1‑79 (MQ=255)
                                                                                             aCCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCt                                                                                                <  2:170038/103‑1 (MQ=255)
                                                                                             aCCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCt                                                                                                >  1:170038/1‑103 (MQ=255)
                                                                                             aCCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCt                                                   >  3:24712/1‑148 (MQ=255)
                                                                                                      aTGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTc                                                                                                                       >  1:378537/1‑71 (MQ=255)
                                                                                                      aTGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTc                                                                                                                       <  2:378537/71‑1 (MQ=255)
                                                                                                           cTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCAt                                                                                                         <  4:164987/80‑1 (MQ=255)
                                                                                                           cTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCAt                                                                                                         >  3:164987/1‑80 (MQ=255)
                                                                                                               ttAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGt                                                                                                             >  3:526317/1‑72 (MQ=255)
                                                                                                               ttAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGt                                                                                                             <  4:526317/72‑1 (MQ=255)
                                                                                                               ttAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGtt                                                           <  1:16390/122‑1 (MQ=255)
                                                                                                               ttAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGtt                                                           >  2:16390/1‑122 (MQ=255)
                                                                                                               ttAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCa                                 >  2:180761/1‑148 (MQ=255)
                                                                                                                tAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTg                                                                                                               >  2:159664/1‑69 (MQ=255)
                                                                                                                 aaGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCaaa                               <  4:24712/148‑1 (MQ=255)
                                                                                                                   gACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTa                                                                                              <  1:1052756/83‑1 (MQ=255)
                                                                                                                   gACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTa                                                                                              >  2:1052756/1‑83 (MQ=255)
                                                                                                                   gACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGt                                                            <  4:414896/117‑1 (MQ=255)
                                                                                                                   gACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGt                                                            >  4:565196/1‑117 (MQ=255)
                                                                                                                   gACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGt                                                            >  3:414896/1‑117 (MQ=255)
                                                                                                                   gACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGt                                                            <  3:565196/117‑1 (MQ=255)
                                                                                                                    aCGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGa                            <  1:180761/148‑1 (MQ=255)
                                                                                                                       aaGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTc                                                                                                                       <  1:803466/54‑1 (MQ=255)
                                                                                                                       aaGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTc                                                                                                                       >  2:803466/1‑54 (MQ=255)
                                                                                                                         gCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGt                       >  1:953949/1‑148 (MQ=255)
                                                                                                                          cATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCaaa                               <  4:261381/139‑1 (MQ=255)
                                                                                                                          cATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCaaa                               >  3:261381/1‑139 (MQ=255)
                                                                                                                               ggCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTg                                                                              <  2:322179/87‑1 (MQ=255)
                                                                                                                               ggCTAACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTg                                                                              >  1:322179/1‑87 (MQ=255)
                                                                                                                                gCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAg                        <  3:679449/140‑1 (MQ=255)
                                                                                                                                  tGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTg              <  2:580303/148‑1 (MQ=255)
                                                                                                                                   gACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTgg             <  4:258425/148‑1 (MQ=255)
                                                                                                                                     cGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGtt                                                           <  1:169364/100‑1 (MQ=255)
                                                                                                                                     cGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGtt                                                           >  2:169364/1‑100 (MQ=255)
                                                                                                                                       tcatcaACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCa                                                                     <  3:617468/88‑1 (MQ=255)
                                                                                                                                       tcatcaACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCa                                                                     >  4:617468/1‑88 (MQ=255)
                                                                                                                                       tcatcaACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTc          >  4:794056/1‑147 (MQ=255)
                                                                                                                                        catcaACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGc                                                                      >  1:60629/1‑86 (MQ=255)
                                                                                                                                          tcaACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACg                                                                                <  1:621236/74‑1 (MQ=255)
                                                                                                                                          tcaACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACg                                                                                >  2:621236/1‑74 (MQ=255)
                                                                                                                                          tcaACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTg              <  2:821065/140‑1 (MQ=255)
                                                                                                                                          tcaACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTg              >  1:821065/1‑140 (MQ=255)
                                                                                                                                           caACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCa                                                                     <  1:1078931/84‑1 (MQ=255)
                                                                                                                                           caACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCa                                                                     >  2:1078931/1‑84 (MQ=255)
                                                                                                                                           caACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAAt                                      >  2:696554/1‑115 (MQ=255)
                                                                                                                                           caACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAAt                                      <  1:696554/115‑1 (MQ=255)
                                                                                                                                             aCTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCg   >  2:395715/1‑148 (MQ=255)
                                                                                                                                                  cTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTc                                                          >  3:765162/1‑88 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
TCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGC  >  NC_000913/2866794‑2867083

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