Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_pgi 4,295,641 +GC intergenic (+587/+55) 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_000913_3_pae_pgi4,295,6401.C98.6% 231.0 / ‑5.6 72intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  major base C (38/32);  minor base T (0/1);  total (38/34)
Fisher's exact test for biased strand distribution p-value = 4.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_000913_3_pae_pgi4,295,6402.G98.6% 235.5 / NA 72intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  major base G (37/33);  minor base T (1/0);  total (38/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 = 1.00e+00

ATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTA‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCC  >  NC_000913_3_pae_pgi/4295534‑4295777
                                                                                                             ||                                                                                                                                         
aTCCGACATTAACGCCTTCTGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACTGCGTTATTGCTTCATTTACCATGGTTATCTACATCGAAGCCGACGCTCTTG‑‑CCCTACCCGACACCCATGCCTGATGCGACGCTTGCCGCGTc                                                                                                  <  2:1147765/150‑1 (MQ=0)
                   tGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTATGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAc                                                                               <  4:512017/150‑1 (MQ=11)
                     agACGCGTTAAGCGTCGCATCAGGCGTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCg                                                                             >  2:220871/2‑150 (MQ=11)
                      gACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGt                                                                            >  3:523965/1‑150 (MQ=12)
                       aCGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGtt                                                                           <  4:690505/150‑1 (MQ=12)
                         gCTTGACGCGTCTTCTCAGGCCTACGCCAGAAAGCGCCATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCAGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTgg                                                                         <  1:305897/150‑1 (MQ=11)
                         gCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGGTAATTCG‑G‑TCACTCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTgg                                                                         >  1:548366/1‑150 (MQ=11)
                         gCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTgg                                                                         >  4:696772/1‑150 (MQ=255)
                          cTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGt                                                                        <  4:495688/150‑1 (MQ=255)
                             aaCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTaa                                                                     >  1:115040/1‑150 (MQ=255)
                              aCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAt                                                                    >  2:1013188/1‑150 (MQ=255)
                               cgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtt                                                                   >  3:625468/1‑150 (MQ=255)
                               cgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtt                                                                   >  1:157274/1‑150 (MQ=255)
                                 cgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACCTTAACCGTTGGTTAAtttt                                                                 >  2:985051/1‑150 (MQ=255)
                                   tCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc                                                               <  3:567847/150‑1 (MQ=255)
                                    cTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCa                                                              <  3:71589/150‑1 (MQ=255)
                                    cTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCa                                                              <  2:826773/150‑1 (MQ=255)
                                    cTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCa                                                              >  4:558488/1‑150 (MQ=255)
                                     ttATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaa                                                             >  2:230128/1‑150 (MQ=255)
                                      tATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa                                                            <  4:208198/150‑1 (MQ=255)
                                       aTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc                                                           >  2:577920/1‑150 (MQ=255)
                                       aTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc                                                           >  1:87536/1‑150 (MQ=255)
                                        tCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAct                                                          >  2:147537/1‑150 (MQ=255)
                                        tCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAct                                                          >  4:194656/1‑150 (MQ=255)
                                        tCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAACTTTTCAAAct                                                          >  1:874556/1‑150 (MQ=255)
                                         cAGGGCTCCCCCCGACAGCGAAATATGCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActc                                                         <  4:537622/150‑1 (MQ=17)
                                          aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGCTAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActct                                                        <  1:425849/150‑1 (MQ=255)
                                          aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc                                                               <  2:1037793/143‑1 (MQ=255)
                                          aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc                                                               >  1:1037793/1‑143 (MQ=255)
                                          aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActct                                                        >  1:1012337/1‑150 (MQ=255)
                                            gCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTc                                                      >  3:473515/1‑150 (MQ=255)
                                            gCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGACGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTc                                                      >  4:8602/1‑150 (MQ=255)
                                             ccTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTca                                                     <  1:382754/150‑1 (MQ=255)
                                                  gCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAccc                                                >  1:334956/1‑150 (MQ=255)
                                                  gCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAccc                                                >  1:809280/1‑150 (MQ=255)
                                                  gCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAccc                                                >  4:500197/1‑150 (MQ=255)
                                                    cAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTCTCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCgg                                              >  1:194725/1‑150 (MQ=255)
                                                     aGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGt                                             >  3:703214/1‑150 (MQ=255)
                                                      gACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTa                                            >  4:551365/1‑150 (MQ=255)
                                                         aGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCa                                         >  2:69335/1‑150 (MQ=255)
                                                         aGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCa                                         >  3:412104/1‑150 (MQ=255)
                                                         aGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCa                                         >  1:567488/1‑150 (MQ=255)
                                                           cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaga                                       >  1:119493/1‑148 (MQ=255)
                                                           cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                       <  1:1138255/150‑1 (MQ=255)
                                                           cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                       <  2:1158554/150‑1 (MQ=255)
                                                                tAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCtt                                  <  3:171052/150‑1 (MQ=255)
                                                                  gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTcc                                <  4:561280/150‑1 (MQ=255)
                                                                   ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCa                               >  1:15886/1‑150 (MQ=255)
                                                                   ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCa                               >  1:15909/1‑150 (MQ=255)
                                                                   ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCa                               <  2:571211/150‑1 (MQ=255)
                                                                   ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCa                               <  1:1173339/150‑1 (MQ=255)
                                                                   ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCa                               <  1:699093/150‑1 (MQ=255)
                                                                   ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCa                               >  1:769644/1‑150 (MQ=255)
                                                                   ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCa                               <  3:764224/150‑1 (MQ=255)
                                                                      gATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATAc                            <  4:14895/150‑1 (MQ=255)
                                                                       aTTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACa                           <  1:499447/150‑1 (MQ=255)
                                                                       aTTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACa                           <  1:69335/150‑1 (MQ=255)
                                                                           aGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCATACTCTTCACACCCGTTATCAAACCCTTCCATACAGCTc                       >  4:154027/1‑150 (MQ=21)
                                                                              gTGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAgg                    >  1:558077/1‑150 (MQ=255)
                                                                               tGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGt                   <  3:606162/150‑1 (MQ=255)
                                                                               tGATTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGt                   <  4:593898/150‑1 (MQ=255)
                                                                                 aTTTTGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTc                 <  2:1161786/150‑1 (MQ=255)
                                                                                   tttGTAGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCaa               <  2:975633/150‑1 (MQ=255)
                                                                                        aGGTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGt          <  1:537798/150‑1 (MQ=255)
                                                                                         ggTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAGACCCTTCCATACAGCTCAGGTTCAACCAGTg         <  2:119493/150‑1 (MQ=255)
                                                                                         ggTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTg         <  4:412483/150‑1 (MQ=255)
                                                                                         ggTCGGATAAGGCG‑T‑GTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTg         <  2:665862/150‑1 (MQ=255)
                                                                                          gTCGGATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGGGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGttc        >  2:1122679/1‑148 (MQ=37)
                                                                                            cGGATAAGGTG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAg      <  1:68172/150‑1 (MQ=255)
                                                                                              gATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                                                       >  2:598448/1‑99 (MQ=18)
                                                                                              gATAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                                                       <  1:598448/99‑1 (MQ=18)
                                                                                                tAAGGCG‑T‑TTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGcc  >  1:1206653/1‑150 (MQ=255)
                                                                                                             ||                                                                                                                                         
ATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑T‑TTA‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCC  >  NC_000913_3_pae_pgi/4295534‑4295777

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