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.9% 292.1 / NA 88intergenic (+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 (41/45);  minor base T (0/1);  total (41/47)
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
*NC_000913_3_pae_pgi4,295,6402.G98.9% 292.9 / ‑3.3 88intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (41/46);  total (41/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑CC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTT  >  NC_000913_3_pae_pgi/4295540‑4295780
                                                                                                     ||                                                                                                                                              
cATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑CCAAGCAGCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCa                                                                                               <  1:1658674/150‑1 (MQ=255)
               cGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCg                                                                                >  4:167256/1‑150 (MQ=9)
                  cGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTg                                                                             <  3:1155936/150‑1 (MQ=255)
                         cgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtt                                                                      <  3:371943/150‑1 (MQ=255)
                          gcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttt                                                                     >  3:559343/1‑150 (MQ=255)
                           cgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtttt                                                                    <  1:1549118/150‑1 (MQ=255)
                            gTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttttt                                                                   <  4:686254/150‑1 (MQ=255)
                             tCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc                                                                  <  4:828665/150‑1 (MQ=255)
                                  tCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAct                                                             >  2:692922/1‑150 (MQ=255)
                                   cAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActc                                                            >  4:314868/1‑150 (MQ=255)
                                   cAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActc                                                            >  2:1835314/1‑150 (MQ=255)
                                     ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                                                          <  4:479716/150‑1 (MQ=255)
                                      gCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCCTCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTc                                                         <  2:1211809/150‑1 (MQ=255)
                                      gCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTc                                                         >  1:1431552/1‑150 (MQ=255)
                                      gCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTc                                                         <  2:1155191/150‑1 (MQ=255)
                                       cgtACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTca                                                        <  2:1953456/148‑1 (MQ=255)
                                       ccTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTca                                                        >  1:1163834/1‑150 (MQ=255)
                                         tACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcaca                                                      >  2:823704/1‑150 (MQ=255)
                                          aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCacac                                                     <  2:445247/150‑1 (MQ=255)
                                           cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc                                                    >  1:188788/1‑150 (MQ=255)
                                           cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc                                                    <  1:1835314/150‑1 (MQ=255)
                                            gCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCTCGGCTTATCATGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAccc                                                   >  1:1486314/1‑150 (MQ=17)
                                             ccAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCg                                                  >  2:1589503/1‑150 (MQ=255)
                                              cAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCgg                                                 >  4:1279602/1‑150 (MQ=255)
                                                gACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTa                                               <  4:1182400/150‑1 (MQ=255)
                                                gACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTa                                               <  3:529519/150‑1 (MQ=255)
                                                 aCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTAt                                              <  3:824838/150‑1 (MQ=255)
                                                 aCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTAt                                              <  2:1265857/150‑1 (MQ=255)
                                                   aGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCa                                            >  2:424103/1‑150 (MQ=255)
                                                    gcgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaa                                           <  2:167086/150‑1 (MQ=255)
                                                      gcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAc                                         >  2:116737/1‑150 (MQ=255)
                                                       cAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTCATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAcc                                        <  4:650856/150‑1 (MQ=255)
                                                        cgcaGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCGCACCCGGTATCAAAccc                                       <  2:1742353/147‑1 (MQ=18)
                                                        aaTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAccc                                       >  4:893873/1‑150 (MQ=255)
                                                        aaTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTCTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTc                                                         >  1:2035946/1‑132 (MQ=17)
                                                        aaTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTCTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTc                                                         <  2:2035946/132‑1 (MQ=17)
                                                         aTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCt                                      >  1:272557/1‑150 (MQ=255)
                                                          tAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCtt                                     <  1:2034205/150‑1 (MQ=255)
                                                          tAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCtt                                     <  3:1068412/150‑1 (MQ=255)
                                                          tAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCtt                                     <  1:2034086/150‑1 (MQ=255)
                                                           aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTc                                    <  2:1163834/150‑1 (MQ=255)
                                                           aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTc                                    <  3:173547/150‑1 (MQ=255)
                                                             ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTg                                                                             >  1:1850447/1‑107 (MQ=255)
                                                             ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTg                                                                             <  2:1850447/107‑1 (MQ=255)
                                                              cTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGtt                                                                              <  1:1456140/105‑1 (MQ=255)
                                                              cTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGtt                                                                              >  2:1456140/1‑105 (MQ=255)
                                                              cTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtttt                                                                    <  2:1987535/115‑1 (MQ=255)
                                                              cTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtttt                                                                    >  1:1987535/1‑115 (MQ=255)
                                                               tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATa                                <  1:462439/150‑1 (MQ=255)
                                                               tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATa                                >  2:295085/1‑150 (MQ=255)
                                                               tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATa                                >  2:1155686/1‑150 (MQ=255)
                                                               tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATa                                <  3:468179/150‑1 (MQ=255)
                                                                gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATAc                               <  3:1279602/150‑1 (MQ=255)
                                                                 aTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACa                              >  2:1698600/1‑150 (MQ=255)
                                                                  tttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAg                             <  3:954939/150‑1 (MQ=255)
                                                                    tAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCt                           <  1:1589503/150‑1 (MQ=255)
                                                                     aGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTc                          <  4:855216/150‑1 (MQ=255)
                                                                      gCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCCGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCa                         <  3:1102110/150‑1 (MQ=255)
                                                                       cGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAg                        >  4:780725/1‑150 (MQ=255)
                                                                       cGTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAg                        >  2:1961652/1‑150 (MQ=255)
                                                                        gTGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAgg                       <  1:1155686/150‑1 (MQ=255)
                                                                         tGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGt                      >  1:1599124/1‑150 (MQ=255)
                                                                         tGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGt                      <  1:714391/150‑1 (MQ=255)
                                                                         tGATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGt                      >  3:460929/1‑150 (MQ=255)
                                                                          gATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTc                                    >  1:1743489/1‑135 (MQ=255)
                                                                          gATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTc                                    <  2:1743489/135‑1 (MQ=255)
                                                                          gATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGtt                     >  2:74444/1‑150 (MQ=255)
                                                                          gATTTTGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGtt                     >  2:199725/1‑150 (MQ=255)
                                                                            ttttGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCa                   >  1:1495886/1‑150 (MQ=255)
                                                                              ttGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAAc                 <  3:167256/150‑1 (MQ=255)
                                                                              ttGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAAc                 <  2:630470/150‑1 (MQ=255)
                                                                               tGTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACACGGTATCAAACCCTTCCATACAGCTCAGGTTCAAcc                >  1:1292404/1‑150 (MQ=255)
                                                                                gTAGGTCGGATAAGGCGTTTATGCC‑‑GCATCCGACATCAACGCCTGATGCGACGCTTAACGCGCCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCa               <  2:1645364/150‑1 (MQ=21)
                                                                                gTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttttt                                                                   >  4:447558/1‑98 (MQ=17)
                                                                                gTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttttt                                                                   <  3:447558/98‑1 (MQ=17)
                                                                                gTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTACCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActct                                                           <  2:1901827/106‑1 (MQ=18)
                                                                                gTAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTACCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActct                                                           >  1:1901827/1‑106 (MQ=18)
                                                                                 tAGGTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAg              <  1:692922/150‑1 (MQ=255)
                                                                                   ggTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAAACAGTg            <  4:908433/150‑1 (MQ=255)
                                                                                    gTCGGATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGc           >  2:2076883/1‑150 (MQ=255)
                                                                                       ggATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAgc        >  2:1812424/1‑150 (MQ=255)
                                                                                        gATAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTACCCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAgcg       <  1:2076883/150‑1 (MQ=255)
                                                                                         aTAAGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAgcgc      >  2:1851857/1‑150 (MQ=255)
                                                                                           aaGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCt    >  2:419649/1‑150 (MQ=255)
                                                                                           aaGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCt    >  2:1755671/1‑150 (MQ=255)
                                                                                            aGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttttt                                                                   <  3:704758/86‑1 (MQ=17)
                                                                                            aGGCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttttt                                                                   >  4:704758/1‑86 (MQ=17)
                                                                                             ggCGTTTACGCC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCttt  >  3:311170/1‑150 (MQ=255)
                                                                                                     ||                                                                                                                                              
CATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑CC‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTT  >  NC_000913_3_pae_pgi/4295540‑4295780

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