Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,391,414 C→G 12.5% R31P (CGA→CCA)  nuoN ← NADH:ubiquinone oxidoreductase, membrane subunit N

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,391,4140CG12.5% 297.9 / 32.5 97R31P (CGA→CCA) nuoNNADH:ubiquinone oxidoreductase, membrane subunit N
Reads supporting (aligned to +/- strand):  ref base C (34/50);  new base G (7/5);  total (42/55)
Fisher's exact test for biased strand distribution p-value = 3.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.42e-01

AGCCCGGTGTAAAGCATGGCGAAACCATCAACGCGCATCAGCGGCGTAACGTCCATAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGCCTTGTAGTAGTAACGGAATTAACAAACCACTGCTGGATATTGC  >  NC_000913/2391270‑2391561
                                                                                                                                                |                                                                                                                                                   
aGCCCGGTGTAAAGCATGGCGAAACCATCAACGCGCATCAGCGGCGTAACGTCCATAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGcc                                                                                                                                                  >  3:133464/1‑148 (MQ=255)
   ccGGTGTAAAGCATGGCGAAACCATCAACGCGCATCAGCGGCGTAACGTCCATAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGc                                                                                                                                              >  2:448363/1‑149 (MQ=255)
     ggTGTGAAGCATGGCGAAACCATCAACGCGCATCAGCGGCGTAACGTCCATAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTGGCGACGCaa                                                                                                                                            <  2:701516/149‑1 (MQ=255)
      gtgtAAAGCATGGCGAAACCATCAACGCGCATCAGCGGCGTAACGTCCATAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTGGCGACGCAAt                                                                                                                                           <  2:820315/149‑1 (MQ=255)
       tgtAAAGCATGGCGAAACCATCAACGCGCATCAGCGGCGTAACGTCCATAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATg                                                                                                                                          <  1:764489/149‑1 (MQ=255)
         tAAAGCATGGCGAAACCATCAACACGCATCAGCGGCGTAACGTCCATAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGga                                                                                                                                        <  2:179426/149‑1 (MQ=255)
           aaGCATGGCGAAACCATCAACGCGCATCAGCGGCGTAACGTCCATAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTATCGTTGAGGAAATGATTGCGTCGCCACGCAATGgag                                                                                                                                       <  4:1002142/148‑1 (MQ=255)
           aaGCATGGCGAAACCATCAACGCGCATCAGCGGCGTAACGTCCATAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGTGTCGCCACGCAATGgaga                                                                                                                                      <  2:281138/149‑1 (MQ=255)
           aaGCATGGCGAAACCATCAACGCGCATCAGCGGCGTAACGTCCATAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTGGCGACGCAAtcattt                                                                                                                                      <  1:757320/149‑6 (MQ=255)
           aaGCATGGCGAAACCATCAACGCGCATCAGCGGCGTAACGTCCATAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGgag                                                                                                                                       <  4:895/148‑1 (MQ=255)
            aGCATGGCGAAACCATCAACGCGCATCAGCGGCGTAACGTCCATAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGagag                                                                                                                                     <  2:292798/149‑1 (MQ=255)
              cATGGCGAAACCATCAACGCGCATCAGCGGCGTAACGTCCATAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTATCCTTGAAGAAATGATCCCGTCGGTCCGCAGTGgagatg                                                                                                                                    <  4:935956/148‑3 (MQ=255)
              cATGGCGAAACCATCAACGCGCATCAGCGGCGTAACGTCCATAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCa                                                                                                                                   <  2:581107/149‑1 (MQ=255)
                  gCGAAACCATCAACGCGCATCAGCGGCGTAACGTCCATAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATcac                                                                                                                               <  2:356549/149‑1 (MQ=255)
                     aaaCCATCAACGCGCATCAGCGGCGTAACGTCCATAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAac                                                                                                                            <  1:1038623/149‑1 (MQ=255)
                      aaCCATCAACGCGCATCAGCGGCGTAACGTCCATAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAac                                                                                                                            <  3:71622/148‑1 (MQ=255)
                            cAACGCGCATCAGCGGCGTAACGTCCATAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAaccaccac                                                                                                                      >  4:993840/1‑148 (MQ=255)
                              aCGCGCATCAGCGGCGTAACGTCCATAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCg                                                                                                                    >  3:398976/1‑148 (MQ=255)
                                      cAGCGGCGTAACGTCCATAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCcg                                                                                                            <  3:457566/148‑1 (MQ=255)
                                      cAGCGGCGTAACGTCCATAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCcg                                                                                                            >  3:291623/1‑148 (MQ=255)
                                         cggcgTAACGTCCATAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGCTAAGCCCAATAACCGAAAGCGTAGCGTTGAGGAGATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCcgacg                                                                                                         >  4:401985/1‑148 (MQ=255)
                                          ggcgTAACGTCCATAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCcgacga                                                                                                        <  4:655853/148‑1 (MQ=255)
                                              tAACGTCCATAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAg                                                                                                    >  3:385720/1‑148 (MQ=255)
                                               aaCGTCCATAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAAt                                                                                                                                           <  1:759313/108‑1 (MQ=255)
                                               aaCGTCCATAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAAt                                                                                                                                           >  2:759313/1‑108 (MQ=255)
                                                 cGTCCATAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCaa                                                                                                 <  4:503614/148‑1 (MQ=255)
                                                  gTCCATAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACg                                                                                               >  2:905006/1‑149 (MQ=255)
                                                   tCCATAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACg                                                                                               >  4:751900/1‑148 (MQ=255)
                                                       tAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAg                                                                                           >  3:279562/1‑148 (MQ=255)
                                                        aGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGc                                                                                          <  4:346239/148‑1 (MQ=255)
                                                        aGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGc                                                                                          <  4:398976/148‑1 (MQ=255)
                                                        aGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCa                                                                                         >  1:649449/1‑149 (MQ=255)
                                                            cccgccTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGc                                                                                     <  1:335757/149‑1 (MQ=255)
                                                            cccgccTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGc                                                                                     >  1:39083/1‑149 (MQ=255)
                                                             ccgccTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGgag                                                                                                                                       <  2:353524/98‑1 (MQ=255)
                                                             ccgccTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGgag                                                                                                                                       >  1:353524/1‑98 (MQ=255)
                                                                 cTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCa                                                                                <  1:524510/149‑1 (MQ=255)
                                                                  tGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAg                                                                               <  1:153894/149‑1 (MQ=255)
                                                                      caacaaACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGtt                                                                            >  3:293130/1‑148 (MQ=255)
                                                                       aacaaACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTGGCGACGCAAtc                                                                                                                                          >  2:524990/1‑84 (MQ=255)
                                                                       aacaaACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTGGCGACGCAAtc                                                                                                                                          <  1:524990/85‑2 (MQ=255)
                                                                       aacaaACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGtttt                                                                          >  1:33425/1‑149 (MQ=255)
                                                                       aacaaACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGtttt                                                                          >  1:192450/1‑149 (MQ=255)
                                                                       aacaaACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGttt                                                                           >  3:248265/1‑148 (MQ=255)
                                                                         caaACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTtgt                                                                        <  2:130347/149‑1 (MQ=255)
                                                                         caaACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTtgt                                                                        <  2:1045472/149‑1 (MQ=255)
                                                                         caaACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTtg                                                                         >  3:870845/1‑148 (MQ=255)
                                                                         caaACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTtg                                                                         <  3:177226/148‑1 (MQ=255)
                                                                         caaACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTtg                                                                         <  4:870351/148‑1 (MQ=255)
                                                                          aaaCCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGTATCACAACCACCTCCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTtgt                                                                        <  4:1007189/148‑1 (MQ=255)
                                                                          aaaCCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTtgtg                                                                       <  1:71671/149‑1 (MQ=255)
                                                                            aCCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTgg                                                                      <  4:459540/148‑1 (MQ=255)
                                                                            aCCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGa                                                                     <  1:480434/149‑1 (MQ=255)
                                                                             ccAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAg                                                                    <  1:448363/149‑1 (MQ=255)
                                                                                gagCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAg                                                                               >  2:948117/1‑135 (MQ=255)
                                                                                gagCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAg                                                                               <  1:948117/135‑1 (MQ=255)
                                                                                  gCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGATATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTAt                                                                >  3:92118/1‑148 (MQ=255)
                                                                                  gCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTAtt                                                               <  2:33425/149‑1 (MQ=255)
                                                                                   cGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTAtt                                                               >  4:570570/1‑148 (MQ=255)
                                                                                    gagatcAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTATCCACGCAATGGAGAGCATGACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTg                                                              >  4:902037/6‑148 (MQ=255)
                                                                                      aaCCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTGGCGACGCAAt                                                                                                                                           <  4:348817/69‑1 (MQ=255)
                                                                                      aaCCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTGGCGACGCAAt                                                                                                                                           >  3:348817/1‑69 (MQ=255)
                                                                                      aaCCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGtttt                                                                          <  2:978530/134‑1 (MQ=255)
                                                                                      aaCCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGtttt                                                                          >  1:978530/1‑134 (MQ=255)
                                                                                       aCCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCAt                                                          >  2:694149/1‑149 (MQ=255)
                                                                                        ccAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATg                                                         >  1:183552/1‑149 (MQ=255)
                                                                                         cAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATg                                                         <  4:114855/148‑1 (MQ=255)
                                                                                          aGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATgg                                                        <  3:579972/148‑1 (MQ=255)
                                                                                            cgccgcGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGa                                                     <  1:119148/149‑1 (MQ=255)
                                                                                                 cgTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGAGATTTAc                                                <  2:573936/149‑1 (MQ=255)
                                                                                                       gCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGCCt                                           <  3:570570/148‑1 (MQ=255)
                                                                                                         ccAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGCCTTgt                                        <  1:694149/149‑1 (MQ=255)
                                                                                                          cAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGCCTTgta                                       <  1:54788/149‑1 (MQ=255)
                                                                                                            aTAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGGCTTgtag                                      <  4:92118/148‑1 (MQ=255)
                                                                                                                ccGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGCCTTgtagtagta                                 >  2:76301/1‑149 (MQ=255)
                                                                                                                  gagaGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGCCTTGTAGTAGTaa                                >  4:625076/1‑148 (MQ=255)
                                                                                                                  gagaGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGCCTTGTAGTAGTaa                                <  4:59159/148‑1 (MQ=255)
                                                                                                                  gagaGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGCCTTGTAGTAGTAAc                               <  1:905006/149‑1 (MQ=255)
                                                                                                                  gagaGCGTAGCGTTGAGGAAATGAATGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGAACAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGCCTTGTAGTAGTAAc                               >  2:689924/1‑149 (MQ=255)
                                                                                                                         tagcgtTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTAGTGGAGTTATTGTCATGGCGATTTACGGCCTTGTAGGAGTAACGGAAtt                         <  4:385720/148‑1 (MQ=255)
                                                                                                                          agcgtTGAGGAAATGATTGCGTCGCAACGCAATGGAGAGCATCACAACCACCACCATCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGCCTTGTAGTAGTAGCGGAATTa                        >  3:656201/1‑148 (MQ=255)
                                                                                                                              ttGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGCCTTGTAGTAGTAACGGAATTAACaa                    <  3:15636/148‑1 (MQ=255)
                                                                                                                                  ggAAATGATTGCGTCGCCACGCAAGGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGGAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGCCTTGTAGTAGTAACGGAATTAACAAACCAc               >  1:544005/1‑149 (MQ=255)
                                                                                                                                   gAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGCCTTGTAGTAGTAACGGAATTAACAAACCAct              >  1:1030785/1‑149 (MQ=255)
                                                                                                                                   gAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGCCTTGTAGTAGTAACGGAATTAACAAACCAct              >  2:943898/1‑149 (MQ=255)
                                                                                                                                    tctccaTTGCGTGGCGACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGCCTTGTAGTAGTAACGGAATTAACAAACCActg             >  2:45516/6‑149 (MQ=255)
                                                                                                                                    aaaTGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGCCTTGTAGTAGTAACGGAATTAACAAACCAct              <  4:248265/148‑1 (MQ=255)
                                                                                                                                     aaTGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGCCTTGTAGTAGTAACGGAATTAACAAACCActg             >  4:597268/1‑148 (MQ=255)
                                                                                                                                      tccaTTGCGTGGCGACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGCCTTGTAGTAGTAACGGAATTAACAAACCActgc            >  3:1047062/4‑148 (MQ=255)
                                                                                                                                        caTTGCGTGGCGACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGCCTTGTAGTAGTAACGGAATTAACAAACCActgctg          >  3:720450/2‑148 (MQ=255)
                                                                                                                                        caTTGCGTGGCGACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGCCTTGTAGTAGTAACGGAATTAACAAACCACTGCTgg         >  1:854361/2‑149 (MQ=255)
                                                                                                                                         aTTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGCCTTGTAGTAGTAACGGAATTAACAAACCACTGCTgg         <  4:293130/148‑1 (MQ=255)
                                                                                                                                         aTTGCGTCACCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGCCTTGTAGTAGTAACGGAATTAACAAACCACTGCTGGa        <  2:881440/149‑1 (MQ=255)
                                                                                                                                          ttGCGTGGCGACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGCCTTGTAGTAGTAACGGAATTAACAAACCACTGCTGGa        >  3:321457/1‑148 (MQ=255)
                                                                                                                                               tCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGCCTTGTAGTAGTAACGGAATTAACAAACCACTGCTGGATATTGc  <  2:39083/149‑1 (MQ=255)
                                                                                                                                               tCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGCCTTGTAGTAGTAACGGAATTAACAAACCACTGCTGGATATTGc  >  2:786161/1‑149 (MQ=255)
                                                                                                                                               tCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGCCTTGTAGTAGTAACGGAATTAACAAACCACTGCTGGATATTGc  <  2:152420/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
AGCCCGGTGTAAAGCATGGCGAAACCATCAACGCGCATCAGCGGCGTAACGTCCATAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGCCTTGTAGTAGTAACGGAATTAACAAACCACTGCTGGATATTGC  >  NC_000913/2391270‑2391561

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