Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,440,441 T→C 5.0% H115R (CAT→CGT)  pfo ← pyruvate‑flavodoxin oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,440,4410TC 5.0% 247.1 / 2.1 79H115R (CAT→CGT) pfopyruvate‑flavodoxin oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (36/39);  new base C (2/2);  total (38/41)
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.74e-01

AATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGATCATCAGCAGCAAACCCTGCGATGACGTAAACGATGTTGAAAGGGCACCCGTCTGCAAAGCGCCA  >  NC_000913/1440293‑1440581
                                                                                                                                                    |                                                                                                                                            
aaTAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCAt                                                                                                                                              <  2:430119/149‑1 (MQ=255)
    aaTGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCAtgtgt                                                                                                                                          <  1:218367/149‑1 (MQ=255)
     aTGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCAtgtgtg                                                                                                                                         >  2:360112/1‑149 (MQ=255)
     aTGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCAtgtgtg                                                                                                                                         >  7:427222/1‑149 (MQ=255)
        gCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGc                                                                                                                                       >  5:295134/1‑148 (MQ=255)
                 ttttCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCAtgtgt                                                                                                                                          >  2:302669/1‑136 (MQ=255)
                 ttttCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCAtgtgt                                                                                                                                          <  1:302669/136‑1 (MQ=255)
                   ttCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGt                                                                                                                             <  2:323816/147‑1 (MQ=255)
                                     gAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACaaa                                                                                                         >  1:398747/1‑149 (MQ=255)
                                     gAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACaaa                                                                                                         >  2:84979/1‑149 (MQ=255)
                                        aTGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACgg                                                                                                      <  6:364985/149‑1 (MQ=255)
                                         tGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGgt                                                                                                     <  2:398747/149‑1 (MQ=255)
                                             agCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTa                                                                                                 <  7:5633/149‑1 (MQ=255)
                                               cAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAAc                                                                                                        <  1:177369/140‑1 (MQ=255)
                                               cAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAAc                                                                                                        >  2:177369/1‑140 (MQ=255)
                                               cAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGt                                                                                               >  1:90043/1‑149 (MQ=255)
                                                 aaGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTc                                                                                             <  6:295134/149‑1 (MQ=255)
                                                   gTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGc                                                                                           <  2:19024/149‑1 (MQ=255)
                                                      ttGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTg                                                                                        >  7:46360/1‑149 (MQ=255)
                                                        gcgcTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCt                                                                                         <  8:427222/146‑1 (MQ=255)
                                                              ccTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTg                                                                                >  8:22516/1‑149 (MQ=255)
                                                               cTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTg                                                                                <  1:84979/148‑1 (MQ=255)
                                                               cTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTg                                                                                >  2:450146/1‑148 (MQ=255)
                                                                tGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGtt                                                                                                  <  5:384712/129‑1 (MQ=255)
                                                                tGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGtt                                                                                                  >  6:384712/1‑129 (MQ=255)
                                                                tGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTa                                                                              >  8:243443/1‑149 (MQ=255)
                                                                 ggACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTa                                                                              >  2:18181/1‑148 (MQ=255)
                                                                 ggACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTACATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGgt                                                                                                     >  2:162511/1‑125 (MQ=255)
                                                                 ggACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTACATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGgt                                                                                                     <  1:162511/125‑1 (MQ=255)
                                                                  gACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTAc                                                                             <  7:22516/148‑1 (MQ=255)
                                                                  gACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACa                                                                            <  1:18181/149‑1 (MQ=255)
                                                                             cacacaACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGat                                                                  >  2:94932/1‑148 (MQ=255)
                                                                             cacacaACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGat                                                                  <  1:360112/148‑1 (MQ=255)
                                                                              acacaACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCgg                                                                     >  2:238426/1‑144 (MQ=255)
                                                                              acacaACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCgg                                                                     <  1:238426/144‑1 (MQ=255)
                                                                               cacaACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGAGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGatcat                                                               <  2:94701/149‑1 (MQ=255)
                                                                               cacaACATCGCGCAACCCGTCTGGCGCACCGCCATAACCTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGGGGTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGatcat                                                               <  5:382880/149‑1 (MQ=255)
                                                                                acaacaTCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGatcatc                                                              <  6:398152/149‑1 (MQ=255)
                                                                                acaacaTCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGatcat                                                               >  1:227556/1‑148 (MQ=255)
                                                                                acaacaTCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGatcat                                                               >  3:320618/1‑148 (MQ=255)
                                                                                  aacaTCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGATCATcag                                                            >  2:3154/1‑149 (MQ=255)
                                                                                  aacaTCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGATCATcag                                                            >  7:439365/1‑149 (MQ=255)
                                                                                  aacaTCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGATCATcag                                                            >  4:387010/1‑149 (MQ=255)
                                                                                    caTCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGATCATcagca                                                          <  1:450146/149‑1 (MQ=255)
                                                                                    caTCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGATCATcag                                                            >  5:31226/1‑147 (MQ=255)
                                                                                     aTCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGATCATcagcag                                                         <  1:3154/149‑1 (MQ=255)
                                                                                     aTCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGGCGGGATCATcagcag                                                         >  3:96749/1‑149 (MQ=255)
                                                                                        gcgcAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGATCATCAGCAGCaa                                                      <  2:469233/149‑1 (MQ=255)
                                                                                         cgcAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGATCATCAGCAGCaaa                                                     <  6:31226/149‑1 (MQ=255)
                                                                                            aaCCCGTCTGGCGCACCGCCATAACGTCGGAATGCTCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTACGCCTGCCAGTTTGTACAGCGTCGGGATCATCAGCAGCAAAccc                                                  >  2:82352/1‑149 (MQ=255)
                                                                                            aaCCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACAGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGATCATCAGCAGCAAAccc                                                  >  2:346423/1‑149 (MQ=255)
                                                                                             aCCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGatcat                                                               >  4:178633/1‑135 (MQ=255)
                                                                                             aCCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGatcat                                                               <  3:178633/135‑1 (MQ=255)
                                                                                              cccGTCTGTCGCACCGCCATAACGTCGTAATGATCGCCAAAAATAGAGAGGTCATGTGTGGCAACGGTACGCGCCGCTACATGCAGGACACACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGATCATCAGCAGCAAACCCTg                                                <  4:48022/149‑1 (MQ=255)
                                                                                               ccGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGATCATCAGCAGCAAACCCTGc                                               >  5:32096/1‑149 (MQ=255)
                                                                                               ccGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGATCATCAGCAGCAAACCCTGc                                               <  8:439365/149‑1 (MQ=255)
                                                                                                  tCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGATCATCAGCAGCAAACCCTGCGAt                                            <  4:320618/149‑1 (MQ=255)
                                                                                                     ggCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGATCATCAGCAGCAAACCCTGCGATGa                                          >  1:83035/1‑148 (MQ=255)
                                                                                                      gcgcACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGATCATCAGCAGCAAACCCTGCGATGa                                          <  2:90043/147‑1 (MQ=255)
                                                                                                        gcACGGCCATAACGTCGGAATGATCGCCTAAAATAGAGAGTGCTTGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGATCATCAGCAGCAAACCCTGCGATg                                           <  2:149991/144‑1 (MQ=255)
                                                                                                        gcACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGATCATCAGCAGCAAACCCTGCGATg                                           >  1:149991/1‑144 (MQ=255)
                                                                                                        gcACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGATCATCAGCAGCAAACCCTGCGATGACGt                                       >  5:299935/1‑148 (MQ=255)
                                                                                                        gcACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTCGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGATCATCAGCAGCAAACCCTGCGATGACGt                                       >  5:384689/1‑148 (MQ=255)
                                                                                                            cGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGATCATCAGCAGCAAACCCTGCGATGACGTAAACg                                  <  6:32096/149‑1 (MQ=255)
                                                                                                               cATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGATCATCAGCAGCAAACCCTGCGATGACGTAAACGATg                               >  7:280901/1‑149 (MQ=255)
                                                                                                                aTAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGATCATCAGCAGCAAACCCTGCGATGACGTAAACGATg                               <  7:97036/148‑1 (MQ=255)
                                                                                                                aTAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGATCATCAGCAGCAAACCCTGCGATGACGTAAACGATg                               >  8:97036/1‑148 (MQ=255)
                                                                                                                       cGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGatcat                                                               >  1:149334/1‑109 (MQ=255)
                                                                                                                       cGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGatcat                                                               <  2:149334/109‑1 (MQ=255)
                                                                                                                       cGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGATCATCAGCAGCAAACCCTGCGATGACGTAAACGATGTTGAAAgg                       <  1:94932/149‑1 (MQ=255)
                                                                                                                         gAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGATCATCAGCAGCAAACCCTGCGATGACGTAAACGATGTTGAAAGGGc                     <  7:243443/149‑1 (MQ=255)
                                                                                                                                cGCCAAAAATAGAGAGTGCACGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGt                                                                        >  8:239657/1‑91 (MQ=255)
                                                                                                                                cGCCAAAAATAGAGAGTGCACGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGt                                                                        >  8:246251/1‑91 (MQ=255)
                                                                                                                                cGCCAAAAATAGAGAGTGCACGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGt                                                                        <  7:246251/91‑1 (MQ=255)
                                                                                                                                cGCCAAAAATAGAGAGTGCACGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGt                                                                        <  7:239657/91‑1 (MQ=255)
                                                                                                                                 gCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGATCATCAGCAGCAAACCCTGCGATGACGTAAACGATGTTGAAAGGGCACCCGTCt             >  6:124843/1‑149 (MQ=255)
                                                                                                                                  ccAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGATCATCAGCAGCAAACCCTGCGATGACGTAAACGATGTTGAAAGGGCACCCGTCTg            <  8:231936/149‑1 (MQ=255)
                                                                                                                                    aaaaaTAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGATCATCAGCAGCAAACCCTGCGATGACGTAAACGATGTTGAAAGGGCACCCGTCTGc           >  8:336056/1‑148 (MQ=255)
                                                                                                                                      aaaTAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTATTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGATCATCAGCAGCAAACCCTGCGATGACGTAAACGATGTTGAAAGGGCACCCGTCTGCaaa        <  4:93633/149‑1 (MQ=255)
                                                                                                                                         tAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGa                                                                   >  6:245974/1‑87 (MQ=255)
                                                                                                                                         tAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGa                                                                   <  5:245974/87‑1 (MQ=255)
                                                                                                                                           gagagTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGATCATCAGCAGCAAACCCTGCGATGACGTAAACGATGTTGAAAGGGCACCCGTCTGCAAAGCGcc   <  3:387010/149‑1 (MQ=255)
                                                                                                                                            agagTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGATCATCAGCAGCAAACCCTGCGATGACGTAAACGATGTTGAAAGGGCACCCGTCTGCAAAGCGCCa  <  6:299935/149‑1 (MQ=255)
                                                                                                                                                tGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGATCATCAGCAGCAAACCCTGCGATGACGTAAACGATGTTGaaa                         >  2:358279/1‑122 (MQ=255)
                                                                                                                                                tGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGATCATCAGCAGCAAACCCTGCGATGACGTAAACGATGTTGaaa                         <  1:358279/122‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
AATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAACGGTGTTAGTTCGCCTGCCAGTTTGTACAGCGTCGGGATCATCAGCAGCAAACCCTGCGATGACGTAAACGATGTTGAAAGGGCACCCGTCTGCAAAGCGCCA  >  NC_000913/1440293‑1440581

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