Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,843,965 G→T 7.0% R303S (CGC→AGC)  ynjI ← inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,843,9650GT 7.0% 340.2 / 9.9 98R303S (CGC→AGC) ynjIinner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (45/44);  new base T (3/4);  total (49/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGAAGCAAGCATTATATAATGCCAGGCCAGTTCTTCTTCAATTGTCCCGTT  >  NC_000913/1843820‑1844103
                                                                                                                                                 |                                                                                                                                          
aTTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGAt                                                                                                                                          <  3:497467/148‑1 (MQ=255)
                    tttAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATg                                                                                                                             <  3:119689/141‑1 (MQ=255)
                    tttAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATg                                                                                                                             >  4:119689/1‑141 (MQ=255)
                                      ccACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATc                                                                                                               >  2:520910/1‑137 (MQ=255)
                                      ccACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATc                                                                                                               <  1:520910/137‑1 (MQ=255)
                                       cACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCAATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTAt                                                                                                   >  4:1028632/1‑148 (MQ=255)
                                        aCGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATg                                                                                                                             <  1:1069914/121‑1 (MQ=255)
                                        aCGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATg                                                                                                                             >  2:1069914/1‑121 (MQ=255)
                                             aCATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGata                                                                                                                               <  4:918303/114‑1 (MQ=255)
                                             aCATTTCTTAGGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGata                                                                                                                               >  3:918303/1‑114 (MQ=255)
                                                tttCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACGCGCTTCAAGTTCAGCATCCGTTAACTTTCTGCTATcgcagaa                                                                                                                                   >  3:953051/1‑105 (MQ=255)
                                                tttCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCTATcgcagaa                                                                                                                                   <  4:953051/107‑3 (MQ=255)
                                                 ttCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATa                                                                                         >  3:501790/1‑148 (MQ=255)
                                                         ggCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCTATcgcagaa                                                                                                                                   <  3:512639/98‑3 (MQ=255)
                                                         ggCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCTATcgcagaa                                                                                                                                   <  1:293910/98‑3 (MQ=255)
                                                         ggCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCTATcgcagaa                                                                                                                                   >  4:512639/1‑96 (MQ=255)
                                                         ggCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCTATcgcagaa                                                                                                                                   >  2:293910/1‑96 (MQ=255)
                                                          gCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCttcttc                                                                                >  2:698436/1‑148 (MQ=255)
                                                            cAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGaaa                                                                                                                 <  1:383261/113‑1 (MQ=255)
                                                            cAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGaaa                                                                                                                 >  2:383261/1‑113 (MQ=255)
                                                             aaCGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCAt                                                                                                           >  3:226069/1‑118 (MQ=255)
                                                             aaCGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCAt                                                                                                           <  4:226069/118‑1 (MQ=255)
                                                              aCGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCtt                                                                            >  3:313703/1‑148 (MQ=255)
                                                              aCGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCtt                                                                            <  4:313703/148‑1 (MQ=255)
                                                                     ttCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAAt                                                                                       >  2:931790/1‑130 (MQ=255)
                                                                     ttCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAAt                                                                                       <  1:931790/130‑1 (MQ=255)
                                                                     ttCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGaaca                                                                     <  4:501790/148‑1 (MQ=255)
                                                                     ttCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCAATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGaaca                                                                     <  3:1028632/148‑1 (MQ=255)
                                                                         tACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGAt                                                                                          >  3:1038052/1‑123 (MQ=255)
                                                                         tACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGAt                                                                                          <  4:1038052/123‑1 (MQ=255)
                                                                            tGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTg                                                              <  1:698436/148‑1 (MQ=255)
                                                                             gTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTgc                                                             >  2:321816/1‑148 (MQ=255)
                                                                              tAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTAttttt                                                                                               <  3:345008/113‑1 (MQ=255)
                                                                              tAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTAttttt                                                                                               >  4:345008/1‑113 (MQ=255)
                                                                                 tGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTg                                                                                            <  4:332181/113‑1 (MQ=255)
                                                                                 tGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTg                                                                                            >  3:332181/1‑113 (MQ=255)
                                                                                     cccAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCTATagc                                                                                                                                       <  4:268266/66‑1 (MQ=255)
                                                                                     cccAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATagc                                                                                                                                       >  3:268266/1‑66 (MQ=255)
                                                                                     cccAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAAt                                                                                                                <  4:1047657/89‑1 (MQ=255)
                                                                                     cccAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAAt                                                                                                                >  3:1047657/1‑89 (MQ=255)
                                                                                     cccAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAAtcttct                                                                                  <  1:427714/119‑1 (MQ=255)
                                                                                     cccAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAAtcttct                                                                                  >  2:427714/1‑119 (MQ=255)
                                                                                     cccAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTg                                                              >  3:139309/1‑139 (MQ=255)
                                                                                     cccAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTg                                                              <  4:139309/139‑1 (MQ=255)
                                                                                       cAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGt                                                                                                                        <  2:577699/79‑1 (MQ=255)
                                                                                       cAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGt                                                                                                                        >  1:577699/1‑79 (MQ=255)
                                                                                       cAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGaa                                                                                                                  >  3:481390/1‑85 (MQ=255)
                                                                                       cAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGaa                                                                                                                  <  4:481390/85‑1 (MQ=255)
                                                                                       cAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCgag                                                   >  1:266041/1‑148 (MQ=255)
                                                                                         aaaCTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCt                                                                             <  2:886861/120‑1 (MQ=255)
                                                                                         aaaCTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCt                                                                             >  1:886861/1‑120 (MQ=255)
                                                                                            cTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGaaa                                                                                                                 >  4:868081/1‑81 (MQ=255)
                                                                                            cTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGaaa                                                                                                                 <  3:868081/81‑1 (MQ=255)
                                                                                             tGCTTGTCGTCACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATa                                                                                         <  1:193186/104‑1 (MQ=255)
                                                                                             tGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATa                                                                                         >  2:193186/1‑104 (MQ=255)
                                                                                              gCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTAtttttt                                                                                              >  2:495346/1‑98 (MQ=255)
                                                                                              gCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTAtttttt                                                                                              <  1:495346/98‑1 (MQ=255)
                                                                                                       aaCGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGAAGCAAGCATtatata                                   >  3:300589/1‑148 (MQ=255)
                                                                                                           aTAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGAagcaagca                                          >  3:62824/1‑137 (MQ=255)
                                                                                                           aTAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGAagcaagca                                          <  4:62824/137‑1 (MQ=255)
                                                                                                           aTAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGAAGCAAGCATTATATAATgc                               >  4:434057/1‑148 (MQ=255)
                                                                                                            tAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGAAGCAAGCATTATATAATgcc                              <  3:434057/148‑1 (MQ=255)
                                                                                                                 cGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGaa                                                                       <  1:133106/102‑1 (MQ=255)
                                                                                                                 cGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGaa                                                                       >  2:133106/1‑102 (MQ=255)
                                                                                                                 cGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGAagcaagca                                          <  3:161834/131‑1 (MQ=255)
                                                                                                                 cGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGAagcaagca                                          >  4:161834/1‑131 (MQ=255)
                                                                                                                    ttCAAGTTCAGCATCCGTTAACTTTCTGCGACAGCAGCAGATATGCCAGTAAAGATATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTTGAGAAGCAAGCATTATATAATGCCAGGCCAGt                      <  1:321816/148‑1 (MQ=255)
                                                                                                                     tCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAAtcttctg                                                                                 >  4:784271/1‑87 (MQ=255)
                                                                                                                     tCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATcttctt                                                                                 <  3:784271/88‑1 (MQ=255)
                                                                                                                     tCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTg                                                              <  3:511778/107‑1 (MQ=255)
                                                                                                                     tCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTg                                                              >  4:511778/1‑107 (MQ=255)
                                                                                                                         gTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGAAGCAAGCATTATATAATGCCAGGCCAGttcttc                 <  2:266041/148‑1 (MQ=255)
                                                                                                                          ttCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAAtcttct                                                                                  >  4:1064495/1‑82 (MQ=255)
                                                                                                                          ttCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAAtcttct                                                                                  <  3:1064495/82‑1 (MQ=255)
                                                                                                                           tCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTAttttt                                                                                               >  2:522377/1‑68 (MQ=255)
                                                                                                                           tCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTAttttt                                                                                               <  1:522377/68‑1 (MQ=255)
                                                                                                                           tCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGAAGCAAGCATTATATAATgcca                             >  4:13515/1‑134 (MQ=255)
                                                                                                                           tCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGAAGCAAGCATTATATAATgcca                             <  3:13515/134‑1 (MQ=255)
                                                                                                                               cATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGAAGCAAGCATTATATAATGCCAGGCCAGTTCTTCTTCAAt           >  4:942879/1‑148 (MQ=255)
                                                                                                                                aTCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGAAGCAAGCATTATATAATgccaggcc                         >  3:288444/1‑133 (MQ=255)
                                                                                                                                aTCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGAAGCAAGCATTATATAATgccaggcc                         <  4:288444/133‑1 (MQ=255)
                                                                                                                                 tCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGAAGCAAGCATTATATAATGCCAGGCCAGtt                     >  2:922058/1‑136 (MQ=255)
                                                                                                                                 tCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGAAGCAAGCATTATATAATGCCAGGCCAGtt                     <  1:922058/136‑1 (MQ=255)
                                                                                                                                 tCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGAAGCAAGCATTATATAATGCCAGGCCAGTTCTTCTTCAATTg         >  3:382835/1‑148 (MQ=255)
                                                                                                                                  ccGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGAAGCAAGCATTATATAATGCCAGGCCAGtt                     <  2:236025/135‑1 (MQ=255)
                                                                                                                                  ccGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGAAGCAAGCATTATATAATGCCAGGCCAGtt                     >  1:236025/1‑135 (MQ=255)
                                                                                                                                    gTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGAAGCAAGCATTATATAATGCCAGGCCAGTTCTTCTTCAATTGTcc      <  4:382835/148‑1 (MQ=255)
                                                                                                                                     ttAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTg                                                                                                        <  1:641677/49‑1 (MQ=255)
                                                                                                                                     ttAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTg                                                                                                        >  2:641677/1‑49 (MQ=255)
                                                                                                                                     ttAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTg                                                                                            >  4:827238/1‑61 (MQ=255)
                                                                                                                                     ttAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTg                                                                                            <  3:827238/61‑1 (MQ=255)
                                                                                                                                     ttAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTc                                                                          <  4:1089016/79‑1 (MQ=255)
                                                                                                                                     ttAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTc                                                                          >  3:1089016/1‑79 (MQ=255)
                                                                                                                                     ttAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGAAGCAAGCATTATATAATGCCAGGCCAGTTCTTCTTCAAt           <  2:171334/142‑1 (MQ=255)
                                                                                                                                     ttAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGAAGCAAGCATTATATAATGCCAGGCCAGTTCTTCTTCAAt           >  1:171334/1‑142 (MQ=255)
                                                                                                                                        aCTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGAAGCAAGCATTATATAATGCCAGGCCAGTTCTTCTTCAATTGTCCCGtt  <  4:300589/148‑1 (MQ=255)
                                                                                                                                           ttCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGaa                                                 >  2:274997/1‑98 (MQ=255)
                                                                                                                                           ttCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGaa                                                 <  1:274997/98‑1 (MQ=255)
                                                                                                                                                 gATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCtt                                                                            <  2:153133/65‑1 (MQ=255)
                                                                                                                                                 gATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCtt                                                                            >  1:153133/1‑65 (MQ=255)
                                                                                                                                                 |                                                                                                                                          
ATTTCAGGCGTTTATGAGTATTTAACGGATGATGCTCCCCACGGAACATTTCTTATGGGCCAACGGCATTTCTTACTGTAGTGCTCCCAAAACTGCTTGTCGTAACGATAACACGCTTCAAGTTCAGCATCCGTTAACTTTCTGCGATAGCAGCAGATATGCCAGTAAAGAAATCCCATTTGACTATTTTTTTGATAATCTTCTTCGCTTTCGAACAACTCGTGCGCCTTTCGAGAAGCAAGCATTATATAATGCCAGGCCAGTTCTTCTTCAATTGTCCCGTT  >  NC_000913/1843820‑1844103

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