Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,400,684 T→A 7.7% Q291L (CAG→CTG)  nuoF ← NADH:ubiquinone oxidoreductase, chain F

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,400,6840TA 7.7% 304.1 / 11.8 91Q291L (CAG→CTG) nuoFNADH:ubiquinone oxidoreductase, chain F
Reads supporting (aligned to +/- strand):  ref base T (42/42);  new base A (4/3);  total (46/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

ACCAGCGACACCATGTTGATCTCATGGTCAACCGCCATCGCCAGCGCCGTACCCAGACGGCTGCCCGCTTTACCGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGCCGAACGGCAGTTCCCACAGTCCCGGATTTTTCACCCGACCGGAGAAGCCCATCAGCTTGGTGCCAGCATCTTTACT  >  NC_000913/2400540‑2400832
                                                                                                                                                |                                                                                                                                                    
aCCAGCGACACCATGTTGATCTCATGGTCAACCGCCATCGCCAGCGCCGTACCCAGACGGCTGCCCGCTTTACCGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCa                                                                                                                                                  >  2:16396/1‑149 (MQ=255)
aCCAGCGACACCATGTTGATCTCATGGTCAACCGCCATCGCCAGCGCCGTACCCAGACGGCTGCCCGCTTTACCGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCa                                                                                                                                                  <  1:16396/149‑1 (MQ=255)
aCCAGCGACACCATGTTGATCTCATGGTCAACCGCCATCGCCAGCGCCGTACCCAGACGGCTGCCCGCTTTACCGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCa                                                                                                                                                  >  1:786566/1‑149 (MQ=255)
 ccAGCGACACCATGTTGATCTCATGGTCAACCGCCATCGCCAGCGCCGTACCCAGACGGCTGCCCGCTTTACCGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAg                                                                                                                                                 >  1:501178/1‑149 (MQ=255)
   aGCGACACCATGTTGATCTCATGGTCAACCGCCATCGCCAGCGCCGTACCCAGACGGCTGCCCGCTTTACCGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGc                                                                                                                                               <  2:477634/149‑1 (MQ=255)
     cGACACCATGTTGATCTCATGGTCAACCGCCATCGCCAGCGCCGTACCCAGACGGCTGCCCGCTTTACCGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCtt                                                                                                                                             <  3:369639/149‑1 (MQ=255)
        cacCATGTTGATCTCATGGTCAACCGCCATCGCCAGCGCCGTACCCAGACGGCTGCCCGCTTTACCGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCAGCCAGGcggcgc                                                                                                                                          <  2:305392/149‑6 (MQ=255)
        cacCATGTTGATCTCATGGTCAACCGCAATCGCCAGCGCCGTACCCAGACGGCTGCCCGCTTTACCGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCAGCCAGGcggcgc                                                                                                                                          <  4:312343/149‑6 (MQ=255)
         acCATGTTGATCTCATGGTCAACCGCCATCGCCAGCGCCGTACCCAGACGGCTGCCCGCTTTACCGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTaaa                                                                                                                                         <  1:295053/149‑1 (MQ=255)
         acCATGTTGATCTCATGGTCAACCGCCATCGCCAGCGCCGTACCCAGACGGCTGCCCGCTTTACCGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTaaa                                                                                                                                         <  2:269305/149‑1 (MQ=255)
         acCATGTTGATCTCATGGTCAACCGCCATCGCCAGCGCCGTACCCAGACGGCTGCCCGCTTTACCGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTaaa                                                                                                                                         >  4:665174/1‑149 (MQ=255)
         acCATGTTGATCTCATGGTCAACCGCCATCGCCAGCGCCGTACCCAGACGGCTGCCCGCTTTACCGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTaaa                                                                                                                                         >  4:16055/1‑149 (MQ=255)
         acCATGTTGATCTCATGGTCAACCGCCATCGCCAGCGCCGTACCCAGACGGCTGCCCGCTTTACCGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTaaa                                                                                                                                         >  2:826188/1‑149 (MQ=255)
         acCATGTTGATCTCATGGTCAACCGCCATCGCCAGCGCCGTACCCAGACGGCTGCCCGCTTTACCGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTaaa                                                                                                                                         <  3:65267/149‑1 (MQ=255)
          ccATGTTGATCTCATGGTCAACCGCCATCGCCAGCGCCGTACCCAGACGGCTGCCCGCTTTACCGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAAt                                                                                                                                        >  4:454028/1‑149 (MQ=255)
          ccATGTTGATCTCATGGTCAACCGCCATCGCCAGCGCCGTACCCAGACGGCTGCCCGCTTTACCGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCACCGCCTGGCAGCCAGGcggcgcgg                                                                                                                                        <  2:655816/149‑8 (MQ=255)
            aTGTTGATCTCATGGTCAACCGCCATCGCCAGCGCCGTACCCAGACGGCTGCCCGCTTTACCGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAAttt                                                                                                                                      <  1:168035/149‑1 (MQ=255)
              gTTGATCTCATGGTCAACCGCCATCGCCAGCGCCGTACCCAGACGGCTGCCCGCTTTACCGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCa                                                                                                                                    >  4:383596/1‑149 (MQ=255)
               ttGATCTCATGGTCAACCGCCATCGCCAGCGCCGTACCCAGACGGCTGCCCGCTTTACCGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAg                                                                                                                                   >  4:315977/1‑149 (MQ=255)
                tGATCTCATGGTCAACCGCCATCGCCAGCGCCGTACCCAGACGGCTGCCCGCTTTACCGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGa                                                                                                                                  <  3:783928/149‑1 (MQ=255)
                 gATCTGATGGCCAACCGCAACAGCCAGCGCCGTAGCCAGATGGCTGCACGCCTTGCCGAGACTTCCGAATTCCATCGGCAGATCAAGGTTCCCTTCGATCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAAGCTTTAAATTTCacac                                                                                                                                 <  2:401380/149‑1 (MQ=255)
                 gATCTCATGGTCAACCGCCATCGCCAGCGCCGTACCCAGACGGCTGCCCGCTTTACCGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCa                                                                                                                                    >  1:12512/1‑146 (MQ=255)
                 gATCTCATGGTCAACCGCCATCGCCAGCGCCGTACCCAGACGGCTGCCCGCTTTACCGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCa                                                                                                                                    <  2:12512/146‑1 (MQ=255)
                     tcATGGTCAACCGCCATCGCCAGCGCCGTACCCAGACGGCTGCCCGCCTTACCGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATc                                                                                                                             >  4:464604/1‑149 (MQ=255)
                         ggTCAACCGCCATCGCCAGCGCCGTACCCAGACGGCTGCCCGCTTTACCGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTTGGTCAGGAAGTCAGTCCCCGCGCCGGCTGGCTGCCAGGCTTTAAATTTCCGACCATCACGc                                                                                                                         >  4:516282/1‑149 (MQ=255)
                         ggTCAACCGCCATCGCCAGCGCCGTACCCAGACGGCTGCCCGCTTTACCGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGc                                                                                                                         >  3:117534/1‑149 (MQ=255)
                          gTCAACCGCCATCGCCAGCGCCGTACCCAGACGGCTGCCCGCTTTACCGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCa                                                                                                                        >  3:413636/1‑149 (MQ=255)
                           tccaCCGCCATCGCCAGCGCCGTACCCAGACGGCTGCCCGCTTTACCGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCAt                                                                                                                       <  4:522040/146‑1 (MQ=255)
                             aaCCGCCATCGCCAGCGCCGTACCCAGACGGCTGCCCGCTTTACCGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATac                                                                                                                     >  1:303894/1‑149 (MQ=255)
                                 gCCATCGCCAGCGCCGTACCCAGACGGCTGCCCGCTTTACCGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATaccacc                                                                                                                 >  2:233762/1‑149 (MQ=255)
                                          aGCGCCGTACCCAGACGGCTGCCCGCTTTACCGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATc                                                                                                        >  4:622242/1‑149 (MQ=255)
                                             gcCGTACCCAGACGGCTGCCCGCTTTACCGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTc                                                                                                     <  1:494054/149‑1 (MQ=255)
                                               cGTACCCAGACGGCTGCCCGCTTTACCGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGa                                                                                                   <  3:16055/149‑1 (MQ=255)
                                               cGTACCCAGACGGCTGCCCGCTTTACCGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGa                                                                                                   <  2:501178/149‑1 (MQ=255)
                                               cGTACCCAGACGGCTGCCCGCTTTACCGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGa                                                                                                   <  2:344402/149‑1 (MQ=255)
                                               cGTACCCAGACGGCTGCCCGCTTTACCGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGa                                                                                                   <  2:293254/149‑1 (MQ=255)
                                               cGTACCCAGACGGCTGCCCGCTTTACCGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGa                                                                                                   <  2:121585/149‑1 (MQ=255)
                                                gTACCCAGACGGCTGCCCGCTTTACCGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAg                                                                                                  <  1:565235/149‑1 (MQ=255)
                                                gTACCCAGACGGCTGCCCGCTTTACCGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAg                                                                                                  <  4:56126/149‑1 (MQ=255)
                                                  aCCCAGACGGCTGCCCGCTTTACCGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCACCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGa                                                                                                >  3:307256/1‑149 (MQ=255)
                                                    ccAGACGGCTGCCCGCTTTACCGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGAtc                                                                                              >  1:374032/1‑149 (MQ=255)
                                                             tGCCCGCTTTACCGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCg                                                                                     >  2:803658/1‑149 (MQ=255)
                                                                  gCTTTACCGATACTTTCTAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCCAGTATCTCGCGTGCgggggt                                                                                >  2:791962/1‑149 (MQ=255)
                                                                  gCTTTACCGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCggtggt                                                                                <  3:34444/149‑1 (MQ=255)
                                                                  gCTTTACCGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCggtggt                                                                                >  4:341726/1‑149 (MQ=255)
                                                                    tttACCGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGc                                                                              <  2:252948/149‑1 (MQ=255)
                                                                         cGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGATCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGCCGAAc                                                                         >  2:250081/1‑149 (MQ=255)
                                                                         cGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGCCGAAc                                                                         >  4:265051/1‑149 (MQ=255)
                                                                         cGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGCCGAAc                                                                         >  2:218859/1‑149 (MQ=255)
                                                                             aCTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGCCGAACGGCa                                                                     <  2:764243/149‑1 (MQ=255)
                                                                             aCTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGCCGAACGGCa                                                                     >  3:102737/1‑149 (MQ=255)
                                                                             aCTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGCCGAACGGCa                                                                     <  1:484768/149‑1 (MQ=255)
                                                                               tttCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATac                                                                                                                     >  3:712828/1‑99 (MQ=255)
                                                                               tttCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATac                                                                                                                     <  4:712828/99‑1 (MQ=255)
                                                                               tttCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGCCGAACGGCAGt                                                                   >  4:341842/1‑149 (MQ=255)
                                                                               tttCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGCCGAACGGCAGt                                                                   >  1:152905/1‑149 (MQ=255)
                                                                               tttCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGCCGAACGGCAGt                                                                   >  1:774580/1‑149 (MQ=255)
                                                                                    aaTTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGCCGAACGGCAGTTCCca                                                              <  4:400753/149‑1 (MQ=255)
                                                                                    aaTTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGCCGAACGGCAGTTCCca                                                              <  4:360082/149‑1 (MQ=255)
                                                                                      ttCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGCCGAACGGCAGTTCCcaca                                                            >  2:560711/1‑149 (MQ=255)
                                                                                         cATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGCCGAACGGCAGTTCCCACAGTc                                                         >  3:561954/1‑149 (MQ=255)
                                                                                           tCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGCCGAACGGCAGTTCCCACAGTccc                                                       <  1:261641/149‑1 (MQ=255)
                                                                                              gCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGCCGAACGGCAGTTCCCACAGTCCCGGa                                                    >  3:215580/1‑149 (MQ=255)
                                                                                              gCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGCCGAACGGCAGTTCCCACAGTCCCGGa                                                    <  4:561954/149‑1 (MQ=255)
                                                                                               cAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGCCGAACGGCAGTTCCCACAGTCCCGGAt                                                   <  4:413636/149‑1 (MQ=255)
                                                                                                 gATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGCCGAACGGCAGTTCCCACAGTCCCGGAttt                                                 <  4:665750/149‑1 (MQ=255)
                                                                                                  aTCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGCCGAACGGCAGTTCCCACAGTCCCGGAtttt                                                >  3:506863/1‑149 (MQ=255)
                                                                                                  aTCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGCCGAACGGCAGTTCCCACAGTCCCGGAtttt                                                <  1:332189/149‑1 (MQ=255)
                                                                                                     aaGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGCCGAACGGCAGTTCCCACAGTCCCGGATTTTTCa                                             <  2:470756/149‑1 (MQ=255)
                                                                                                     aaGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGCCGAACGGCAGTTCCCACAGTCCCGGATTTTTCa                                             <  2:429677/149‑1 (MQ=255)
                                                                                                        gTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGCCGAACGGCAGTTCCCACAGTCCCGGATTTTTCAccc                                          <  4:215580/149‑1 (MQ=255)
                                                                                                         tGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGCCGAACGGCAGTTCCCACAGTCCCGGATTTTTCACCCg                                         >  4:561343/1‑149 (MQ=255)
                                                                                                              ttCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGCCGAACGGCAGTTCCCACAGTCCCGGATTTTTCACCCGACCgg                                    >  4:226948/1‑149 (MQ=255)
                                                                                                              ttCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGCCGAACGGCAGTTCCCACAGTCCCGGATTTTTCACCCGACCgg                                    >  3:164875/1‑149 (MQ=255)
                                                                                                               tCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGCCGAACGGCAGTTCCCACAGTCCCGGATTTTTCACCCGACCGga                                   <  4:134302/149‑1 (MQ=255)
                                                                                                                     aGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGCCGAACGGCAGTTCCCACAGTCCCGGATTTTTCACCCGACCGGAGAAGcc                             >  1:474107/1‑149 (MQ=255)
                                                                                                                      ggAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGCCGAACGGCAGTTCCCACAGTCCCGGATTTTTCACCCGACCGGAGAAGccc                            >  2:786432/1‑149 (MQ=255)
                                                                                                                       gAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGCCGAACGGCAGTTccc                                                               >  4:701709/1‑113 (MQ=255)
                                                                                                                       gAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGCCGAACGGCAGTTccc                                                               <  3:701709/113‑1 (MQ=255)
                                                                                                                           tcagtcCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGCCGAACGGCAGTTCCCACAGTCCCGGATTTTTCACCCGACCGGAGAAGCCCATCAg                       >  1:825125/1‑149 (MQ=255)
                                                                                                                             agtcCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGCCGAACGGCAGTTCCCACAGTCCCGGATTTTTCACCCGACCGGAGAAGCCCATCAgat                     <  2:819599/149‑3 (MQ=255)
                                                                                                                                ccccGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGCCGAACGGCAGTTCCCACAGTCCCGGATTTTTCACCCGACCGGAGAAGCCCATCAGCTTgg                  <  1:299795/149‑1 (MQ=255)
                                                                                                                                 cccGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGCCGAACGGCAGTTCCCACAGTCCCGGATTTTTCACCCGACCGGAGAAGCCCATCAGCTTGGt                 >  3:355843/1‑149 (MQ=255)
                                                                                                                                   tttaaagccTGGCAGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGCCGAACGGCAGTTCCCACAGTCCCGGATTTTTCACCCGACCGGAGAAGCCCATCAGCTTGGTGc               >  1:415288/7‑149 (MQ=255)
                                                                                                                                    ttaaagccTGGCAGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGCCGAACGGCAGTTCCCACAGTCCCGGCTTTTTCACCCGACCGGAGAAGCCCATCAGCTTGGTGcc              >  4:70482/6‑149 (MQ=255)
                                                                                                                                     taaagccTGGCAGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGCCGAACGGCAGTTCCCACAGTCCCGGATTTTTCACCCGACCGGAGAAGCCCATCAGCTTGGTGCCa             >  1:229886/5‑149 (MQ=255)
                                                                                                                                     cgccgcCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGCCGAACGGCAGTTCCCACAGTCCCGGATTTTTCACCCGACCGGAGAAGCCCATCAGCTTGGTGCCa             <  4:555774/149‑1 (MQ=255)
                                                                                                                                      aaagccTGGCAGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGCCGAACGGCAGTTCCCACAGTCCCGGATTTTTCACCCGACCGGAGAAGCCCATCAGCTTGGTGCCAg            >  3:768382/4‑149 (MQ=255)
                                                                                                                                         gccTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGCCGAACGGCAGTTCCCACAGTCCCGGATTTTTCACCCGACCGGAGAAGCCCATCAGCTTGGTGCCAGCAt         <  3:315977/149‑1 (MQ=255)
                                                                                                                                              gCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGCCGAACGGCAGTTCCCACAGTCCCGGATTTTTCACCCGACCGGAGAAGCCCATCAGCTTGGTGCCAGCATCTTTa    <  2:285458/149‑1 (MQ=255)
                                                                                                                                               cTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGCCGAACGGCAGTTCCCACAGTCCCGGATTTTTCACCCGACCGGAGAAGCCCATCAGCTTGGTGCCAGCCTCTTTAc   <  2:774580/149‑1 (MQ=255)
                                                                                                                                                tGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGCCGAACGGCAGTTCCCACAGTCCCGGATTTTTCACCCGACCGGAGAAGCCCATCAGCTTGGTGCCAGCATCTTTACt  <  3:244271/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
ACCAGCGACACCATGTTGATCTCATGGTCAACCGCCATCGCCAGCGCCGTACCCAGACGGCTGCCCGCTTTACCGATACTTTCGAATTCCATCGGCAGATCAAGGTGCGCTTCGGTCAGGAAGTCAGTCCCCGCGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGCCGAACGGCAGTTCCCACAGTCCCGGATTTTTCACCCGACCGGAGAAGCCCATCAGCTTGGTGCCAGCATCTTTACT  >  NC_000913/2400540‑2400832

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