Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,296,060 C→T 36.7% intergenic (+266/+376) gltP → / ← yjcO glutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,296,0600CT36.7% 53.3 / 51.0 73intergenic (+266/+376)gltP/yjcOglutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO
Reads supporting (aligned to +/- strand):  ref base C (30/15);  new base T (14/12);  total (45/28)
Fisher's exact test for biased strand distribution p-value = 3.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATC  >  NC_000913/4295919‑4296202
                                                                                                                                             |                                                                                                                                              
aaGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAGCg                                                                                                                                               >  1:99926/1‑143 (MQ=17)
     gTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCGTCCGTCATCAATGCCTGa                                                                                                                                          >  1:601377/1‑143 (MQ=12)
     gTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGa                                                                                                                                          >  1:466458/1‑143 (MQ=32)
     gTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCCGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGa                                                                                                                                          >  1:221124/1‑143 (MQ=17)
      tttATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCCGGACTACGCCAGAACGCGTGAAACCCTGATCTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCTGACATCAATCCCTGAt                                                                                                                                         >  1:496731/1‑143 (MQ=11)
      tttATGCCGCATCCGACATCAACACCTGAGGCGGCTCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGAt                                                                                                                                         >  1:159847/1‑143 (MQ=11)
            ccGCGTCCGACATCAACGCCGGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAGTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACg                                                                                                                                   <  1:558327/143‑1 (MQ=12)
            ccGCATCCGACATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCGTGATTTAGCGTGATTGTGTAGGTCGGGTAACGCGTTTACGCCGGCTCCGACATCAAGGGCTGCTGCGACg                                                                                                                                   >  2:70142/1‑143 (MQ=11)
            ccGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTCCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGCCCTCCATGCCTGATGCGACg                                                                                                                                   >  1:26835/1‑143 (MQ=11)
            ccGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACCTCAATGCCTGATGCGCCg                                                                                                                                   >  2:61141/1‑143 (MQ=12)
             cGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGTTACGGTGTTGACGCCGCATCCGACATCGATGCCTGATGCTACGc                                                                                                                                  >  1:639277/1‑143 (MQ=11)
                 tcggagaTCAACGCCGGATGCGACGTTTTACGCGTCTTATCAGGCCTACGCCAGACAGCTCAGTAGCCTGATTTAGCGTGATTTGGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGc                                                                                                                              <  2:466458/137‑1 (MQ=11)
                  ccGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGcc                                                                                                                             >  2:85202/1‑143 (MQ=21)
                    gACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCcgc                                                                                                                           >  2:505715/1‑143 (MQ=25)
                       aTCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGAGTTAGCGTGATTTTGTAGGGCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGCCGCTTGcgcctta                                                                                                                        >  2:102840/1‑137 (MQ=11)
                        tCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCt                                                                                                                       >  1:507351/1‑143 (MQ=25)
                        tCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCt                                                                                                                       >  1:661022/1‑143 (MQ=21)
                        tCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCt                                                                                                                       >  1:589917/1‑143 (MQ=21)
                        tCAACGCCTGACGCGACGCTGACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCt                                                                                                                       <  2:633684/143‑1 (MQ=255)
                          aaCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCCATGCCTGATGCGACGCTTGCCGCGTCttc                                                                                                                     >  2:21214/1‑142 (MQ=25)
                           aCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                                                    >  1:326288/1‑143 (MQ=25)
                             gCCTGATGTGACGCTCAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCa                                                                                                                  <  1:186217/143‑1 (MQ=255)
                             gCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGTTTGCCGCGTCTTATCa                                                                                                                  <  2:15527/143‑1 (MQ=25)
                             gCCTGATAAGACGCGTTAAGCGTCGCATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACACCAACGCCTGATGCGACGCTTGCCGCGTCTtctcc                                                                                                                  >  1:149173/1‑142 (MQ=11)
                               ctgctgCGACGCGTTCCGCGTCGTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                                                                                                <  2:173738/143‑1 (MQ=11)
                               ctcatacGACGCGTGTCGCGTCTCATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                                                                                                <  2:149173/137‑1 (MQ=11)
                                tGATGCGGCGGTGAGCGGGGCTTATCGGGCCTAGGGCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                               <  2:221124/143‑1 (MQ=11)
                                tGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGAt                                                                                                                                         <  1:599241/117‑1 (MQ=18)
                                tGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAt                                                                                                                                         >  2:599241/1‑117 (MQ=32)
                                   tGCGACGCTGAAGGTGGGTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTa                                                                                                            <  2:163820/143‑1 (MQ=14)
                                    gCGACGCTTAACGCGTCTTATCAGGCCTCCTCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTat                                                                                                           <  1:596810/143‑2 (MQ=255)
                                    gCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTat                                                                                                           <  1:571811/143‑2 (MQ=21)
                                    gCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAACAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAc                                                                                                           <  1:537286/143‑1 (MQ=255)
                                     cGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACg                                                                                                          <  2:470735/143‑1 (MQ=18)
                                      gACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatgg                                                                                                         >  1:666828/1‑140 (MQ=38)
                                       aCGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGcc                                                                                                        >  1:176181/1‑143 (MQ=38)
                                             aaCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAg                                                                                                  <  1:93440/143‑1 (MQ=38)
                                                      tagcAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGc                                                                                         <  1:229960/140‑1 (MQ=255)
                                                       aTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGGCAGCCGGCGGCATAGcc                                                                                        >  2:563179/1‑143 (MQ=11)
                                                          aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGCAGCTGGCCCCGTCTTATCAGGCCTACGCCAGACACCGCCATCGCCTGa                                                                                     >  1:670568/1‑143 (MQ=11)
                                                          aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                                                                                                >  1:198077/1‑116 (MQ=40)
                                                          aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCTAGCCGGCGCAATAGCCttc                                                                                     >  2:246451/1‑141 (MQ=255)
                                                               tACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGCCGGCGCAATGGCCTCAtttcg                                                                                >  2:364206/1‑141 (MQ=255)
                                                               tACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGACTGCTTTAg                                                                                >  1:366547/1‑143 (MQ=255)
                                                               tACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGGTTACGCCGAATACGACATCATCGCCTGATGCGACGCTTACCGCGTCTTATCCGGCCTACGCCAGCCGGCGCGTTTGCCTGGtttcg                                                                                >  2:409245/1‑141 (MQ=11)
                                                                 cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                              >  2:455050/1‑111 (MQ=25)
                                                                 cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                              <  1:455050/111‑1 (MQ=25)
                                                                  gCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACATCGCAAAAGCCTTCTTTATCGt                                                                             >  2:176950/1‑143 (MQ=11)
                                                                    cAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCACCGCCTGATGCGTCGCTTGCCGCGTCTTCTCAGGCCTGCGCCAGGCCTCGCAACAGCCTGCTCTAGCgcgg                                                                           >  1:417782/1‑140 (MQ=11)
                                                                           cgcAATAGCCTGCTTGAGCGTGATTTTGTAGGTCGGATAAGGCGTATACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTa                                                                    <  2:408610/143‑1 (MQ=255)
                                                                           caacatggCCTGATTTCGGGTGATTTTGTAGGGCGTGTAGGGCGTTTACTCCGCGTCCGGCATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTa                                                                    <  2:310827/136‑1 (MQ=11)
                                                                                             cGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCGTCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTa                                                  <  1:59234/143‑1 (MQ=255)
                                                                                               tGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTCTg                                                >  1:463860/1‑143 (MQ=255)
                                                                                               cggtcttGTAGGGCGGACAAGGCATTTACGCCGCGTCAGACATGAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATg                                                <  2:326915/138‑1 (MQ=11)
                                                                                                  ttttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCtt                                                                                                                      <  2:255164/70‑1 (MQ=255)
                                                                                                  ttttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCtt                                                                                                                      >  1:255164/1‑70 (MQ=255)
                                                                                                      gTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTa                                                  <  2:219687/134‑1 (MQ=40)
                                                                                                      gTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTa                                                  >  1:219687/1‑134 (MQ=40)
                                                                                                         ggTCGGTTAAGGCGGTTCCGCCGCACCCGACATCAAGGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGACGGATAAGGCGTTTATGCCGCATCCGa                                      <  1:31162/143‑1 (MQ=11)
                                                                                                           tCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACa                                    >  1:271191/1‑143 (MQ=25)
                                                                                                                tAAGGCGTTTACGCCGCTTTCGACATCAACTCCTGTTGCGTCGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAAc                               <  2:99926/143‑1 (MQ=11)
                                                                                                                tAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAGGCGACTCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCCATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGCTTATGCCGCATCCGACATCAAc                               >  1:465795/1‑143 (MQ=255)
                                                                                                                  aGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCTCGACGGCGCAATAGCCTGATTTCGCGTGACTTTGTAGGTCGGATACGGCGTTTCTGCCGCATCCGACCCCAAccc                             >  2:365044/1‑141 (MQ=14)
                                                                                                                  aGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGc                             <  2:271191/143‑1 (MQ=21)
                                                                                                                       tttACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCcgcg                                                                                                                          <  1:507909/45‑1 (MQ=255)
                                                                                                                       tttACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCcgcg                                                                                                                          >  2:507909/1‑45 (MQ=255)
                                                                                                                                 aTCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATATCCTGATTTAGCGTGATTTTGTAGGTCGGATCAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCttg              >  1:185534/1‑142 (MQ=255)
                                                                                                                                     gACATCAACGCCTGATGCGACGCTTGCCGCGTCTTCTCAGGCCTACGCCAGACAGCGCCATAGCCTGATTTAGCGCGATTTTGTAGGTCGGCTCCGGCGTTTGTGCCGCTTCCGACATCCACGCCCGATGCGTCGCtttccgc          >  1:575282/1‑138 (MQ=11)
                                                                                                                                     gACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGc                                                                               >  2:45690/1‑74 (MQ=255)
                                                                                                                                     gACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAAcgc          >  2:521394/1‑143 (MQ=35)
                                                                                                                                        aTCAACGCCTGATGCGACGCTTGCCGCGTCTTACCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATACGGCGTTTATGCCGCATCAGACTTCACCTCCCGACGCGACGCTTCACGCGTc       >  1:422476/1‑143 (MQ=11)
                                                                                                                                          cAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTCGCGTGATTTTGTAGGTCGGATACGGCGTTTATGCCGCATCCGACATCAACGCCTGATGGCACGCTTCACCCGTCtt     >  2:285431/1‑143 (MQ=11)
                                                                                                                                            aCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGa                                      <  1:3578/108‑1 (MQ=31)
                                                                                                                                             cGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTCGCGCGATCCTTTAGGTCGGATAAGGCGCTTATGCAGCACCCGACATCAACGCCTGATGCGACGCTTACCGCGTCTTATc  >  2:590707/1‑143 (MQ=11)
                                                                                                                                             |                                                                                                                                              
AAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATC  >  NC_000913/4295919‑4296202

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: