Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_pae_tpiA 4,296,296:1 +CG 100% intergenic (+586/+56) 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_tpiA4,296,2961.C93.7% 156.6 / ‑1.4 60intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (0/2);  major base C (32/20);  minor base T (2/2);  total (34/26)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
*NC_000913_3_pae_tpiA4,296,2962.G95.0% 180.1 / 1.6 60intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (0/3);  new base G (34/23);  total (34/26)
Fisher's exact test for biased strand distribution p-value = 7.60e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑A‑‑CCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCG  >  NC_000913_3_pae_tpiA/4296210‑4296431
                                                                                           ||                                                                                                                                        
gCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAc                                                                             <  1:1099949/149‑1 (MQ=2)
ccGACGCTTAACGCGTCTTATCAGGCCTCCGCCAGCCAGCGCACTAGCCTGATTTAGCGTGATTTTGTGCATTGTATACCGCACGAGG‑‑G‑‑GCGCATC‑CGATATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAc                                                                             <  1:1038993/148‑1 (MQ=2)
    cGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGtt                                                                         <  2:499204/149‑1 (MQ=255)
       ttAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCTTTGGt                                                                      >  2:1102227/1‑149 (MQ=255)
          ccgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCGGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTaa                                                                   <  1:985287/148‑1 (MQ=255)
            gcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtt                                                                 >  2:1117244/1‑149 (MQ=255)
             cgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGGGTCTTATCAGGCATATCTTAACCGTTGGTTAAttt                                                                >  2:1353328/1‑149 (MQ=14)
              gTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ATGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTATCTTAGCCGTCGGATAAtttt                                                               >  2:558589/1‑149 (MQ=11)
              gTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ATGCCGCATC‑CGACATCAACGCCTGATGCGACGCTTCACGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtttt                                                               >  2:558609/1‑149 (MQ=11)
              gTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtttt                                                               >  2:1417205/1‑149 (MQ=255)
               tCTTATCAGGCCTACGCCAGACAGCGCAATCGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGCATCTTATCAGGCCTATCTTACCGTTTGGTTAAtttat                                                              >  1:234798/1‑147 (MQ=11)
               tCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTCTCTTAACCGTTGGTTACttttt                                                              >  2:343118/1‑149 (MQ=11)
                 ttATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTAGTTAATTTTTCa                                                            >  2:102507/1‑149 (MQ=255)
                  tATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTcca                                                           >  2:260253/1‑147 (MQ=255)
                  tATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaa                                                           <  1:133881/149‑1 (MQ=255)
                  tATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaa                                                           >  2:217672/1‑149 (MQ=255)
                   aTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGGAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa                                                          <  1:382420/149‑1 (MQ=255)
                    tCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc                                                         >  1:753253/1‑149 (MQ=255)
                      aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ATGCCGCATC‑CGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActc                                                       <  1:558609/149‑1 (MQ=255)
                      aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActc                                                       <  1:558589/149‑1 (MQ=255)
                      aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGACCTATCTTAACCGTTGGTTAATTTTTCAAActc                                                       >  2:751120/1‑149 (MQ=255)
                      aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTCGTAGGTCGGATAAGG‑‑CGGTC‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGCCGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActc                                                       >  2:898081/1‑149 (MQ=14)
                      aGGCCTACGCCAGACAGCGCAAGAGCCTTATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActc                                                       >  2:1076299/1‑149 (MQ=14)
                        gCCTACGCCAGACGCCGCAATAGCCCGACTTAGCGTGATTTAATAGCTTGGATAAGG‑‑CGTTT‑‑AG‑GCAGTTCATTACATCAATGAATACTGCTACTCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                                                     <  1:700891/149‑1 (MQ=12)
                          cTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTca                                                   >  2:582112/1‑149 (MQ=255)
                           tACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGGCGCTTGCCGCGTCTTATCAGGCGTATCTTAACCGTTGGTTAATTTTTCAAACGCTTcac                                                  >  1:359013/1‑149 (MQ=17)
                              gCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc                                               >  2:1025397/1‑149 (MQ=255)
                              gCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc                                               >  2:1238780/1‑149 (MQ=255)
                                cAGTCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGCTCGGATAAGG‑‑CGTTT‑‑ATGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCg                                             <  2:618152/149‑1 (MQ=18)
                                   aCAGCGCAATGTCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGT‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTa                                          <  2:195902/149‑1 (MQ=18)
                                    catggcAATGGCCTGAGTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTAt                                         <  2:203916/145‑1 (MQ=18)
                                       cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaa                                      >  1:211823/1‑149 (MQ=255)
                                        gcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                     <  1:245273/149‑1 (MQ=255)
                                        gcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                     <  1:1305565/149‑1 (MQ=255)
                                          aaTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAcc                                   >  1:845053/1‑149 (MQ=255)
                                              gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTc                               <  1:495953/149‑1 (MQ=255)
                                              gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTc                               <  1:1003731/149‑1 (MQ=255)
                                               ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTcc                              >  2:354632/1‑149 (MQ=255)
                                               ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTcc                              >  2:109764/1‑149 (MQ=255)
                                                ttGACTTACCGTGATTCTGTAGTTCCGATAAGG‑‑CGGTT‑‑AGGCCGAACC‑CGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCa                             <  2:379513/148‑1 (MQ=21)
                                                cTGATTTGGCGTGATTTTGTAGGTCGGCTAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCa                             <  2:825493/149‑1 (MQ=21)
                                                cTGAGTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCa                             <  1:1151167/149‑1 (MQ=255)
                                                   aTTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATAc                          >  2:806069/1‑149 (MQ=255)
                                                    tttAGTCTGATTTCGTTGGGCGGCTAGGG‑‑CGTTTCCG‑‑CCGCATC‑CGACACCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGTCCTATCTTATCCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACa                         <  2:180319/149‑1 (MQ=21)
                                                     ttAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGAGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAg                        >  2:79295/1‑149 (MQ=255)
                                                       aGCGTGATTTTGGAGGTCGGATAAGG‑‑CGTTG‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACCCCCGGTATCAAACCCTTCCATACAGCt                      <  2:736430/149‑1 (MQ=21)
                                                        gCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCAGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTc                     >  2:963170/1‑149 (MQ=255)
                                                         cGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCa                    >  2:117936/1‑149 (MQ=255)
                                                          gTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAg                   <  1:542908/149‑1 (MQ=255)
                                                          gTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAg                   <  1:1274826/149‑1 (MQ=255)
                                                          gTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAg                   <  1:1186710/149‑1 (MQ=255)
                                                           tGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAgg                  <  1:1224546/149‑1 (MQ=255)
                                                              ttttGTAGGTCGGATCCGA‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGACGATTGCCGCGTCTTAGCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTc               >  2:532743/1‑149 (MQ=25)
                                                                     tgTCTGATAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGCCATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTAAACCatg        <  2:538679/148‑3 (MQ=25)
                                                                      gTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATCCAGCTCAGGTTCAACCAGTg       >  2:219123/1‑149 (MQ=255)
                                                                      gTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGTTTCACCAAGTg       >  1:64851/1‑149 (MQ=25)
                                                                         ggATAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAg    >  1:802535/1‑149 (MQ=255)
                                                                         ggATAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAg    >  1:1107501/1‑149 (MQ=255)
                                                                          gATAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCCGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCaac   >  1:661528/1‑147 (MQ=37)
                                                                           aTAAGG‑‑CGTTT‑‑ACGCCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTGTCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAgcg  >  1:812603/1‑149 (MQ=255)
                                                                                           ||                                                                                                                                        
GCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑A‑‑CCGCATC‑CGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCG  >  NC_000913_3_pae_tpiA/4296210‑4296431

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