Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,117,263 C→T 6.1% N136N (AAC→AAT)  iorA‑II → isoquinoline 1‑oxidoreductase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,117,2630CT 6.1% 196.9 / 2.4 65N136N (AAC→AAT) iorA‑IIisoquinoline 1‑oxidoreductase subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (33/28);  new base T (2/2);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

GTCGAGCACCAGGTGGCTCAATGCGGCTATTGCCAGTCTGGCCAGGTCATGGCAGCCACGGCGCTGCTCAAGCACACGCCCAAACCCAGCGATGCGCAGATCGATGCAGCCATGGTCAACCTGTGCCGCTGCGGCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAGGGGGAGCAGGCCTGATGAACACGCCAGTCGATACACCCGCCGAACTGCTCAGGCTGCAGCGGGGCGAAACGGTCAACCTGTCACGTCGGCGCTTTCTTGCCGGCACCGCAGTCG  >  NC_002947/4117120‑4117410
                                                                                                                                               |                                                                                                                                                   
gTCGAGCACCAGGTGGCTCAATGCGGCTATTGCCAGTCTGGCCAGGTCATGGCAGCCACGGCGCTGCTCAAGCACACGCCCAAACCCAGCGATGCGCAGATCGATGCAGCCATGGTCAACCTGTGCCGCTGCGGCACCTACAACGccat                                                                                                                                                <  2:43642/149‑1 (MQ=255)
 tCGAGCACCAGGTGGCTCAATGCGGCTATTGCCAGTCTGGCCAGGTCATGGCAGCCACGGCGCTGCTCAAGCACACGCCCAAACCCAGCGATGCGCAGATCGATGCAGCCATGGTCAACCTGTGCCGCTGCGGCACCTACAACGccatc                                                                                                                                               >  5:41247/1‑149 (MQ=255)
  cGAGCACCAGGTGGCTCAATGCGGCTATTGCCAGTCTGGCCAGGTCATGGCAGCCACGGCGCTGCTCAAGCACACGCCCAAACCCAGCGATGCGCAGATCGATGCAGCCATGGTCAACCTGTGCCGCTGCGGCACCTACAACGccatcc                                                                                                                                              <  7:117498/149‑1 (MQ=255)
         cAGGTGGCTCAATGCGGCTATTGCCAGTCTGGCCAGGTCATGGCAGCCACGGCGCTGCTCAAGCACACGCCCAAACCCAGCGATGCGCAGATCGATGCAGCCATGGTCAACCTGTGCCGCTGCGGCACCTACAACGCCATCCATGCCGc                                                                                                                                       >  4:329506/1‑149 (MQ=255)
            gTGGCTCAATGCGGCTATTGCCAGTCTGGCCAGGTCATGGCAGCCACGGCGCTGCTCAAGCACACGCCCAAACCCAGCGATGCGCAGATCGATGCAGCCATGGTCAACCTGTGCCGCTGCGGCACCTACAACGCCATCCATGCCGCtgt                                                                                                                                    >  7:250313/1‑149 (MQ=255)
                cTCAATGCGGCTATTGCCAGTCTGGCCAGGTCATGGCAGCCACGGCGCTGCTCAAGCACACGCCCAAACCCAGCGATGCGCAGATCGATGCAGCCATGGTCAACCTGTGCCGCTGCGGCACCTACAACGCCATCCATGCCGCTGTGCAt                                                                                                                                <  4:358382/149‑1 (MQ=255)
                 tCAATGCGGCTATTGCCAGTCTGGCCAGGTCATGGCAGCCACGGCGCTGCTCAAGCACACGCCCAAACCCAGCGATGCGCAGATCGATGCAGCCATGGTCAACCTGTGCCGCTGCGGCACCTACAAc                                                                                                                                                     <  5:237315/127‑1 (MQ=255)
                 tCAATGCGGCTATTGCCAGTCTGGCCAGGTCATGGCAGCCACGGCGCTGCTCAAGCACACGCCCAAACCCAGCGATGCGCAGATCGATGCAGCCATGGTCAACCTGTGCCGCTGCGGCACCTACAAc                                                                                                                                                     >  6:237315/1‑127 (MQ=255)
                             ttGCCAGTCTGGCCAGGTCATGGCAGCCACGGCGCTGCTCAAGCACACGCCCAATCCCAGCGATGCGCAGATCGATGCAGCCATGGTCAACCTGTGCCGCTGCGGCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAgg                                                                                                                   <  1:264369/149‑1 (MQ=255)
                              tGCCAGTCTGGCCAGGTCATGGCAGCCACGGCGCTGCTCAAGCACACGCCCAAACCCAGCGATGCGCAGATCGATGCAGCCATGGTCAACCTGTGCCGCTGCGGCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAgg                                                                                                                   <  6:41247/148‑1 (MQ=255)
                                ccAGTCTGGCCAGGTCATGGCAGCCACGGCGCTGCTCAAGCACACGCCCAAACCCAGCGATGCGCAGATCGATGCAGCCATGGTCAACCTGTGCCGCTGCGGCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAgggg                                                                                                                 >  8:235213/1‑148 (MQ=255)
                                      tGGCCAGGTCATGGCAGCCACGGCGCTGCTCAAGCACACGCCCAAACCCAGCGATGCGCAGATCGATGCAGCCATGGTCAACCTGTGCCCCTCCGGCACCTACAACGCCATCCATCCCGCAGTGCATGACCTGGCCCAGGGGGAGCAgg                                                                                                          >  7:333794/1‑149 (MQ=255)
                                        gCCAGGTCATGGCAGCCACGGCGCTGCTCAAGCACACGCCCAAACCCAGCGATGCGCAGATCGATGCAGCCATGGTCAACCTGTGCCGCTGCGGCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAGGGGGAGCAGGcc                                                                                                        >  7:31820/1‑149 (MQ=255)
                                         ccAGGTCATGGCAGCCACGGCGCTGCTCAAGCACACGCCCAAACCCAGCGATGCGCAGATCGATGCAGCCATGGTCAACCTGTGCCGCTGCGGCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAGGGGGAGCAGGCCt                                                                                                       <  7:235213/149‑1 (MQ=255)
                                           aGGTCATGGCAGCCACGGCGCTGCTCAAGCACACGCCCAAACCCAGCGATGCGCAGATCGATGCAGCCATGGTCAACCTGTGCCGCTGCGGCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAGGGGGAGCAGGCCtga                                                                                                     <  6:192941/149‑1 (MQ=255)
                                            ggTCATGGCAGCCACGGCGCTGCTCAAGCACACGCCCAAACCCAGCGATGCGCAGATCGATGCAGCCATGGTCAACCTGTGCCGCTGCGGCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAGGGGGAGCAGGCCtgat                                                                                                    <  8:250313/149‑1 (MQ=255)
                                             gTCATGGCAGCCACGGCGCTGCTCAAGCACACGCCCAAACCCAGCGATGCGCAGATCGATGCAGCCATGGTCAACCTGTGCCGCTGCGGCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAGGGGGAGCAGGCCtgat                                                                                                    >  8:70621/1‑148 (MQ=255)
                                             gTCATGGCAGCCACGGCGCTGCTCAAGCACACGCCCAAACCCAGCGATGCGCAGATCGATGCAGCCATGGTCAACCTGTGCCGCTGCGGCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAGGGGGAGCAGGCCtgat                                                                                                    >  5:90625/1‑148 (MQ=255)
                                                       ccACGGCGCTGCTCAAGCACACGCCCAAACCCAGCGATGCGCAGATCGATGCAGCCATGGTCAACCTGTGCCGCTGCGGCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAGGGGGAGCAGGCCtgatga                                                                                                  <  3:394472/140‑1 (MQ=255)
                                                       ccACGGCGCTGCTCAAGCACACGCCCAAACCCAGCGATGCGCAGATCGATGCAGCCATGGTCAACCTGTGCCGCTGCGGCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAGGGGGAGCAGGCCtgatga                                                                                                  >  4:394472/1‑140 (MQ=255)
                                                                 gctCAAGCACACGCCCAAACCCAGCGATGCGCAGATCGATGCAGCCATGGTCAACCTGTGCCGCTGCGGCACCTACAATGccatccat                                                                                                                                            <  7:83933/88‑1 (MQ=255)
                                                                 gctCAAGCACACGCCCAAACCCAGCGATGCGCAGATCGATGCAGCCATGGTCAACCTGTGCCGCTGCGGCACCTACAATGccatccat                                                                                                                                            <  7:89393/88‑1 (MQ=255)
                                                                 gctCAAGCACACGCCCAAACCCAGCGATGCGCAGATCGATGCAGCCATGGTCAACCTGTGCCGCTGCGGCACCTACAATGccatccat                                                                                                                                            >  8:83933/1‑88 (MQ=255)
                                                                 gctCAAGCACACGCCCAAACCCAGCGATGCGCAGATCGATGCAGCCATGGTCAACCTGTGCCGCTGCGGCACCTACAATGccatccat                                                                                                                                            >  8:89393/1‑88 (MQ=255)
                                                                        cacaCGCCCAAACCCAGCGATGCGCAGATCGATGCAGCCATGGTCAACCTGTGCCGCTGCGGCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAGGGGGAGCAGGCCTGATGAACACGCCAGTCGATACACCCGCCGaa                                                                        <  6:90625/149‑1 (MQ=255)
                                                                         acacGCCCAAACCCAGCGATGCGCAGATCGATGCAGCCATGGTCAACCTGTGCCGCTGCGGCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAGGGGGAGCAGGCCTGATGAACACGCCAGTCGATACACCCGCCGAAc                                                                       <  8:31820/149‑1 (MQ=255)
                                                                         acacGCCCAAACCCAGCGATGCGCAGATCGATGCAGCCATGGTCAACCTGTCCCGCTGCGGCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAGGGGGAGCAGGCCTGATGAACACGCCAGTCGATACACCCGCCGaa                                                                        >  6:281577/1‑148 (MQ=255)
                                                                           acGCCCAAACCCAGCGATGCGCAGATCGATGCAGCCATGGTCAACCTGTGCCGCTGCGGCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAGGGGGAGCAGGCCTGATGAACACGCCAGTCGATACACCCGCCGAACTg                                                                     >  6:378273/1‑149 (MQ=255)
                                                                             gCCCAAACCCAGCGATGCGCAGATCGATGCAGCCATGGTCAACCTGTGCCGCTGCGGCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAGGGGGAGCAGGCCTGATGAACACGCCAGTCGATACACCCGCCGAACTGCt                                                                   >  3:185321/1‑149 (MQ=255)
                                                                                cAAACCCAGCGATGCGCAGATCGATGCAGCCATGGTCAACCTGTGCCGCTGCGGCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAGGGGGAGCAGGCCTGATGAACACGCCAGTCGATACACCCGCCGAACTGCTCAg                                                                <  5:378273/149‑1 (MQ=255)
                                                                                  aaCCCAGCGATGCGCAGATCGATGCAGCCATGGTCAACCTGTGCCGCTGCGGCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAggggg                                                                                                                <  1:38351/99‑1 (MQ=255)
                                                                                  aaCCCAGCGATGCGCAGATCGATGCAGCCATGGTCAACCTGTGCCGCTGCGGCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAggggg                                                                                                                >  2:38351/1‑99 (MQ=255)
                                                                                     ccAGCGATGCGCAGATCGATGCAGCCATGGTCAACCTGTGCCGCTGCGGCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAGGGGGAGCTGGCCTGATGAACACGCCAGTCGATACACCCGCCGAACTGCTCAGGCTg                                                            >  8:122388/1‑148 (MQ=255)
                                                                                      cAGCGATGCGCAGATCGATGCAGCCATGGTCAACCTGTGCCGCTGCGGCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAGGGGGAGCAGGCCTGATGAACACGCCAGTCGATACACCCGCCGAACTGCTCAGGCTGCa                                                          <  2:81268/149‑1 (MQ=255)
                                                                                          gATGCGAAGATCGATGCAGCCATGGTCTACCTGTGCCGCTGCGTCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAGGGGGAGCAGGCCTGATGAACACGCCAGTCGATACACCCGCCGAACTGCTCAGGCTGCAGCg                                                       <  8:333794/148‑1 (MQ=255)
                                                                                           aTGCGCAGATCGATGCAGCCATGGTCAACCTGTGCCGCTGCGGCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAGGGGGAGCAGGCCTGATGAACACGCCAGTCGATACACCCGCCGAACTGCTCAGGCTGCAGCggg                                                     <  3:174508/149‑1 (MQ=255)
                                                                                             gcgcAGATCGATGCAGCCATGGTCAACCTGTGCCGCTGCGGCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAGGGGGAGCAGGCCTGATGAACACGCCAGTCGATACACCCGCCGAACTGCTCAGGCTGCAGCGGGGc                                                   >  5:329000/1‑149 (MQ=255)
                                                                                               gcAGATCGATGCAGCCATGGTCAACCTGTGCCGCTGCGGCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAGGGGGAGCAGGCCTGATGAACACGCCAGTc                                                                                       >  4:395153/1‑111 (MQ=255)
                                                                                               gcAGATCGATGCAGCCATGGTCAACCTGTGCCGCTGCGGCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAGGGGGAGCAGGCCTGATGAACACGCCAGTc                                                                                       <  3:395153/111‑1 (MQ=255)
                                                                                               gcAGATCGATGCAGCCATGGTCAACCTGTGCCGCTGCGGCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAGGGGGAGCAGGCCTGATGAACACGCCAGTCGATACACCCGCCGAACTGCTCAGGCTGCAGCGGGGCg                                                  >  5:322522/1‑148 (MQ=255)
                                                                                                   aTCGATGCAGCCATGGTCAACCTGTGCCGCTGCGGCACCTACAACGCCATCCATGCCGCTGTGCATTACCTGGCCCAGGGGGAGCAGGCCTGATGAACACGCCAGTCGATACACCCGCCGAACTGCTCAGGCTGCAGCGGGGCGAAACg                                             <  3:329506/149‑1 (MQ=255)
                                                                                                      gATGCAGCCATGGTCAACCTGTGCCGCTGCGGCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAGGGGGAGCAGGCCTGATGAACACGCCAGTCGATACACCCGCCGAACTGCTCAGGCTGCAGCGGGGCGAAACGGTc                                          <  5:281577/149‑1 (MQ=255)
                                                                                                           aGCCATGGTCAACCTGTGCCGCTGCGGCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAGGGGGAGCAGGCCTGATGAACACGCCAGTCGATACACCCGCCGAACTGCTCAGGCTGCAGCGGGGCGAAACGGTCAACCt                                     <  4:185321/149‑1 (MQ=255)
                                                                                                             ccATGGTCAACCTGTGCCGCTGCGGCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAGGGGGAGCAGGCCTGATGAACACGCCAGTCGATACACCCGCCGAACTGCTCAGGCTGCAGCgggg                                                    >  8:348472/1‑132 (MQ=255)
                                                                                                             ccATGGTCAACCTGTGCCGCTGCGGCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAGGGGGAGCAGGCCTGATGAACACGCCAGTCGATACACCCGCCGAACTGCTCAGGCTGCAGCgggg                                                    <  7:348472/132‑1 (MQ=255)
                                                                                                             ccATGGTCAACCTGTGCCGCTGCGGCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAGGGGGAGCAGGCCTGATGAACACGCCAGTCGATACACCCGCCGAACTGCTCAGGCTGCAGCGGGGCGAAACGGTCAACCTg                                    >  7:349063/1‑148 (MQ=255)
                                                                                                              cATGGTCAACCTGTGCCGCTGCGGCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAGGGGGAGCAGGCCTGATGAACACGCCAGTCGATACACCCGCCGAACTGCTCAGGCTGCAGCGGGGCGAAACGGTCAACCTGt                                   >  5:87425/1‑148 (MQ=255)
                                                                                                                  gTCAACCTGTGCCGCTGCGGCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAGGGGGAGCAGGCCTGATGAACACGCCAGTCGATACACCCGCCGAACTGCTCAGGCTGCAGCGGGGCGAAACGGTCAACCTGTCACGt                              >  4:146293/1‑149 (MQ=255)
                                                                                                                   tCAACCTGTGCCGCTGCGGCACCTACAACGCCATCCATGCCGCTGTGCATGAc                                                                                                                             <  5:29941/53‑1 (MQ=255)
                                                                                                                   tCAACCTGTGCCGCTGCGGCACCTACAACGCCATCCATGCCGCTGTGCATGAc                                                                                                                             >  6:29941/1‑53 (MQ=255)
                                                                                                                   tCAACCTGTGCCGCTGCGGCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAGGGGGAGCAGGCCTGATGAACACGCCAGTCGATACACCCGCCGAACTGCTCAGGCTGCAGCGGGGCGAAACGGTCAACCTGTCACGTc                             >  7:391490/1‑149 (MQ=255)
                                                                                                                   tCAACCTGTGCCGCTGCGGCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAGGGGGAGCAGGCCTGATGAACACGCCAGTCGATACACCCGCCGAACTGCTCAGGCTGCAGCGGGGCGAAACGGTCAACCTGTCACGTc                             >  5:114619/1‑149 (MQ=255)
                                                                                                                    ttgtaggtgCCGCAGCGGCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAGGGGGAGCAGGCCTGATGAACACGCCAGTCGATACACCCGCCGAACTGCTCAGGCTGCAGCGGGGCGAAACGGTCAACCTGTCACGTCg                            >  1:389242/7‑149 (MQ=255)
                                                                                                                       ccTGTGCCGCTGCGGCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAGGGGGAGCAGGCCTGATGAACACGCCAGTCGATACACCCGCCGAACTGCTCAGGCTGCAGCGGGGCGAAACGGTCaa                                        >  5:36085/1‑134 (MQ=255)
                                                                                                                       ccTGTGCCGCTGCGGCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAGGGGGAGCAGGCCTGATGAACACGCCAGTCGATACACCCGCCGAACTGCTCAGGCTGCAGCGGGGCGAAACGGTCaa                                        <  6:36085/134‑1 (MQ=255)
                                                                                                                       ccTGTGCCGCTGCGGCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAGGGGGAGCAGGCCTGATGAACACGCCAGTCGATACACCCGCCGAACTGCTCAGGCTGCAGCGGGGCGAAACGGTCAACCTGTCACGTCGgcg                         >  7:52419/1‑149 (MQ=255)
                                                                                                                        cTGTGCCGCTGCGGCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAGGGGGAGCAGGCCTGATGAACACGCCAGTc                                                                                       >  2:300478/1‑86 (MQ=255)
                                                                                                                        cTGTGCCGCTGCGGCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAGGGGGAGCAGGCCTGATGAACACGCCAGTc                                                                                       <  1:300478/86‑1 (MQ=255)
                                                                                                                              cGCTGCGGCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAGGGGGAGCAGGCCTGATGAACACGCCAGTCGATACACCCGCCGAACTGCTCAGGCTGCAGCGGGGCGAAACGGTCAACCTGTCACGTCGGCGCTTTCtt                  >  6:193203/1‑149 (MQ=255)
                                                                                                                                cTGCGGCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAGGGGGAGCAGGCCTGATGAACACGCCAGTCGATACACCCGCCGAACTGCTCAGGCTGCAGCGGGGCGAAACGGTCAACCTGTCACGTCGGCGCTTTCTTGc                <  6:329000/149‑1 (MQ=255)
                                                                                                                                     gCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAGGGGGAGCAGGCCTGATGAACACGCCAGTCGATACACCCGCCGAACTGCTCAGGCTGCAGCGGGGCGAAACGGTCAACCTGTCACGTCGGCGCTTTCTTGc                <  3:179442/144‑1 (MQ=255)
                                                                                                                                     gCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAGGGGGAGCAGGCCTGATGAACACGCCAGTCGATACACCCGCCGAACTGCTCAGGCTGCAGCGGGGCGAAACGGTCAACCTGTCACGTCGGCGCTTTCTTGc                >  4:179442/1‑144 (MQ=255)
                                                                                                                                      cACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAGGGGGAGCAGGCCTGATGAACACGCCAGTCGATACACCCGCCGAACTGCTCAGGCTGCAGCGGGGCGAAACGGTCAACCTGTCACGTCGGCGCTTTCTTGCCGGCAc          <  6:322522/149‑1 (MQ=255)
                                                                                                                                          tACAACGCCATCCATGCCGCTGTTCATGACCTGGCCCAGGGGGAGCAGGCCTGATGAACACTCCAGTCGATACACCCGCCGAACTTCACAGGCTGCAGCGGGGCGAAACGGTCAACCTGTCACGTCGGCGCTTTCtt                  <  7:259079/137‑1 (MQ=255)
                                                                                                                                          tACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAGGGGGAGCAGGCCTTATGAACACGCCAGTAGATACACCCGCCGAACTGCACAGTCTGCAGCGGGGCGAAACGGTAAAACTGTAAGGTCGGCGGTTTCtg                  >  8:259079/1‑136 (MQ=255)
                                                                                                                                          tACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAGGGGGAGCAGGCCTGATGAACACGCCAGTCGATACACCCGCCGAACTGCTCAGGCTGCAGCGGGGCGAAACGGTCAACCTGTCACGTCGGCGCTTTCTTGCCGGCACCg        >  7:265409/1‑147 (MQ=255)
                                                                                                                                              aCGCCATCCATGCCGCTGTGCATGACCTGGCCCAGGGGGAGCAGGCCTGATGAACACGCCAGTCGATACACCCGCCGAACTGCTCAGGCTGCAGCGGGGCGAAACGGTCAACCTGTCACGTCGGCGCTTTCTTGCCGGCACCGCAGTCg  >  2:272104/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
GTCGAGCACCAGGTGGCTCAATGCGGCTATTGCCAGTCTGGCCAGGTCATGGCAGCCACGGCGCTGCTCAAGCACACGCCCAAACCCAGCGATGCGCAGATCGATGCAGCCATGGTCAACCTGTGCCGCTGCGGCACCTACAACGCCATCCATGCCGCTGTGCATGACCTGGCCCAGGGGGAGCAGGCCTGATGAACACGCCAGTCGATACACCCGCCGAACTGCTCAGGCTGCAGCGGGGCGAAACGGTCAACCTGTCACGTCGGCGCTTTCTTGCCGGCACCGCAGTCG  >  NC_002947/4117120‑4117410

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