Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,296,060 C→T 35.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,0600CT35.7% 53.3 / 56.4 70intergenic (+266/+376)gltP/yjcOglutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO
Reads supporting (aligned to +/- strand):  ref base C (25/20);  new base T (11/14);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

GGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTAT  >  NC_000913/4295921‑4296201
                                                                                                                                            |                                                                                                                                             
ggCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGcc                                                                                                                                            >  1:135328/1‑143 (MQ=17)
 gCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACCTCAACGCCt                                                                                                                                           >  1:953436/1‑143 (MQ=17)
 gCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCt                                                                                                                                           <  1:491508/143‑1 (MQ=32)
 gCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCt                                                                                                                                           >  1:229355/1‑143 (MQ=32)
   gTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATTCCTGa                                                                                                                                         >  1:590576/1‑143 (MQ=17)
   gTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTgc                                                                                                                                         >  1:451919/1‑142 (MQ=17)
    tttATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAt                                                                                                                                        >  1:59289/1‑143 (MQ=32)
         gCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTCTAGGTCGGACAAGGCGTTTCCGCCGCATCCCACCTCCATGCCTGCTGCCAc                                                                                                                                   >  2:831407/1‑143 (MQ=11)
           cGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGCCATCAATGCCTGATGCGACGc                                                                                                                                 >  1:771137/1‑143 (MQ=17)
           cGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGc                                                                                                                                 >  1:698161/1‑143 (MQ=17)
           cGCACCCGACATAAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGc                                                                                                                                 <  1:343143/143‑1 (MQ=12)
             cacccGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCtt                                                                                                                               <  1:803334/140‑1 (MQ=17)
             cATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCtt                                                                                                                               <  1:549503/143‑1 (MQ=17)
                    catccccGCCTGATACGACGCATAACGCGGCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGt                                                                                                                        <  1:1016367/137‑1 (MQ=11)
                       cAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCtt                                                                                                                     <  1:605225/143‑1 (MQ=21)
                        aaCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTa                                                                                                                    >  1:46284/1‑143 (MQ=25)
                            ccTGATGCGTGGCTCGGCGCGTCTTATCAGGCCTACGCCAGACGGCGCACTAGCCTGATTTGGCGTGATTTTGTAGGCCGGATAAGGCGGTTACGCCGCATCCTGCATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAg                                                                                                                <  2:1034348/143‑1 (MQ=11)
                             cTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTg                                                                                                                              >  1:114442/1‑128 (MQ=33)
                             cTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGACTTTGTAGGTCGGATAAGGCGTGTACGCCGCATCCGACAGCAACGCCTGATGCGACGCGTGCCGCGTCCTATCAgg                                                                                                               >  2:501909/1‑143 (MQ=11)
                             cTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATCTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTg                                                                                                                              <  2:114442/128‑1 (MQ=18)
                              ttttGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATACGGCGCTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                              <  1:930483/140‑1 (MQ=255)
                              tGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                              <  1:279301/143‑1 (MQ=38)
                                 tGCGACGCTTGCGCGGTCTTAGCGGGCGTGCGCCAGCCCGCGCAATAGCGTGATGTAGCGGGATGTGGGAGGGCGGATAAGGCGATTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTa                                                                                                           <  2:1107410/143‑1 (MQ=2)
                                  tcGACGCCTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTat                                                                                                          <  1:679601/142‑2 (MQ=16)
                                  gCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGc                                                                                                                                             >  2:421416/1‑108 (MQ=17)
                                  gCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGc                                                                                                                                             <  1:421416/108‑1 (MQ=18)
                                  gCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTat                                                                                                          <  1:322645/143‑2 (MQ=18)
                                  gCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTat                                                                                                          >  1:67853/1‑142 (MQ=18)
                                  gCGACGCTTAACGCGTCTTAGCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTat                                                                                                          <  1:791231/143‑2 (MQ=16)
                                    gccGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCtctct                                                                                                        <  1:814040/141‑5 (MQ=255)
                                    gACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAGATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCtctct                                                                                                        >  1:750369/1‑139 (MQ=255)
                                    gACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCAGATGCGACGCTTGCCGCGTCTTGTCAGGCCTAAGc                                                                                                        >  1:354622/1‑143 (MQ=255)
                                     aCGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGGTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatctt                                                                                                       >  1:474426/1‑139 (MQ=9)
                                     aCGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATACGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCTGGccactcct                                                                                                       >  2:1005413/1‑137 (MQ=11)
                                      cGCTTAACGCGTCTTATCAGGCCTACGCTAGACAGCGCAATAGCTTGATTTAGCGTGATTTTGTAGGTCTGATAAGGTGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCa                                                                                                      <  1:111877/143‑1 (MQ=255)
                                         ttAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGctcaac                                                                                                   >  2:888759/1‑138 (MQ=255)
                                             cgcgTCTTATCAGGCCTGCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAgcg                                                                                               <  1:128703/143‑1 (MQ=38)
                                                 gTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGc                                                                                                        <  1:766813/131‑1 (MQ=40)
                                                     tATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGc                                                                                        <  1:815931/143‑1 (MQ=38)
                                                      aTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGcc                                                                                       >  2:832641/1‑143 (MQ=37)
                                                       tCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCa                                                                                                                 <  2:1089049/116‑1 (MQ=25)
                                                       tCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCa                                                                                                                 >  1:1089049/1‑116 (MQ=40)
                                                            ccTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGa                                                                                                                                         >  1:664316/1‑87 (MQ=33)
                                                            ccTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAg                                                                                                 >  2:67129/1‑127 (MQ=38)
                                                                  ccAGACGGCGCAAGAGCCTGATTTAGCGTGATTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTg                                                                           <  2:451919/143‑1 (MQ=255)
                                                                        aGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGAt                                                                                                                                        <  1:879874/76‑1 (MQ=25)
                                                                          ctcAATAGCCTGACTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTCCGCCGCATCCTACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTa                                                                   <  1:1037797/141‑1 (MQ=255)
                                                                          cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGATCGCGCGATCGTCTGGTGTAGCgagagtgtcta                                                                   >  2:814040/1‑136 (MQ=11)
                                                                             aaTAGCCGGATTTCGGGTGATTTTGTGGGTCGGGCAGGGCGTTTACGCCGCGTCCGACATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGt                                                                <  2:410121/143‑1 (MQ=11)
                                                                              aTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCg                                                                                                                                     <  2:231139/73‑1 (MQ=25)
                                                                              aTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCg                                                                                                                                     >  1:231139/1‑73 (MQ=35)
                                                                                           gCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGt                                                                                                                        >  2:660560/1‑73 (MQ=19)
                                                                                           gCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGt                                                                                                                        <  1:660560/73‑1 (MQ=255)
                                                                                            cGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATGAGGCGTTTa                                                 >  1:835279/1‑143 (MQ=38)
                                                                                               gATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGc                                                                                                                                 <  1:1711/60‑1 (MQ=38)
                                                                                                 ttttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGc                                                                                                                                      <  1:656019/53‑1 (MQ=255)
                                                                                                   ttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTg                                                                           >  2:43986/1‑110 (MQ=255)
                                                                                                   ttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTg                                                                           <  1:43986/110‑1 (MQ=25)
                                                                                                        ggTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGGCGCTTGACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTCATTTAGCGTGATTTTGGCGGTCGGATCGGGCGTTTCAGCCgagccaga                                     >  2:213135/1‑136 (MQ=11)
                                                                                                        cgTCCGAGAAGGGGTTGACGCCCGATCCGACATCAACGCCTTATGCCACGCTTGCCGCGTCTTATCAGCCCTACGCGAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTATGTCGGATAATGCGTTTATGCCGCATCCGa                                     <  1:865572/142‑1 (MQ=11)
                                                                                                          tCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCCGACAGCGCGATAGCCTGCTTTAGCGGGATTTTGTAGGTCGGATAAGGCGCTTATTCCGCATCCGCCa                                   >  1:135338/1‑143 (MQ=11)
                                                                                                          tCGGACAAGGCGTTTCCGCCGCACCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACa                                   <  1:395457/143‑1 (MQ=255)
                                                                                                             gATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCGTGATTTAGCGTGATTTTGTAGGCCGGATAAGGCGTTTATGCGGCATCCGACATCa                                >  1:1075750/1‑143 (MQ=255)
                                                                                                              aTAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCaa                               >  1:1008233/1‑143 (MQ=25)
                                                                                                                          cgccgcCTCCGGCATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGa                   <  2:209058/143‑1 (MQ=25)
                                                                                                                            ccgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACg                 >  2:75077/1‑143 (MQ=18)
                                                                                                                                    gACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCTAGACAGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATTACCTCCAGATGCGACGCTTAAcgc         >  2:1010702/1‑143 (MQ=11)
                                                                                                                                    gACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATGTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAAcgc         >  2:549503/1‑143 (MQ=25)
                                                                                                                                     aCATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCGGCATCCGACAGGGAGGCCTGAGTCGAGGCTGCAcgcg        >  1:222688/1‑143 (MQ=11)
                                                                                                                                           aCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTAt  <  1:404730/143‑1 (MQ=31)
                                                                                                                                            |                                                                                                                                             
GGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTAT  >  NC_000913/4295921‑4296201

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: