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‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCG  >  NC_000913_3_pae_tpiA/4296210‑4296431
                                                                                           ||                                                                                                                                        
gCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAc                                                                             <  1:1099949/149‑1 (MQ=2)
ccGACGCTTAACGCGTCTTATCAGGCCTCCGCCAGCCAGCGCACTAGCCTGATTTAGCGTGATTTTGTGCATTGTATACCGCACGAGG‑‑G‑‑GCGCATCCGATATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAc                                                                             <  1:1038993/148‑1 (MQ=2)
    cGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGtt                                                                         <  2:499204/149‑1 (MQ=255)
       ttAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCTTTGGt                                                                      >  2:1102227/1‑149 (MQ=255)
          ccgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCGGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTaa                                                                   <  1:985287/148‑1 (MQ=255)
            gcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtt                                                                 >  2:1117244/1‑149 (MQ=255)
             cgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGGGTCTTATCAGGCATATCTTAACCGTTGGTTAAttt                                                                >  2:1353328/1‑149 (MQ=14)
              gTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ATGCCGCATCCGACATCAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTATCTTAGCCGTCGGATAAtttt                                                               >  2:558589/1‑149 (MQ=11)
              gTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ATGCCGCATCCGACATCAACGCCTGATGCGACGCTTCACGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtttt                                                               >  2:558609/1‑149 (MQ=11)
              gTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtttt                                                               >  2:1417205/1‑149 (MQ=255)
               tCTTATCAGGCCTACGCCAGACAGCGCAATCGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCATCTTATCAGGCCTATCTTACCGTTTGGTTAAtttat                                                              >  1:234798/1‑147 (MQ=11)
               tCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTCTCTTAACCGTTGGTTACttttt                                                              >  2:343118/1‑149 (MQ=11)
                 ttATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTAGTTAATTTTTCa                                                            >  2:102507/1‑149 (MQ=255)
                  tATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTcca                                                           >  2:260253/1‑147 (MQ=255)
                  tATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaa                                                           <  1:133881/149‑1 (MQ=255)
                  tATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaa                                                           >  2:217672/1‑149 (MQ=255)
                   aTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGGAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa                                                          <  1:382420/149‑1 (MQ=255)
                    tCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc                                                         >  1:753253/1‑149 (MQ=255)
                      aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ATGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActc                                                       <  1:558609/149‑1 (MQ=255)
                      aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActc                                                       <  1:558589/149‑1 (MQ=255)
                      aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGACCTATCTTAACCGTTGGTTAATTTTTCAAActc                                                       >  2:751120/1‑149 (MQ=255)
                      aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTCGTAGGTCGGATAAGG‑‑CGGTC‑‑ACGCCGCATCCGACATCAATGCCTGATGCGCCGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActc                                                       >  2:898081/1‑149 (MQ=14)
                      aGGCCTACGCCAGACAGCGCAAGAGCCTTATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActc                                                       >  2:1076299/1‑149 (MQ=14)
                        gCCTACGCCAGACGCCGCAATAGCCCGACTTAGCGTGATTTAATAGCTTGGATAAGG‑‑CGTTT‑‑AGGCAGTTCATTACATCAATGAATACTGCTACTCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                                                     <  1:700891/149‑1 (MQ=12)
                          cTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTca                                                   >  2:582112/1‑149 (MQ=255)
                           tACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGGCGCTTGCCGCGTCTTATCAGGCGTATCTTAACCGTTGGTTAATTTTTCAAACGCTTcac                                                  >  1:359013/1‑149 (MQ=17)
                              gCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc                                               >  2:1025397/1‑149 (MQ=255)
                              gCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc                                               >  2:1238780/1‑149 (MQ=255)
                                cAGTCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGCTCGGATAAGG‑‑CGTTT‑‑ATGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCg                                             <  2:618152/149‑1 (MQ=18)
                                   aCAGCGCAATGTCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGT‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTa                                          <  2:195902/149‑1 (MQ=18)
                                    catggcAATGGCCTGAGTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTAt                                         <  2:203916/145‑1 (MQ=18)
                                       cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaa                                      >  1:211823/1‑149 (MQ=255)
                                        gcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                     <  1:245273/149‑1 (MQ=255)
                                        gcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                     <  1:1305565/149‑1 (MQ=255)
                                          aaTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAcc                                   >  1:845053/1‑149 (MQ=255)
                                              gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTc                               <  1:495953/149‑1 (MQ=255)
                                              gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTc                               <  1:1003731/149‑1 (MQ=255)
                                               ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTcc                              >  2:354632/1‑149 (MQ=255)
                                               ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTcc                              >  2:109764/1‑149 (MQ=255)
                                                ttGACTTACCGTGATTCTGTAGTTCCGATAAGG‑‑CGGTT‑‑AGGCCGAACCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCa                             <  2:379513/148‑1 (MQ=21)
                                                cTGATTTGGCGTGATTTTGTAGGTCGGCTAAGG‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCa                             <  2:825493/149‑1 (MQ=21)
                                                cTGAGTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCa                             <  1:1151167/149‑1 (MQ=255)
                                                   aTTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATAc                          >  2:806069/1‑149 (MQ=255)
                                                    tttAGTCTGATTTCGTTGGGCGGCTAGGG‑‑CGTTTCCG‑‑CCGCATCCGACACCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGTCCTATCTTATCCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACa                         <  2:180319/149‑1 (MQ=21)
                                                     ttAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGAGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAg                        >  2:79295/1‑149 (MQ=255)
                                                       aGCGTGATTTTGGAGGTCGGATAAGG‑‑CGTTG‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACCCCCGGTATCAAACCCTTCCATACAGCt                      <  2:736430/149‑1 (MQ=21)
                                                        gCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCAGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTc                     >  2:963170/1‑149 (MQ=255)
                                                         cGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCa                    >  2:117936/1‑149 (MQ=255)
                                                          gTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAg                   <  1:542908/149‑1 (MQ=255)
                                                          gTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAg                   <  1:1274826/149‑1 (MQ=255)
                                                          gTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAg                   <  1:1186710/149‑1 (MQ=255)
                                                           tGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAgg                  <  1:1224546/149‑1 (MQ=255)
                                                              ttttGTAGGTCGGATCCGA‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGATTGCCGCGTCTTAGCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTc               >  2:532743/1‑149 (MQ=25)
                                                                     tgTCTGATAAGG‑‑CGTTT‑‑ACGCCGCATCCGCCATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTAAACCatg        <  2:538679/148‑3 (MQ=25)
                                                                      gTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATCCAGCTCAGGTTCAACCAGTg       >  2:219123/1‑149 (MQ=255)
                                                                      gTCGGATAAGG‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGTTTCACCAAGTg       >  1:64851/1‑149 (MQ=25)
                                                                         ggATAAGG‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAg    >  1:802535/1‑149 (MQ=255)
                                                                         ggATAAGG‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAg    >  1:1107501/1‑149 (MQ=255)
                                                                          gATAAGG‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCCGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCaac   >  1:661528/1‑147 (MQ=37)
                                                                           aTAAGG‑‑CGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTGTCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAgcg  >  1:812603/1‑149 (MQ=255)
                                                                                           ||                                                                                                                                        
GCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGG‑‑CGTTT‑‑A‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCG  >  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: