Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_tpiA 3,215,447 C→A G9G (GGG→GGT)  mug ← G/U mismatch‑specific DNA glycosylase; xanthine DNA glycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA3,215,4470CA98.5% 230.9 / ‑5.1 66G9G (GGG→GGT) mugG/U mismatch‑specific DNA glycosylase; xanthine DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (29/36);  total (30/36)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

GTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCGGCCTGATATATCACCTTCCAGAAGCGATTTGCCGGATGAGCAAAGGGAAAACCAGTCCCGGCGGATGAAAGCCCAGGGTTGATACCGCAAAACACGACCCGTAACCCTGGAGCCAAAATATCCTCAACCATATTTACTCCTTTCTAACATCACCAGGAAAAGTATAACGGATTGAAAATACATTGTTTATAAAAACAGCAGGCGCGCGGTAATGGCTGGATTGCGACACGGAGTTACTTTATA  >  NC_000913_3_pae_tpiA/3215305‑3215585
                                                                                                                                              |                                                                                                                                          
gtgCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCGGCCTGATATATCACCTTCCAGAAGCGATTTGCCGGATGAGCAAAGGGAAAACCAGTCCCGGCGGATGAAAGCCCAGGGTTGATACCGCAAAACACGACCCGTAAACCTGGAg                                                                                                                                     >  1:742335/1‑150 (MQ=255)
         gCGGCTTCAACTGACGGTCGGTAAACCCGGCCTGATATATCACCTTCCAGAAGCGATTTGCCGGATGAGCAAAGGGAAAACCAGTCCCGGCGGATGAAAGCCCAGGGTTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAatat                                                                                                                            >  2:355043/1‑150 (MQ=255)
           ggCTTCAACTGACGGTCGGTAAACCCGGCCTGATATATCACCTTCCAGAAGCGATTTGCCGGATGAGCAAAGGGAAAACCAGTCCCGGCGGATGAAAGCCCAGGGTTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATcc                                                                                                                          <  4:490361/150‑1 (MQ=255)
           ggCTTCAACTGACGGTCGGTAAACCCGGCCTGATATATCACCTTCCAGAAGCGATTTGCCGGATGAGCAAAGGGAAAACCAGTCCCGGCGGATGAAAGCCCAGGGTTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATcc                                                                                                                          <  2:667223/150‑1 (MQ=255)
                 aaCTGACGGTCGGTAAACCCGGCCTGATATATCACCTTCCAGAAGCGATTTGCCGGATGAGCAAAGGGAAAACCAGTCCCGGCGGATGAAAGCCCAGGGTTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAAcc                                                                                                                    >  4:359249/1‑150 (MQ=255)
                 aaCTGACGGTCGGTAAACCCGGCCTGATATATCACCTTCCAGAAGCGATTTGCCGGATGAGCAAAGGGAAAACCAGTCCCGGCGGATGAAAGCCCAGGGTTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAAcc                                                                                                                    >  1:437590/1‑150 (MQ=255)
                      aCGGTCGGTAAACCCGGCCTGATATATCACCTTCCAGAAGCGATTTGCCGGATGAGCAAAGGGAAAACCAGTCCCGGCGGATGAAAGCCCAGGGTTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATAtt                                                                                                               <  2:34635/150‑1 (MQ=255)
                        ggtcggtAAACCCGGCCTGATATATCACCTTCCAGAAGCGATTTGCCGGATGAGCAAAGGGAAAACCAGTCCCGGCGGATGAAAGCCCAGGGTTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTa                                                                                                             <  1:355043/150‑1 (MQ=255)
                         gtcggtAAACCCGGCCTGATATATCACCTTCCAGAAGCGATTTGCCGGATGAGCAAAGGGAAAACCAGTCCCGGCGGATGAAAGCCCAGGGTTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTAc                                                                                                            >  2:4992/1‑150 (MQ=255)
                         gtcggtAAACCCGGCCTGATATATCACCTTCCAGAAGCGATTTGCCGGATGAGCAAAGGGAAAACCAGTCCCGGCGGATGAAAGCCCAGGGTTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTAc                                                                                                            <  3:11486/150‑1 (MQ=255)
                         gtcggtAAACCCGGCCTGATATATCACCTTCCAGAAGCGATTTGCCGGATGAGCAAAGGGAAAACCAGTCCCGGCGGATGAAAGCCCAGGGTTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTAc                                                                                                            >  2:161148/1‑150 (MQ=255)
                            ggtAAACCCGGCCTGATATATCACCTTCCAGAAGCGATTTGCCGGATGAGCAAAGGGAAAACCAGTCCCGGCGGATGAAAGCCCAGGGTTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACTcc                                                                                                         >  1:760880/1‑150 (MQ=255)
                             gtAAACCCGGCCTGATATATCACCTTCCAGAAGCGATTTGCCGGATGAGCAAAGGGAAAACCAGTCCCGGCGGATGAAAGCCCAGGGTTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACTCCt                                                                                                        >  1:728817/1‑150 (MQ=255)
                             gtAAACCCGGCCTGATATATCACCTTCCAGAAGCGATTTGCCGGATGAGCAAAGGGAAAACCAGTCCCGGCGGATGAAAGCCCAGGGTTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACTCCt                                                                                                        >  1:818274/1‑150 (MQ=255)
                             gtAAACCCGGCCTGATATATCACCTTCCAGAAGCGATTTGCCGGATGAGCAAAGGGAAAACCAGTCCCGGCGGATGAAAGCCCAGGGTTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACTCCt                                                                                                        <  1:287544/150‑1 (MQ=255)
                                 gctaggacTGTTATAGCACCTTCCAGGAGCGATTTTGCAGATAAGCTAAGGGGAAACCAATCCCTGCGGATTAAAGACCAGGTTTTATACCGTAAAACACGTCGCGTAAACCTGAAGCCAAAATTTCCTCAACCATATTTACTCCTTTCt                                                                                                    <  2:635705/143‑1 (MQ=255)
                                   ccGGCCTGATATATCACCTTCCAGAAGCGATTTGCCGGATGAGCAAAGGGAAAACCAGTCCCGGCGGATGAAAGCCCAGGGTTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACTCCTTTCTaa                                                                                                  >  1:292206/1‑150 (MQ=255)
                                   ccGGCCTGATATATCACCTTCCAGAAGCGATTTGCCGGATGAGCAAAGGGAAAACCAGTCCCGGCGGATGAAAGCCCAGGGTTGACACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACTCCTTTCTaa                                                                                                  >  2:589055/1‑150 (MQ=255)
                                        cTGATATATCACCTTCCAGAAGCGATTTGCCGGATGAGCAAAGGGAAAACCAGTCCCGGCGGATGAAAGCCCAGGGTTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACTCCTTTCTAACATCa                                                                                             >  1:398975/1‑150 (MQ=255)
                                           atatatCACCTTCCAGAAGCGATTTGCCGGATGAGCAAAGGGAAAACCAGTCCCGGCGGATGAAAGCCCAGGGTTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACTCCTTTCTAACATCACCa                                                                                          <  2:247389/150‑1 (MQ=255)
                                           atatatCACCTTCCAGAAGCGATTTGCCGGATGAGCAAAGGGAAAACCAGTCCCGGCGGATGAAAGCCCAGGGTTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACTCCTTTCTAACATCACCa                                                                                          <  4:434333/150‑1 (MQ=255)
                                                  aCCTTCCAGAAGCGATTTGCCGGATGACCAAAGGGAAAACCAGTCCCGGCGGATGAAAGCCCAGGGGTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACTCCTTTCTAACATCACCAGGAAAAg                                                                                   >  2:14302/1‑150 (MQ=255)
                                                       ccAGAAGCGATTTGCCGGATGAGCAAAGGGAAAACCAGTCCCGGCGGATGAAAGCCCAGGGTTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACTCCTTTCTAACATCACCAGGAAAAGTATaa                                                                              >  1:260558/1‑150 (MQ=255)
                                                         aGAAGCGATTTGCCGGATGAGCAAAGGGAAAACCAGTCCCGGCGGATGAAAGCCCAGGGTTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACTCCTTTCTAACATCACCAGGAAAAGTATAACg                                                                            <  2:818274/150‑1 (MQ=255)
                                                         aGAAGCGATTTGCCGGATGAGCAAAGGGAAAACCAGTCCCGGCGGATGAAAGCCCAGGGTTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACTCCTTTCTAACATCACCAGGAAAAGTATAACg                                                                            <  1:161148/150‑1 (MQ=255)
                                                          gAAGCGATTTGCCGGATGAGCAAAGGGAAAACCAGTCCCGGCGGATGAAAGCCCAGGGTTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACTCCTTTCTAACATCACCAGGAAAAGTATAACgg                                                                           >  4:302406/1‑150 (MQ=255)
                                                              cGATTTGCCGGATGAGCAAAGGGAAAACCAGTCCCGGCGGATGAAAGCCCAGGGTTGATACCGCAAAACACGACCCGAAAACCAGGAGCCAAAATATCCTCAACCATATTTACTCCTTTCTAACATCACCAGGAAAAGTATAACGGATTg                                                                       >  3:49064/1‑150 (MQ=255)
                                                                 tttGCCGGATGAGCAAAGGGAAAACCAGTCCCGGCGGATGAAAGCCCAGGGTTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATc                                                                                                                           >  4:40937/1‑95 (MQ=255)
                                                                 tttGCCGGATGAGCAAAGGGAAAACCAGTCCCGGCGGATGAAAGCCCAGGGTTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATc                                                                                                                           <  3:40937/95‑1 (MQ=255)
                                                                   taaccactGAGCAAAGGGAAAACCAGTCCCGGCGGATGAAAGCCCAGGGTTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACTCCTTTCTAACATCACCAGGAAAAGTATAACGGATTGAAAAt                                                                  >  4:76822/8‑150 (MQ=255)
                                                                     ccGGATGAGCAAAGGGAAAACCAGTCCCGGCGGATGAAAGCCCAGGGTTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACTCCTTTCTAACATCACCAGGAAAAGTATAACGGATTGAAAATAc                                                                <  2:141654/150‑1 (MQ=255)
                                                                       ggATGAGCAAAGGGAAAACCAGTCCCGGCGGATGAAAGCCCAGGGTTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACTCCTTTCTAACATCACCAGGAAAAGTATAACGGATTGAAAATACAt                                                              <  2:128958/150‑1 (MQ=255)
                                                                        gATGAGCAAAGGGAAAACCAGTCCCGGCGGATGAAAGCCCAGGGTTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACTCCTTTCTAACATCTCCAGGAAAAGTATAACGGATTGAAAATACAtt                                                             <  4:335540/150‑1 (MQ=255)
                                                                         aTGAGCAAAGGGAAAACCAGTCCCGGCGGATGAAAGCCCAGGGTTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACTCCTTTCTAACATCACCAGGAAAAGTATAACGGATTGAAAATACATTg                                                            <  3:513893/150‑1 (MQ=255)
                                                                           gAGCAAAGGGAAAACCAGTCCCGGCGGATGAAAGCCCAGGGTTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACTCCTTTCTAACATCACCAGGAAAAGTATAACGGATTGAAAATACATTGtt                                                          <  4:149317/150‑1 (MQ=255)
                                                                            aGCAAAGGGAAAACCAGTCCCGGCGGATGAAAGCCCAGGGTTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACTCCTTTCTAACATCACCAGGAAAAGTATAACGGATTGAAAATACATTGttt                                                         <  4:378237/150‑1 (MQ=255)
                                                                                 aGGGAAAACCAGTCCCGGCGGATGAAAGCCCAGGGTTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACTCCTTTCTAACATCACCAGGAAAAGTATAACGGATTGAAAATACATTGTTTATaaa                                                    <  2:312879/150‑1 (MQ=255)
                                                                                   ggAAAACCAGTCCCGGCGGATGAAAGCCCAGGGTTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACTCCTTTCTAACATCACCAGGAAAAGTATAACGGATTGAAAATACATTGTTTATaaaaa                                                  >  1:804589/1‑150 (MQ=255)
                                                                                    gAAAACCAGTCCCGGCGGATGAAAGCCCAGGGTTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACTCCTTTCTAACATCACCAGGAAAAGTATAACGGATTGAAAATACATTGTTTATAAAAAc                                                 <  2:458329/150‑1 (MQ=255)
                                                                                               ccGGCGGATGAAAGCCCAGGGTTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACTCCTTTCTAACATCACCAGGAAAAGTATAACGGATTGAAAATACATTGTTTATAAAAAcag                                               <  2:152899/141‑1 (MQ=255)
                                                                                               ccGGCGGATGAAAGCCCAGGGTTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACTCCTTTCTAACATCACCAGGAAAAGTATAACGGATTGAAAATACATTGTTTATAAAAAcag                                               >  1:152899/1‑141 (MQ=255)
                                                                                                cggcggATGAAAGCCCAGGGTTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACTCCTTTCTAACATCACCAGGAAAAGTATAACGGATTGAAAATACATTGTTTATAAAAACAGCAGGcgcgcg                                     >  1:150090/1‑150 (MQ=255)
                                                                                                       tGAAAGCCCAGGGTTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACTCCTTTCTAACATCACCAGGAAAAGTATAACGGATTGAAAATACATTGTTTATAAAAACAGCAGGCGCGCGGTAATgg                              <  3:302406/150‑1 (MQ=255)
                                                                                                        gAAAGCCCAGGGTTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACt                                                                                                           <  2:590338/72‑1 (MQ=255)
                                                                                                        gAAAGCCCAGGGTTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACt                                                                                                           >  1:590338/1‑72 (MQ=255)
                                                                                                        gAAAGCCCAGGGTTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACTCCTTTCTAACATCACCAGGAAAAGTATAACGGATTGAAAATACATTGTTtat                                                       >  4:242246/1‑124 (MQ=255)
                                                                                                        gAAAGCCCAGGGTTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACTCCTTTCTAACATCACCAGGAAAAGTATAACGGATTGAAAATACATTGTTtat                                                       <  3:242246/124‑1 (MQ=255)
                                                                                                          aaGCCCAGGGTTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACTCCTTTCTAACATCACCAGGAAAAGTATAACGGATTGAAAATACATTGTTTATAAAAACAGCAGGCGCGCGGTAATGGCTg                           <  3:359249/150‑1 (MQ=255)
                                                                                                          aaGCCCAGGGTTGATACCGCAAAACACAACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACTCCTTTCTAACATCAc                                                                                            >  1:432589/1‑85 (MQ=255)
                                                                                                          aaGCCCAGGGTTGATACCGCAAAACACAACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACTCCTTTCTAACATCAc                                                                                            <  2:432589/85‑1 (MQ=255)
                                                                                                              ccAGGGTTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACTCCTTTCTAACATCACCAGGAAAAGTATAACGGATTGAAAATACATTGTTTATAAAAACAGCAGGCGCGCGGTAATGGCTGGAtt                       >  4:390932/1‑150 (MQ=255)
                                                                                                              ccAGGGTTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACTCCTTTCTAACATCACCAGGAAAAGTATAACGGATTGAAAATACATTGTTTATAAAAACAGCAGGCGCGCGGTAATGGCTGGAtt                       <  1:281467/150‑1 (MQ=255)
                                                                                                                aGGGTTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACTCCTTTCTAACATCACCAGGAAAAGTATAACGGATTGAAAATACATTGTTTATAAAAACAGCAGGCGCGCGGTAATGGCTGGATTGc                     <  3:134206/150‑1 (MQ=255)
                                                                                                                 gggTTGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACTCCTTTCTAACATCAACAGGAAAAGTATAATGGATTGAAAATACATTGTTTATAAAAACAGCAGGCGCGCGGTAATGGCTGGATTGCg                    <  2:437590/150‑1 (MQ=255)
                                                                                                                    ttGATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACTCCTTTCTAACATCACCAGGAAAAGTATAACGGATTGAAAATACATTGTTTATAAAAACAGCAGGCGCGCGGTAATGGCTGGATTGCGaca                 <  3:236105/150‑1 (MQ=255)
                                                                                                                      gATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACTCCTTTCTAACATCAcc                                                                                           <  2:100930/74‑1 (MQ=255)
                                                                                                                      gATACCGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACTCCTTTCTAACATCAcc                                                                                           >  1:100930/1‑74 (MQ=255)
                                                                                                                        tACCGCAAAACACGACCCGTAAACCTGAAGCCAAAATATCCTCAACCATATTTACTCCTTTCTAACATCACCAGGAAAAGTATAACGGATTGAAAATACATTGTTTATAAAAACAGCAGGCGCGCGGTAATGGCTGGATTGCGACACGGa             <  4:49064/150‑1 (MQ=255)
                                                                                                                           cGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACTCCTTTCTAACATCACCAGGAAAAGTATAACGGATTGAAAATACATTGTTTATAAAAACAGCAGGCGCGCGGTAATGGCTGGATTGCGACACGGAGtt          >  3:167452/1‑150 (MQ=255)
                                                                                                                           cGCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACTCCTTTCTAACATCACCAGGAAAAGTATAACGGATTGAAAATACATTGTTTATAAAAACAGCAGGCGCGCGGTAATGGCTGGATTGCGACACGGAGtt          >  2:716378/1‑150 (MQ=255)
                                                                                                                            gCAAAACACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACTCCTTTCTAACATCACCAGGAAAAGTATAACGGATTGAAAATACATTGTTTATAAAAACAGCAGGCGCGCGGTAATGGCTGGATTGCGACACGGAGTTa         <  2:804589/150‑1 (MQ=255)
                                                                                                                              aaaaCACGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACTCCTTTCTAACATCACCAGGAAAAGTATAACGGATTGAAAATACATTGTTTATAAAAACAGCAGGCGCGCGGTAATGGCTGGATTGCGACACGGAGTTACt       <  2:150090/150‑1 (MQ=255)
                                                                                                                                 acacGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACTCCTTTCTAACATCACCAGGAAAAGTATAACGGATTGAAAATACATTGTTTATAAAAACAGCAGGCGCGCGGTAATGGTTGGATTGCGACACGGAGTTACTTta    <  4:299440/150‑1 (MQ=255)
                                                                                                                                 acacGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACTCCTTTCTAACATCACCAGGAAAAGTATAACGGATTGAAAATACATTGTTTATAAAAACAGCAGGCGCGCGGTAATGGCTGGATTGCGACACGGAGTTACTTta    <  4:279308/150‑1 (MQ=255)
                                                                                                                                 acacGACCCGTAAACCTGGAGCCAAAATATCCTCAACCATATTTACTCCTTTCTAACATCACCAGGAAAAGTATAACGGATTGAAAATACATTGTTTATAAAAACAGCAGGCGCGCGGTAATGGCTGGATTGCGACACGGAGTTACTTta    >  2:61223/1‑150 (MQ=255)
                                                                                                                                   acGACCCGTAACCCTGGAGCCAAAATATCCTCAACCATATTTACTCCTTTCTAACATCACCAGGAAAAGTATAACGGATTGAAAATACATTGATTATAAAAACAGCAGGCGCGCGGTAATGGCTGGATTGCGACACGGAGTTACTTtata  >  1:810747/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                          
GTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCGGCCTGATATATCACCTTCCAGAAGCGATTTGCCGGATGAGCAAAGGGAAAACCAGTCCCGGCGGATGAAAGCCCAGGGTTGATACCGCAAAACACGACCCGTAACCCTGGAGCCAAAATATCCTCAACCATATTTACTCCTTTCTAACATCACCAGGAAAAGTATAACGGATTGAAAATACATTGTTTATAAAAACAGCAGGCGCGCGGTAATGGCTGGATTGCGACACGGAGTTACTTTATA  >  NC_000913_3_pae_tpiA/3215305‑3215585

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