Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,295,972 T→G 7.6% intergenic (+178/+464) gltP → / ← yjcO glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,295,9720TG 7.6% 142.2 / 5.0 69intergenic (+178/+464)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base T (35/28);  new base G (2/3);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

CAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCT  >  NC_000913/4295823‑4296117
                                                                                                                                                     |                                                                                                                                                 
cAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCt                                                                                                                                                   <  2:1609462/150‑1 (MQ=39)
     tgtgAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGGTTAACTCGTCTTATct                                                                                                                                              >  1:1226687/1‑149 (MQ=25)
     tgtgAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTTTCa                                                                                                                                              >  1:1157200/1‑150 (MQ=25)
     tgtgAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCa                                                                                                                                              <  1:1828341/150‑1 (MQ=35)
      gtgAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATTGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAg                                                                                                                                             <  2:1196417/150‑1 (MQ=21)
      gtgAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGATTATGCCGCCAACGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAg                                                                                                                                             >  2:1095741/1‑150 (MQ=255)
       tgAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAgg                                                                                                                                            <  1:1309685/150‑1 (MQ=33)
        gAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGc                                                                                                                                           <  2:733427/150‑1 (MQ=32)
        gAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATAAGACGCGTTAAGCGTCGCATCAGGc                                                                                                                                           <  1:1087934/150‑1 (MQ=11)
          aTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATAAGACGCGTTAAGCGTCGCATCAGGcgt                                                                                                                                         <  1:499389/150‑3 (MQ=11)
          aTGCCTGATGCGACGCTGGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCTGATATGGCGTTTATGCCGCATCCGACATCAACGCCTGATAAGACGCGTTAAGCGTCGCATCAGGcgt                                                                                                                                         <  1:1202480/150‑3 (MQ=11)
                                                                        ttAGCGTGATTTTGTATGTCGGATTTGGCGTTTATTCCGCATCCGTCATCAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAcgc                                                                           <  1:421735/150‑1 (MQ=11)
                                                                        ttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAcgc                                                                           <  1:572460/150‑1 (MQ=31)
                                                                         tAGCGTGATTTTGTAGTGCTTATACGGCGTTTATGCCGCTTCCGTCTTCAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAcgcc                                                                          <  2:1443498/150‑1 (MQ=11)
                                                                         tAGCGTGATTTTGTAGGTCTGATAAGGGGTTTAAGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAcgcc                                                                          <  1:493128/150‑1 (MQ=12)
                                                                         tAGCGTGATTTTGTAGGGCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAcgcc                                                                          <  1:1266309/150‑1 (MQ=17)
                                                                          aGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGACTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAcgccg                                                                         >  2:393402/1‑150 (MQ=17)
                                                                           gCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAcgccgc                                                                        >  2:577517/1‑150 (MQ=32)
                                                                             gTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGCCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAt                                                                      >  2:1776909/1‑150 (MQ=17)
                                                                                aTTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCg                                                                   >  1:1401226/1‑150 (MQ=32)
                                                                                 ttttGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGa                                                                  >  1:460791/1‑150 (MQ=32)
                                                                                    tGTAGCTCGGATAAGGCTTTTATGCCTCATCCGACTTTAACGACTTATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACAt                                                               <  2:1322842/150‑1 (MQ=11)
                                                                                        ggTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCaat                                                           >  2:1795422/1‑149 (MQ=17)
                                                                                         gTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCaatg                                                          >  1:761957/1‑148 (MQ=17)
                                                                                         gTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCaatg                                                          >  2:999896/1‑148 (MQ=17)
                                                                                         atCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTAAGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTATGTCGGATAAGGCGATTACGACGCATCCGACATCaaag                                                          >  1:1317796/2‑148 (MQ=11)
                                                                                          tCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGc                                                         >  2:819976/1‑150 (MQ=17)
                                                                                            ggATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCt                                                       >  1:1819015/1‑150 (MQ=17)
                                                                                             gATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCATAAGCCTGATATAGCGTGAATTAGTAGATCGGATGAGGCGTTTACGCCGCATCGGCAATCAATGCCTg                                                      >  2:213066/1‑150 (MQ=11)
                                                                                              ttAAGGCGTTTATGCCGCATCAGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGa                                                     <  2:1048347/149‑1 (MQ=17)
                                                                                              aaaaGGCGTTTATGCCGCATCCGAAATCAACGCCTGATGCGACGCTTCACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGa                                                     <  2:1018339/148‑1 (MQ=12)
                                                                                                 aGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGc                                                  >  2:1802797/1‑150 (MQ=17)
                                                                                                  ggCGTTTATGCCGCTTCCGACATGAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCg                                                 <  2:441004/150‑1 (MQ=17)
                                                                                                          tGCCGCATCCGACATCAACGCCTGATGCGACGCTTCCCGCGTCCTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGc                                         <  2:1285352/150‑1 (MQ=12)
                                                                                                           gCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGcc                                        >  1:1648679/1‑150 (MQ=21)
                                                                                                           gCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGACTGATTTAGCGTGATTTTGTAGGTCGGATACGGCGTTAACGACGTATCCGAAATAAGCACCAGATGCGTCGCtggaa                                        >  2:1143697/1‑146 (MQ=11)
                                                                                                             cGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCcgc                                      >  2:1414944/1‑150 (MQ=25)
                                                                                                             cGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGACTGATGCGACGCTTGCcgc                                      >  2:1422003/1‑150 (MQ=25)
                                                                                                              gCATCCGTCATCAACGCCTGATGCGACGCTTCACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATATCCTGATTTAGCGTGATTTTTTAGGTCTGTTAAGGCTTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCcgcg                                     <  2:1233966/150‑1 (MQ=11)
                                                                                                                aTCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTc                                   <  1:518730/150‑1 (MQ=25)
                                                                                                                aTCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTc                                   >  2:787354/1‑150 (MQ=25)
                                                                                                                  ccGACATCAACGCCTGATGCGACGCTTCACGCGTCTTATCAGGCCTACGCCAGACATCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCtt                                 <  2:4905/150‑1 (MQ=255)
                                                                                                                      cATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCGAGACGGCGCTATAGCGTGATTTAGCGTGAGTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGCCATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCa                             >  1:165800/1‑150 (MQ=11)
                                                                                                                        tCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                           >  1:1758838/1‑150 (MQ=37)
                                                                                                                         cAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                          >  1:1073076/1‑150 (MQ=25)
                                                                                                                         cAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGAAAAAGCGTGAATATCTAGGTCGGATAAGGCGTTTACGCCGCATCCCACATCAACGCCAGATGCGACGCTTGACGCGTCTTATCAgcc                          >  2:1166996/1‑148 (MQ=11)
                                                                                                                         cAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGACAGACAGCGCAATAGCCTGATTTAGCGTGATAAAGTAGTTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGATTGCAGCGTCTTAACAGGc                          >  2:1201070/1‑150 (MQ=11)
                                                                                                                          aatgCCTGATGCTACGCTTAACTGGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATATGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                         <  2:1492210/147‑1 (MQ=11)
                                                                                                                          aaCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                         >  2:966218/1‑150 (MQ=25)
                                                                                                                          aaCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                         <  1:393402/150‑1 (MQ=25)
                                                                                                                            ggCCTGATAAGACGCGTTAAGCGTCGCATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTa                       >  1:1851177/2‑150 (MQ=11)
                                                                                                                            cGCCTGATTCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACTGCTCAATTGCCTGATTTAGCGTGATTTTGTAGGTCGGTTAAGGCGTTTACGCCGCATCCGTCATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTa                       <  1:578883/150‑1 (MQ=11)
                                                                                                                             gCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTat                      >  1:965988/1‑149 (MQ=21)
                                                                                                                              ccTGATAAGACGCGTTAAGCGTCGCATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatc                     >  2:233762/1‑148 (MQ=11)
                                                                                                                               cTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatct                    <  2:1607848/150‑4 (MQ=17)
                                                                                                                               cTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatct                    <  1:1467471/150‑4 (MQ=17)
                                                                                                                                tGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCt                                                       <  1:414467/114‑1 (MQ=17)
                                                                                                                                tGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCt                                                       >  2:414467/1‑114 (MQ=17)
                                                                                                                                tGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatctt                   >  2:1483709/1‑146 (MQ=17)
                                                                                                                                tGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatctt                   >  2:648678/1‑146 (MQ=17)
                                                                                                                                tGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGcc                   <  1:1850236/150‑1 (MQ=37)
                                                                                                                                 gATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatctta                  >  2:1407846/1‑145 (MQ=9)
                                                                                                                                 gATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatctta                  >  2:243921/1‑145 (MQ=9)
                                                                                                                                  aTGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAg                 <  1:1870132/150‑1 (MQ=17)
                                                                                                                                  aTGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAg                 <  1:1631514/150‑1 (MQ=37)
                                                                                                                                              aCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGc                                                         <  2:1502825/98‑1 (MQ=18)
                                                                                                                                              aCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGc                                                         >  1:1502825/1‑98 (MQ=18)
                                                                                                                                                 ggtattATCAGGCCTACTCAAGACAGCTCAATAGCCTGATTTAGCGTGATTTTGTAGGTCTGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATATCCt  <  2:1668746/146‑1 (MQ=11)
                                                                                                                                                 cgTCTTATCAGGCCTACGCCAGACAGCTCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCt  <  1:1403543/150‑1 (MQ=37)
                                                                                                                                                    cTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGc                                                         <  1:92726/92‑1 (MQ=18)
                                                                                                                                                    cTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGc                                                         >  2:92726/1‑92 (MQ=18)
                                                                                                                                                     ttATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGAc                                               <  2:1018740/101‑1 (MQ=17)
                                                                                                                                                     ttATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGAc                                               >  1:1018740/1‑101 (MQ=17)
                                                                                                                                                     |                                                                                                                                                 
CAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCT  >  NC_000913/4295823‑4296117

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