Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,394,256 T→C 5.6% Y264C (TAC→TGC)  nuoL ← NADH:ubiquinone oxidoreductase, membrane subunit L

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,394,2560TC 5.6% 221.7 / 2.8 71Y264C (TAC→TGC) nuoLNADH:ubiquinone oxidoreductase, membrane subunit L
Reads supporting (aligned to +/- strand):  ref base T (41/26);  new base C (2/2);  total (43/28)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

GCGAGAACACGTTTGATGTCGGTCTGTACCAGCGCGGCAAAACCGGCCAGCAGCAGCGTAACCGCCCCGACAATACCCACCAGATGCAGAACTTCCGGCGTCATCAGGAACAGGCCGTGGGTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCGGGCCCGCCATCGCGTCGGCAAGCCATGTCTGCAACGGCAACTGCGCAGATTTACCGACCGCACCGCCCAGCAGCATCAGCGTCGCCCACAT  >  NC_000913/2394121‑2394401
                                                                                                                                       |                                                                                                                                                 
gCGAGAACACGTTTGATGTCGGTCTGTACCAGCGCGGCAAAACCGGCCAGCAGCAGCGTAACCGCCCCGACAATACCCAACTGCAGCAGAACTTCCGGCGTCATCAGGAACAGGCCGTGGGTACTGGCGATCAGGTAGACACCCGCGGt                                                                                                                                      >  1:387034/1‑149 (MQ=255)
          gTTTGATGTCGGTCTGTACCAGCGCGGCAAAACCGGCCAGCAGCAGCGTAACCGCCCCGACAATACCCACCAGATGCAGAACTTCCGGCGTCATCAGGAACAGGCCGTGGGTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGtt                                                                                                                            <  1:271300/149‑1 (MQ=255)
          gTTTGATGTCGGTCTGTACCAGCGCGGCAAAACCGGCCAGCAGCAGCGTAACCGCCCCGACAATACCCACCAGATGCAGAACTTCCGGCGTCATCAGGAACAGGCCGTGGGTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGtt                                                                                                                            <  8:149569/149‑1 (MQ=255)
              gATGTCCGTCTGTACCTGCGCGGCAAAACCGGCCAGCAGCAGCGTAACCGCCCCGACAATACCCACCAGATGCAGAACTTCCGGCGTCATCAGGAACAGGCCGTGGGTACGGGCGATCAGGTAGACTCCCGCGGTTACCATGGTTgcg                                                                                                                         <  4:292145/148‑1 (MQ=255)
                  tCGGTCTGTACCAGCGCGGCAAACCCGGCCAGCAGCATCATAACCGCCCCGAACATCCACAGCATCTGCAGAACTTCCGTCGTCATCAGGAACTGTCAGTGGGTACAGCCGAGCAGGTATACACCCGCGGATCCCATGATTgcgccg                                                                                                                      >  2:91469/1‑147 (MQ=255)
                         gTACCAGCGCGGCAAAACCGGCCAGCAGCAGCGTAACCGCCCCGACAATACCCACCAGATGCAGAACTTCCGGCGTCATCAGGAACAGGCCGTGGGTACGGGCGATCAGGTAGACACCCGCGGt                                                                                                                                      <  7:315217/124‑1 (MQ=255)
                         gTACCAGCGCGGCAAAACCGGCCAGCAGCAGCGTAACCGCCCCGACAATACCCACCAGATGCAGAACTTCCGGCGTCATCAGGAACAGGCCGTGGGTACGGGCGATCAGGTAGACACCCGCGGt                                                                                                                                      >  8:315217/1‑124 (MQ=255)
                            ccAGCGCGGCAAAACCGGCCAGCAGCAGCGTAACCGCCCCGACAATACCCTCCGGAGGCAGAACTGCCGTCGTCAGCAGGACCAGGCGGTCGGTACGGGCGATCAGGTTGACACCCGCGGTTACCATGGTTGCGGCGTGGAGCAgcggg                                                                                                          >  1:70202/1‑147 (MQ=255)
                            ccAGCGCGGCAAAACCGGCCAGCAGCAGCGTAACCGCCCCGACAATACCCACCAGATGCAGAACTTCCGGCGTCATCAGGAACAGGCCGTGGGTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAgcgc                                                                                                           <  6:242751/148‑1 (MQ=255)
                             cAGCGCGGCAAAACCGGCCAGCAGCAGCGTAACCGCCCCGACAATACCCACCAGATGCAGAACTTCCGGCGTCATCAGGAACAGGCCGTGGGTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCgg                                                                                                         <  5:357256/149‑1 (MQ=255)
                                  cgGCAAAACCGGCCAGCAGCAGCGTAACCGCCCCGACAATACCCAGCAGATGCAGAACTTCCGGCGGCATCAGGAACAGGCCGTGGGTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACa                                                                                                    <  1:311118/149‑1 (MQ=255)
                                     cAAAACCGGCCAGCAGCATCGTAAACGCCCACAGAATAACCAGAAGATCGAAAACTAGACGAGTCATCAAGAACAGTGCGTCGGAAAGGGCGAGAAGGTAGACACCCGCGCTTAACAAGAATGAGGGGTGGATCAGCGCGGAGACAGGc                                                                                                 >  6:66550/1‑149 (MQ=255)
                                        aaCCGGCCAGCAGCAGCGTAACCGCCCCGACAATACCCACCAGATGCAGAACTTCCGGCGTCATCAGGAACAGGCCGTGGGTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTc                                                                                              >  3:167092/1‑149 (MQ=255)
                                         aCCGGCCAGCAGCAGCGTAACCGCCCCGACAATACCCACCAGATGCAGAACTTCCGGCGTCATCAGGAACAGGCCGTGGGTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCg                                                                                             >  8:241997/1‑149 (MQ=255)
                                            ggCCAGCAGCAGCGTAACCGCCCCGACAATACCCACCAGATGCAGAACTTCCGGCGTCATCAGGAACAGGCCGTGGGTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCgg                                                                                            >  1:353310/1‑147 (MQ=255)
                                            ggCCAGCAGCAGCGTAACCGCCCCGACAATACCCACCAGATGCAGAACTTCCGGCGTCATCAGGAACAGGCCGTGGGTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCgg                                                                                            <  2:353310/147‑1 (MQ=255)
                                             gCCAGCAGCAGCGTAACCGCCCCGACAATACCCACCAGATCCAGAACTTCCGGCGTCATCAGGAACAGGCCGTGGGTACGGGCCATCAGGTAGACCCCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCGGGcc                                                                                         >  4:317912/1‑149 (MQ=255)
                                                 gcagcagcGTAACCGCCCTGACAATACCCACCAGATGCAGAACTTCCGGCGTCATCAGGAACAGGCCGTGGGTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCGGGCCCGc                                                                                      >  2:371824/1‑148 (MQ=255)
                                                     cagcGTAACCGCCCCGACAATACCCACCAGATGCAGAACTTCCGGCGTCATCAGGAACAGGCCGTGGGTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGcgcg                                                                                                          >  2:227897/1‑124 (MQ=255)
                                                     cagcGTAACCGCCCCGACAATACCCACCAGATGCAGAACTTCCGGCGTCATCAGGAACAGGCCGTGGGTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCGGGCCCGCCATcg                                                                                 >  4:137859/1‑149 (MQ=255)
                                                        cGTAACCGCCCCGACAATACCCACCAGATGCAGAACTTCCGGCGTCATCAGGAACAGGCCGTGGGTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCGGGCCCGCCATCGCGt                                                                              <  5:133546/149‑1 (MQ=255)
                                                             ccGCCCCGACAATACCCACCAGATGCAGAACTTCCGGCGTCATCAGGAACAGGCCGTGGGTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCGGGCCCGCCATCGCGTCGGc                                                                          >  2:273123/1‑148 (MQ=255)
                                                                 cccGACAATACCCACCAGATGCAGAACTTCCGGCGTCATCAGGAACAGGCCGTGGGTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCGGGCCCGCCATCGCGTCGGCAAGcc                                                                     <  5:371271/149‑1 (MQ=255)
                                                                 cccGACAATACCCACCAGATGCAGAACTTCCGGCGTCATCAGGAACAGGCCGTGGGTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCGGGCCCGCCATCGCGTCGGCAAGc                                                                      >  6:371271/1‑148 (MQ=255)
                                                                            ccaccaGATGCAGAACTTCCGGCGTCATCAGGAACAGGCCGTGGTTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTGGCGGCTTGGATCCGCGCGGAGACCGGCGTCGGGCCCGCCATGCCTTCGGCACGCCATGTCTCCAAc                                                          >  6:5995/1‑149 (MQ=255)
                                                                              accaGATGCAGAACTTCCGGCGTCATCAGGAACAGGCCGTGGGTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGc                                                                                                 >  3:188059/1‑108 (MQ=255)
                                                                              accaGATGCAGAACTTCCGGCGTCATCAGGAACAGGCCGTGGGTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGc                                                                                                 <  4:188059/108‑1 (MQ=255)
                                                                              accaGATGCAGAACTTCCGGCGTCATCAGGAACAGGCCGTGGGTACGGGCGATCAGGCAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCggg                                                                                           <  6:217748/114‑1 (MQ=255)
                                                                              accaGATGCAGAACTTCCGGCGTCATCAGGAACAGGCCGTGGGTACGGGCGATCAGGCAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCggg                                                                                           <  6:212261/114‑1 (MQ=255)
                                                                              accaGATGCAGAACTTCCGGCGTCATCAGGAACAGGCCGTGGGTACGGGCGATCAGGCAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCggg                                                                                           >  5:217748/1‑114 (MQ=255)
                                                                              accaGATGCAGAACTTCCGGCGTCATCAGGAACAGGCCGTGGGTACGGGCGATCAGGCAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCggg                                                                                           >  5:212261/1‑114 (MQ=255)
                                                                                  gATGCAGAACTTCCGGCGTCATCAGGAACAGGCCGTGGGTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCGGGCCCGCCATCGCGTCGGCAAGCCAt                                                                   >  3:287790/1‑134 (MQ=255)
                                                                                  gATGCAGAACTTCCGGCGTCATCAGGAACAGGCCGTGGGTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCGGGCCCGCCATCGCGTCGGCAAGCCAt                                                                   <  4:287790/134‑1 (MQ=255)
                                                                                   aTGCAGAACTTCCGGCGTCATCAGGAACAGGCCGTGGGTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCGGGCCCGCCATCGCGTCGGCAAGCCATGTCTGCAACGGCAACt                                                   <  7:241997/149‑1 (MQ=255)
                                                                                       aGAACTTCCGGCGTCATCAGGAACAGGCCGTGGGTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCGGGCCCGCCAt                                                                                   >  2:377132/1‑113 (MQ=255)
                                                                                       aGAACTTCCGGCGTCATCAGGAACAGGCCGTGGGTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCGGGCCCGCCAt                                                                                   <  1:377132/113‑1 (MQ=255)
                                                                                       aGAACTTCCGGCGTCATCAGGAACAGGCCGTGGGTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCGGGCCCGCCATc                                                                                  >  7:104602/1‑114 (MQ=255)
                                                                                       aGAACTTCCGGCGTCATCAGGAACAGGCCGTGGGTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCGGGCCCGCCATc                                                                                  <  8:104602/114‑1 (MQ=255)
                                                                                       aGAACTTCCGGCGTCATCAGGAACAGGCCGTGGGTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCGGGCCCGCCATCGCGTCGGCAAGCCATGTCTGCAACGGCAACTgcg                                                >  4:347387/1‑148 (MQ=255)
                                                                                       aGAACTTCCGGCGTCATCAGGAACAGGCCGTGGGTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCGGGCCCGCCATCGCGTCGGCAACCCATGTCTGCAACGGCAACTgcg                                                >  8:145981/1‑148 (MQ=255)
                                                                                        gAACTTCCGGCGTCATCAGGAACAGGCCGTGGGTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCGGGCCCGCCATCGCGTCGGCAAGCCATGTCTGCAACGGCAACTgcg                                                >  6:308870/1‑147 (MQ=255)
                                                                                            ttCCGGCGTCATCAGGAACAGGCCGTGGGTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCGGGCCCGCCATCGCGTCGGCAAGCCATGTCTGCAACGGCAACTGCGCAGa                                            >  6:87975/1‑147 (MQ=255)
                                                                                            ttCCGGCGTCATCAGGAACAGGCCGTGGGTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCGGGCCCGCCATCGCGTCGGCAAGCCATGTCTGCAACGGCAACTGCGCAGa                                            <  5:87975/147‑1 (MQ=255)
                                                                                                ggCGTCATCAGGAACAGGCCGTGGGTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCGGGCCCGCCATCGCGTCGGCAAGCCATGTCTGCAACGGCAACTGCGCAGATTTacc                                      >  2:212289/1‑149 (MQ=255)
                                                                                                ggCGTCATCAGGAACAGGCCGTGGGTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCGGGCCCGCCATCGCGTCGGCAAGCCATGTCTGCAACGGCAACTGCGCAGATTTacc                                      >  2:118054/1‑149 (MQ=255)
                                                                                                    tcatcaGGAACAGGCCGTGGGTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCGGGCCCGCCATc                                                                                  >  4:91267/1‑101 (MQ=255)
                                                                                                    tcatcaGGAACAGGCCGTGGGTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCGGGCCCGCCATc                                                                                  <  3:91267/101‑1 (MQ=255)
                                                                                                       tcaGGAACAGGCCGTGGGTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCGGGCCCGCCATCGCGTCGGCAAGCCATGTCTGCAACGGCAACTGCGCAGATTTACCGaccgca                               >  3:92777/1‑149 (MQ=255)
                                                                                                       tcaGGAACAGGCCGTGGGTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCGGGCCCGCCATCGCGTCGGCAAGCCATGTCTGCAACGGCAACTGCGCAGATTTACCGaccgca                               >  3:155037/1‑149 (MQ=255)
                                                                                                       tcaGGAACAGGCCGTGGGTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCGGGCCCGCCATCGCGTCGGCAAGCCATGTCTGCAACGGCAACTGCGCAGATTTACCGaccgca                               >  3:74875/1‑149 (MQ=255)
                                                                                                         aGGAACAGGCCGTGGGTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCGGGCCCGCCATCGCGTCGGCAAGCCATGTCTGCAACGGCAACTGCGCAGATTTACCGaccgcacc                             >  3:96842/1‑149 (MQ=255)
                                                                                                             aCAGGGCGTGGGTACGGGCGATCAGGTAGACACCCGCGCTTACAATGCTTGCGGCGTGGATCAGAGAGGAGACAGGCGTCGGGCCCGCCATCGCGTCGGCAAGCAATGAATGAAACGCCAACTGCGCAGAATTACCGACAGCACCGccc                         >  4:46835/1‑149 (MQ=255)
                                                                                                             aCAGGCCGTGGGTACGGGCGATCAGGTAGACAACCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCGGGCCCGCCATCGCGTCGGCAAGCCATGTCTGCACCGGCAAATGCGCAGATTTACCGACCGCACCGcc                          >  2:313317/1‑148 (MQ=255)
                                                                                                                 gCCGTGGGTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCTCGGAGACAGGCGTCGGGCCCGCCATCGCGTCGGCAAGCCATGTCTGCAACg                                                         >  6:70600/1‑113 (MQ=255)
                                                                                                                 gCCGTGGGTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCTCGGAGACAGGCGTCGGGCCCGCCATCGCGTCGGCAAGCCATGTCTGCAACg                                                         <  5:70600/113‑1 (MQ=255)
                                                                                                                 gCCGTGGGTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCGGGCCCGCCATCGCGTCGGCAAGCCATGTCTGCAACGGCAACTGCGCAGATTTACCGACCGCACCGCCcagca                     <  3:317912/149‑1 (MQ=255)
                                                                                                                 gCCGTGGGTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCGGGCCCGCCATCGCGTCGGCAAGCCATGTCTGCAACGGCAACTGCGCAGATTTACCGACCGCACCGCCcagca                     >  2:359192/1‑149 (MQ=255)
                                                                                                                 gCCGTGGGTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCGGGCCCGCCATCGCGTCGGCAAGCCATGTCTGCAACGGCAACTGCGCAGATTTACCGACCGCACCGCCcagca                     <  5:5995/149‑1 (MQ=255)
                                                                                                                 gCCGTGGCTACGGGCGATCACGTAGACACCCGCGGTTACCATGGTTGCGTCGTTGATCAGCGCGGAGAAAGGCGTCGCGCCCGCCAAAGACTCGGCAAGCCATGTCTGCAACGGCAACTCCGCAGATTTACCTACCGCACAGCCcagca                     >  8:356503/1‑149 (MQ=255)
                                                                                                                  ccGTGGGTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCGGGCCCGCCATCGCGTCGGCAAGCCATGTCTGCAACGGCAACTGCGCAGATTTACCGACCGCACCGCCcagcag                    <  3:137859/149‑1 (MQ=255)
                                                                                                                       ggTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCGGGCCCGCCATCGCGTCGGCAAGCCATGTCTGCAACGGCAACTGCGCAGATTTACCGACCGCACCGCCCAGCAGCATCa               <  5:308870/149‑1 (MQ=255)
                                                                                                                       ggTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCGGGCCCGCCATCGCGTCGGCAAGCCATGTCTGCAACGGCAACTGCGCAGATTTACCGACCGCACCGCCCAGCAGCATCa               <  1:212289/149‑1 (MQ=255)
                                                                                                                         tACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCGGGCCCGCCATCGCGTCGGCAAGCCATGTCTGCAACGGCAACTGCGCa                                              >  4:59682/1‑116 (MQ=255)
                                                                                                                         tACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCGGGCCCGCCATCGCGTCGGCAAGCCATGTCTGCAACGGCAACTGCGCAGATTTACCGACCGCACCGCCCAGCAGCATCAGc             <  1:273123/149‑1 (MQ=255)
                                                                                                                            gggCGATCAGGTATACACACGAGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCGGGCACGACATCGAGTAGGCAAGCCATGTCTGCAACGGCAACTGCGCAGATTTACCGACCGCACCGCCCAGCAGCATCATCGTc          >  6:340421/1‑149 (MQ=255)
                                                                                                                            gggCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGGGTGGATCAGCGCGGAGACAGGCGTCGGGCCCGCCATCGCGTCGGCAAGCCATGTCTGCAACGGCAACTGCGCAGATTTACCGACCGCACCGCCCAGCAGCATCAGCGTc          >  3:167601/1‑149 (MQ=255)
                                                                                                                             ggCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCGGGCCCGCCATCGCGTCGGCAAGCCATGTCTGCAACGGCAACTGCGCAGATTTACCGACCGCACCGCCCAGCAGCATCAGCGTCg         <  1:359192/149‑1 (MQ=255)
                                                                                                                               cGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCGGGCCCGCCATCGCTTCGGCAAGCCATGTCTGCAACGGCAACTGCGCAGATTTACCGACATCACCGCCCAGCAGCATCAGCGTCGcc       <  4:167092/149‑1 (MQ=255)
                                                                                                                               cGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCGGGCCCGCCATCGCGTCGGCAAGCCATGTCTGCAACGGCAACTGCGCAGATTTACCGACCGCACCGCCCAGCAGCATCAGCGTCGcc       >  4:12413/1‑149 (MQ=255)
                                                                                                                                      gTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCGGGCCCGCCATCGCGTCGGCAAGCCATGTCTGCAACGGCAACTGCGCAGATTTACCGACCGCACCGCCCAGCAGCATCAGCGTCGCCCACAt  >  7:57949/1‑147 (MQ=255)
                                                                                                                                      gTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCGGGCCCGCCATCGCGTCGGCAAGCCATGTCTGCAACGGCAACTGCGCAGATTTACCGACCGCACCGCCCAGCAGCATCAGCGTCGCCCACAt  >  8:310372/1‑147 (MQ=255)
                                                                                                                                       |                                                                                                                                                 
GCGAGAACACGTTTGATGTCGGTCTGTACCAGCGCGGCAAAACCGGCCAGCAGCAGCGTAACCGCCCCGACAATACCCACCAGATGCAGAACTTCCGGCGTCATCAGGAACAGGCCGTGGGTACGGGCGATCAGGTAGACACCCGCGGTTACCATGGTTGCGGCGTGGATCAGCGCGGAGACAGGCGTCGGGCCCGCCATCGCGTCGGCAAGCCATGTCTGCAACGGCAACTGCGCAGATTTACCGACCGCACCGCCCAGCAGCATCAGCGTCGCCCACAT  >  NC_000913/2394121‑2394401

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