Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,202,309 T→C 6.4% S677S (AGT→AGC)  yehI → DUF4132 domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,202,3090TC 6.4% 194.3 / 3.9 62S677S (AGT→AGC) yehIDUF4132 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (26/32);  new base C (2/2);  total (28/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.87e-01

ACATCGTCCAAAAGAAACATTTCCGTTAATTTTAAATTTTGGCGCAACAGAACTGGCGCTGCCTGTTGCCCACGTCTGGCGACGTTTCGCGGCGCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTACAAACCGTGGCAAACCGCTGGCAGCGTACAGATGTATGGTCTGCC  >  NC_000913/2202162‑2202453
                                                                                                                                                   |                                                                                                                                                
acaTCGTCCAAAAGAAACATTTCCGTTAATTTTAAATTTTCGCGCAACCGAACTGTCGCTGCCTGTTGCCCACGTCTGGCGACGTTTCGCGGCGCAGCGTGATCTGGCTAGCCAGTGGATTTTACAATGGCCTGAACATACGGCTAGt                                                                                                                                                  <  6:144805/148‑1 (MQ=255)
      tCCAAAAGAAACATTTCCGTTAATTTTAAATTTTGGCGCAACAGAACTGGCGCTGCCTGTTGCCCACGTCTGGCGACGTTTCGCGGCGCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACtt                                                                                                                                            <  7:262781/148‑1 (MQ=255)
       ccAAAAGAAACATTTACGTTAATTTTAAATTTTGTCGCAACAGAACTGGCGCTGACTGTTGCCCACGTCATTCGAATTATCGCGGCGCAGCGTGATCTGGATCTCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACtt                                                                                                                                            >  5:96687/1‑147 (MQ=255)
        cAAAAGAAACATTTCCGTTAATTTTAAATTTTGGCGCAACAGAACTGGCGCTGCCTGTTGCCCACGTCTGGCGACGTTTCGCGGCGCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTAt                                                                                                                                          <  2:201091/148‑1 (MQ=255)
           aaGCAACATTTCCGTTAATTTTAAATTTTGGCGCAACAGAACTGGCGCTGCCTGTTGCCCACGTCTGGCGACGTTTCGCGGCGCAGCGTGATCTGGCTCGCCATTGGAATTTACAATGGCCGGACCATACGGCTAGTGCACTTATccc                                                                                                                                       >  2:259842/1‑148 (MQ=255)
           aaGAAACATTTCCGTTAATTTTAAATTTTGGCGCAACAGAACTGGCGCTGCCTGTTGCCCACGTCTGGCGACGTTTCGCGGCGCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATccc                                                                                                                                       >  7:203838/1‑148 (MQ=255)
             gAAACATTTCCGTTAATTTTAAATTTTGGCGCAACAGAACTGGCGCTGCCTGTTGCCCACGTCTGGCGACGTTTCGCGGCGCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCctc                                                                                                                                     <  3:118434/148‑1 (MQ=255)
              aaccATTTCCGTTAATTTTAAATTTTGGCGCAACAGAACTGGCGCTGCCTGTTGCCCACGTCTGGCGACGTTTCGCGGCGCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCctct                                                                                                                                    <  7:297923/145‑1 (MQ=255)
              aaaCATTTCCGTTAATTTTAAATTTTGGCGCAACAGAACTGGCGCTGCCTGTTGCCCACGTCTGGCGACGTTTCGCGGCGCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCctct                                                                                                                                    <  2:231696/148‑1 (MQ=255)
                    ttCCGTTAATTTTAAATTTTGGCGCAACAGAACTGGCGCTGCCTGTTGCCCACGTCTGGCGACGTTTCGCGGCGCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCtt                                                                                                                              >  7:122499/1‑148 (MQ=255)
                     tacGTTAATTTTAAATTTTGGCTCAACAGAACTGGCGCTGACTGTTGCCCACGTCTGGCGACGTTTCGCGGCGCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATAAGGCTAGTGCACTTATCCCTCTTGTCtt                                                                                                                              <  1:259842/145‑1 (MQ=255)
                          tAATTTTAAATTTTGGCGCAACAGAACTGGCGCTGCCTGTTGCCCACGTCTGGCGACGTTTCGCGGCGCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCa                                                                                                                         >  1:317849/1‑147 (MQ=255)
                           aaTTTTAAATTTTGGCGCAACAGAACTGGCGCTGCCTGTTGCCCACGTCTGGCGACGTTTCGCGGCGCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCaa                                                                                                                        <  7:107970/147‑1 (MQ=255)
                            aTTTTAAATTTTGGCGCAACAGAACTGGCGCTGCCTGTTGCCCACGTCTGGCGACGTTTCGCGGCGCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATGCGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAAc                                                                                                                      <  8:292329/148‑1 (MQ=255)
                             ttttAAATTTTGGCGCAACAGAACTGGCGCTGCCTGTTGCCCACGTCTGGCGACGTTTCGCGGCGCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCctct                                                                                                                                    <  4:91828/133‑1 (MQ=255)
                             ttttAAATTTTGGCGCAACAGAACTGGCGCTGCCTGTTGCCCACGTCTGGCGACGTTTCGCGGCGCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCctct                                                                                                                                    >  3:91828/1‑133 (MQ=255)
                              tttAAATTTTGGCGCAACAGAACTGGCGCTGCCTGTTGCCCACGTCTGGCGACGTTTCGCGGCGCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAAccc                                                                                                                    >  5:70866/1‑148 (MQ=255)
                               ttAAATTTTGGCGCAACAGAACTGGCGCTGCCTGTTGCCCACGTCTGGCGACGTTTCGCGGCGCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAAcc                                                                                                                     >  3:222822/1‑146 (MQ=255)
                                              aCAGAACTGGCGCTGCCTGTTGCCCACGTCTGGCGACGTTTCGCGGCGCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAg                                                                                                    >  4:22717/1‑148 (MQ=255)
                                              aCAGAACTCGCGCTGCCTGTTGACCACGTCTGGCGACGTATCGCGGCGCAGCGTGATCTGGATCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCAAGCGATAATATCGAAg                                                                                                    >  6:107347/1‑148 (MQ=255)
                                                      ggCGCTGCCTGTTGCCCACGTCTGGCGACGTTTCGCGGCGCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGCGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGat                                                                                                              >  5:193006/1‑130 (MQ=255)
                                                      ggCGCTGCCTGTTGCCCACGTCTGGCGACGTTTCGCGGCGCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGCGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGat                                                                                                              >  5:188656/1‑130 (MQ=255)
                                                      ggCGCTGCCTGTTGCCCACGTCTGGCGACGTTTCGCGGCGCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGCGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGat                                                                                                              <  6:188656/130‑1 (MQ=255)
                                                      ggCGCTGCCTGTTGCCCACGTCTGGCGACGTTTCGCGGCGCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGCGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGat                                                                                                              <  6:193006/130‑1 (MQ=255)
                                                       gcgcTGCCTGTTGCCCACGTCTGGCGACGTTTCGCGGCGCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTa                                                                                            <  2:317849/147‑1 (MQ=255)
                                                          cTGCCTGTTGCCCACGTCTGGCGACGTTTCGCGGCGCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTg                                                                                        >  1:48602/1‑148 (MQ=255)
                                                           tGCCTGTTGCCCACGTCTGGCGACGTTTCGCGGCGCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGc                                                                                       >  2:24265/1‑148 (MQ=255)
                                                             ccTGTTGCCCACGTCTGGCGACGTTTCGCGGCGCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGccc                                                                                     <  8:203838/148‑1 (MQ=255)
                                                             ccTGTTGCCCACGTCTGGCGACGTTTCGCGGCGCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGcc                                                                                      <  5:307177/147‑1 (MQ=255)
                                                             ccTGTTGCCAACGTCTGGCGTCGTTTCGCGGCGCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTTCTTGcc                                                                                      <  6:96687/147‑1 (MQ=255)
                                                                gTTGCCCACGTCTGGCGACGTTTCGCGGCGCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGc                                                                                  >  8:37468/1‑148 (MQ=255)
                                                                  tGCCCACGTCTGGCGACGTTTCGCGGCGCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGt                                                                                >  5:57425/1‑148 (MQ=255)
                                                                              gCGACGTTTCGCGGCGCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGa                                                                     <  4:222822/147‑1 (MQ=255)
                                                                               cGACGTTTCGCGGCGCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAAc                                                                   <  1:6702/148‑1 (MQ=255)
                                                                                 acatttCGCGGCGCAGCGTGATGTGGCTCGCCTGTGGATTTTTCAATGTCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCTAACCCAGCGATATTTGCGAAGCCGCATTACTTGCCCGGCGTTTACTGTACTAACat                                                                 <  4:208997/145‑2 (MQ=255)
                                                                                  cGTTTCGCGGCGCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAgg                                                                <  1:24265/148‑1 (MQ=255)
                                                                                   gTTTCGCGGCGCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAggg                                                               <  6:57425/148‑1 (MQ=255)
                                                                                      tCGCGGCGCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGAc                                                             <  8:1889/147‑1 (MQ=255)
                                                                                      tCGCGGCGCAGCGTGATCTGGCTCGCCAGTGGAATTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGAc                                                             <  1:76149/147‑1 (MQ=255)
                                                                                          ggcgCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGc                                                        >  8:229936/1‑148 (MQ=255)
                                                                                          ggcgCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGc                                                        >  2:7551/1‑148 (MQ=255)
                                                                                            cgcAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGa                                                      <  7:37468/148‑1 (MQ=255)
                                                                                                  gTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCt                                                <  7:20110/148‑1 (MQ=255)
                                                                                                    gATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTa                                               >  2:193348/1‑147 (MQ=255)
                                                                                                       cTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTACaaa                                           >  8:136039/1‑148 (MQ=255)
                                                                                                          gCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTACAAACCg                                        <  3:22717/148‑1 (MQ=255)
                                                                                                            tCGCCAGTGGATTTTACTTTGGCCGGAACATACGGCTTGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCAGTACTTGGCGTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTACAAACCGTg                                      <  5:107347/148‑1 (MQ=255)
                                                                                                              gCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTACAAACCGTgg                                     <  2:48602/147‑1 (MQ=255)
                                                                                                                    ggATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTACAAACCGTGGCAAACCg                              >  7:271143/1‑148 (MQ=255)
                                                                                                                           aCAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTACAAACCGTGGCAAACCGCTGGCa                        <  8:122499/147‑1 (MQ=255)
                                                                                                                            cAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTACAAACCGTGGCAAACCGCTGGCAGc                      >  7:244524/1‑148 (MQ=255)
                                                                                                                             aaTGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTACAAACCGTGGCAAACCGCTGGCAGCg                     <  1:193348/148‑1 (MQ=255)
                                                                                                                             aaTGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTACAAACCGTGGCAAACCGCTGGCAGCg                     >  3:95524/1‑148 (MQ=255)
                                                                                                                              aTGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTACAAACCGTGGCAAACCGCTGGCAGCg                     >  3:223462/1‑147 (MQ=255)
                                                                                                                              aTGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTACAAACCGTGGCAAACCGCTGGCAGCg                     >  3:219537/1‑147 (MQ=255)
                                                                                                                                   cGGAACATACGGCTAGTGCACTTATCCATCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTACAAACCGTGGCAAACCGCTGGCAGCGTACAGa               <  7:136039/148‑1 (MQ=255)
                                                                                                                                    ggAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTACAAACCGTGGCAAACCGCTGGCAGCGTACAGAt              >  3:122300/1‑148 (MQ=255)
                                                                                                                                           aCGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTACAAACCGTGGCAAACCGCTGGCAGCGTACAGATGTATGGt       >  1:107571/1‑148 (MQ=255)
                                                                                                                                            cGGCTAGTGCACTTCTCCCTCTTGTCTTTATCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTACAAACCGTGGCAAACCGCTGGCAGCGTACAGATGTATGGTc      <  5:17168/148‑1 (MQ=255)
                                                                                                                                                 aGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAAc                                                                   >  7:103418/1‑82 (MQ=255)
                                                                                                                                                 aGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAAc                                                                   <  8:103418/82‑1 (MQ=255)
                                                                                                                                                  gTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTACAAACCGTGGCAAACCGCTGGCAGCGTACAGATGTATGGTCTGcc  <  7:229936/146‑1 (MQ=255)
                                                                                                                                                  gTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTACAAACCGTGGCAAACCGCTGGCAGCGTACAGATGTATGGGCTGcc  <  1:7551/146‑1 (MQ=255)
                                                                                                                                                   tGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTACAAACCGTGGCAAACCGc                             >  7:44442/1‑118 (MQ=255)
                                                                                                                                                   tGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTACAAACCGTGGCAAACCGc                             <  8:44442/118‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
ACATCGTCCAAAAGAAACATTTCCGTTAATTTTAAATTTTGGCGCAACAGAACTGGCGCTGCCTGTTGCCCACGTCTGGCGACGTTTCGCGGCGCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTACAAACCGTGGCAAACCGCTGGCAGCGTACAGATGTATGGTCTGCC  >  NC_000913/2202162‑2202453

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑