Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,399,878 C→G 5.6% M96I (ATG→ATC) ‡ nuoG ← NADH:ubiquinone oxidoreductase, chain G

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,399,8780CG 5.6% 200.8 / 4.1 71M96I (ATG→ATC) ‡nuoGNADH:ubiquinone oxidoreductase, chain G
Reads supporting (aligned to +/- strand):  ref base C (35/32);  new base G (2/2);  total (37/34)
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.16e-01

GAGATGAATGGCCCCAAATCCTGATTACGGTGGGTACGTTTGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCACCAGGCGACCACGCGTGTCTTCCGCGTTTTGGTATTGCTTCACCGCACAC  >  NC_000913/2399732‑2400025
                                                                                                                                                  |                                                                                                                                                   
gagaTGAATGGCCCCAAATCCTCATTACGGTGGGAACGTTTGGTGAACCGTTAGCCACGGAAGCTGTGACCGCTCATCACAGTCATAACCTGAACATGGCAGTTACCGCCCGCTTCACATCCCGGCCCGTCCTGCGGTTCGTTGGtcat                                                                                                                                                   >  2:674846/1‑149 (MQ=255)
   aTGAATGGCCCCAAATCCTGATTACGGTGGGTACGTTTGGTGAATCGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCTCATACCGGACAGTCGTGCGGGTGGTTGGTCATCaa                                                                                                                                                <  2:319423/149‑1 (MQ=255)
   aTGAATGGCCCCAAATCCTGATTACGGTGGGTACGTTTGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCATCTTCACATACCGGACAGTCGTGCGGTTGTTTGGTCATCaa                                                                                                                                                >  2:699830/1‑149 (MQ=255)
      aaTGGCCCCAAATCCTGATTACGGTGGGTACGTTTGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAAccc                                                                                                                                             >  2:1045827/1‑148 (MQ=255)
         ggCCCCAAATCCTGATTACGGTGGGTACGTTTGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTc                                                                                                                                          >  1:45090/1‑149 (MQ=255)
            cccAAATCCTGATTACGGTGGGTACGTTTGGTGAAACGGTAGCGACGGAAGCTGGGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTCACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGAc                                                                                                                                       <  2:689042/149‑1 (MQ=255)
            cccAAATCCTGAATACGGTGGGTACGTTTGGTGAAACGGTAGCGACGGAAGCTGTGTCAGGTCATCACAGTCATATACTGAAGATGGCAGTTACCGCCCTCTTCACAAACAGGACAGTCGTCCGGCTGCTTGGTCATCTCCCCCTCGAc                                                                                                                                       >  2:519434/1‑149 (MQ=255)
             ccAAATCCTGATTACGGTGGGTACGTTTGGTGAAACGGTAGCGACGGAAGCTGTGGCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGAGTGGTTGATGACCAACCACCCgcac                                                                                                                                      <  1:558962/149‑4 (MQ=255)
             ccAAATCCTGATTACGGGGGGTACGGTTGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGGCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGATGACCAACCACCCgcac                                                                                                                                      <  2:1004993/149‑4 (MQ=255)
              cAAATCCTGATAACGGTGGGTACGTTTGGTCAAACGGAAGCGACTGAAGCTGTGTCCGGTCATCACACTCATATACTGAAGATGGCAGATACCGCCCACTTCACATACCGCACACTCGTGCGCGATCTTGGTCATCAACCACTCGAccc                                                                                                                                     >  2:600173/1‑148 (MQ=255)
                aaTCCTGATTTCGGTGGGTACGTTTGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACg                                                                                                                                   <  2:943314/149‑1 (MQ=255)
                aaTCCTGATTACGGTGGGTACGTTTGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACg                                                                                                                                   <  2:938786/149‑1 (MQ=255)
                    cTGATTACGGTGGGTACGTTTGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGGTGGTCATCAACCACTCGACCACGCttt                                                                                                                               <  2:621418/149‑1 (MQ=255)
                     tGATTACGGTGGGTACGTTTGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTc                                                                                                                              <  1:248833/149‑1 (MQ=255)
                     tGATTACGGTGGGTACGTTTGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTGCACATACAGGACAGTCGTGCGGGTGGATGGACATCAACCCCTCGACCCCGCTTTc                                                                                                                              >  2:1105133/1‑149 (MQ=255)
                           cGGTGGGTACGTTTGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGaa                                                                                                                        <  1:485315/149‑1 (MQ=255)
                            ggTGGGTAAGTTTGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAAc                                                                                                                       <  2:828545/149‑1 (MQ=255)
                              tGGGTACGTTTGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTg                                                                                                                     >  1:743929/1‑149 (MQ=255)
                              ggggTACGTTTGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTg                                                                                                                     <  1:1075843/148‑1 (MQ=255)
                                ggTACGTTTGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATACTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGtt                                                                                                                   >  1:967388/1‑149 (MQ=255)
                                  tACGTTTGGTGAAACGGTAGCGACGGCAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTc                                                                                                                 <  1:1149238/149‑1 (MQ=255)
                                    cgctctctGAAACGGTAGCGACGGAAGCTGTCGCCTGTCCTCACAGTCATCTCCTGCAGAAGGCAGTTACCGCCCTCTTCACATACCGTCCAGGCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGc                                                                                                               <  2:820685/142‑1 (MQ=255)
                                    cGTTTGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGGCTTCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGc                                                                                                               <  1:192462/149‑1 (MQ=255)
                                     gTTTGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGct                                                                                                              >  2:275438/1‑149 (MQ=255)
                                           tGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTtcg                                                                                                        >  2:406830/1‑149 (MQ=255)
                                            gAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTtcgt                                                                                                       <  2:1141868/149‑1 (MQ=255)
                                            gAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCACTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTtcgt                                                                                                       >  1:1033735/1‑149 (MQ=255)
                                            gAAACGGTAGCGACGAAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTtcgt                                                                                                       <  1:306865/149‑1 (MQ=255)
                                              accGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTAACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGCCCACGCTTTCACGGAACTGTTTCGCTTCTtcgtcg                                                                                                     >  1:933136/3‑149 (MQ=255)
                                              aaCGGTAGCGACGGAAGCTGTGTACTGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACCTACCGGACAGTCGTGCGTGTAGTTGGTCATCCAACACTCGCCCAATCTTTCAACGAACTGTTTCCCTTCTgcgacg                                                                                                     >  2:1165326/1‑149 (MQ=255)
                                                   tAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAAt                                                                                                <  1:730890/149‑1 (MQ=255)
                                                      cGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCGTCTTCGTCGTCAATGGa                                                                                             >  1:191833/1‑149 (MQ=255)
                                                            aaGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATaaa                                                                                       <  2:899849/149‑1 (MQ=255)
                                                             aGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAg                                                                                      <  1:1045827/149‑1 (MQ=255)
                                                             aGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAg                                                                                      <  1:1158362/149‑1 (MQ=255)
                                                              gCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAgg                                                                                     >  1:624734/1‑149 (MQ=255)
                                                              cctgggcccGGTCCTCACACTCATACCCTGAAGCTGCCAGCACCCGCCCTCTTCACATCCCCGACAGGCGTGCGGGGGGTAGGCCATCACCCACTCGACCACGCTTTCACGGAACTGCTTCGCTTCTGCGTCGTCAATGGAAATAAAgg                                                                                     <  2:868105/142‑1 (MQ=255)
                                                                  tgtCCGGTCATCCCAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACCTACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCAAGGAACTGTTGCGCTTATTAGTCGTCAATGGAAATAAAGGTgac                                                                                 >  1:908069/1‑147 (MQ=255)
                                                                  tgtCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACGGCCCTCATCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGTACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGcc                                                                                 <  1:600173/149‑1 (MQ=255)
                                                                         tcatcaCAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATCCCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGaa                                                                          >  2:1058009/1‑149 (MQ=255)
                                                                          catcaCAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCCCTCGACCACGCTTTCACGGAACGGTTTCGCTTCTTCGGCGGCAATGGAACTAAAGGTGCCCTCGGaac                                                                         >  1:781761/1‑148 (MQ=255)
                                                                                aGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTAGACCCCGCTTTCACGGACCTGTTTCGCTTCTTCGGCGTCAATGGCAATAAAGGGGCCATCGGAAGCCgggg                                                                   >  2:802180/1‑147 (MQ=255)
                                                                                  tCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCAGTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTc                                                                 <  1:956912/149‑1 (MQ=255)
                                                                                  tCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCATTGGAAATAAAGGTGCCATCGGAAGCCGGTGTc                                                                 >  1:318989/1‑149 (MQ=255)
                                                                                      atCCTGAAGATGGCAGTTACCGCCATCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATAc                                                             <  1:344352/149‑1 (MQ=255)
                                                                                        tgtGAAGAGGGCGGTAACTGCCATCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAg                                                           >  1:626863/3‑149 (MQ=255)
                                                                                        ccTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAg                                                           <  2:923320/149‑1 (MQ=255)
                                                                                         cTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAGCCACTCGACCCCGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAAGGGAAATTAAGGTGCCAGCGGAAGCCGGTGTGATACAgg                                                          >  1:148452/1‑149 (MQ=255)
                                                                                          tGAAGATGGCAGTTACCGCCCTCTTCACATAACGGACAGTCGTGCGGGTGGTTGGTCATCCACCACTCGACCACGCGTTCACGGAACGGGTTCGCTTATTCGTCTTCAATGGAAATAAAGGTGCGATCGGAAGCCGGTGTCATACAGGa                                                         >  2:867727/1‑149 (MQ=255)
                                                                                          tGAAGAGGGCGGTAACTGCCATCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGa                                                         >  1:785006/1‑149 (MQ=255)
                                                                                             aGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACAt                                                      <  1:515705/149‑1 (MQ=255)
                                                                                                 ggCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCAcc                                                  >  2:917281/1‑149 (MQ=255)
                                                                                                   cAGTGACCGCCCTCTTCACATACCGGACAGTCGTGCGGGGGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACATGACATCACCAg                                                <  2:839673/149‑1 (MQ=255)
                                                                                                         ccGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCACCAGGCGAcc                                          >  1:574956/1‑149 (MQ=255)
                                                                                                         ccGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACACCGCTTTCACGGAACTGTTTCGATTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCACCAGGCGAcc                                          >  2:18057/1‑149 (MQ=255)
                                                                                                          cGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCACCAGGCGACCa                                         >  1:320514/1‑149 (MQ=255)
                                                                                                            cccTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTATCACGGAACTGTTTCGCTTCTTCGTCGTCAAAGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGCCATCACCAGGCGACCAcg                                       >  2:358363/1‑149 (MQ=255)
                                                                                                             ccTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCACCAGGCGACCAcgc                                      <  2:967388/149‑1 (MQ=255)
                                                                                                                ctTCACAAACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCAGGTGTCATACAGGACATCACCAGGCGACCACGCgtg                                   >  2:921056/1‑149 (MQ=255)
                                                                                                                          cGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCACCAGGCGACCACGCGTGTCTTCCGCGt                         <  2:785006/149‑1 (MQ=255)
                                                                                                                              cAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCACCAGGCGACCACGCGTGTCTTCCGCGTTTTg                     <  1:406830/149‑1 (MQ=255)
                                                                                                                                 tCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCACCAGGCGACCACGCGTGTCTTCCGCGTTTTGGTa                  >  2:47975/1‑149 (MQ=255)
                                                                                                                                 tCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCACCAGGCGACCACGCGTGTCTTCCGCGTTTTGGTa                  >  1:881972/1‑149 (MQ=255)
                                                                                                                                  tggtcgagtggtTGATGACCAACCACCCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCACCAGGCGACCACGCGTGTCTTCCGCGTTTTGGTAt                 >  1:334008/8‑149 (MQ=255)
                                                                                                                                   gTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCACCAGGCGACCACGCGTGTCTTCCGCGTTTTGGTAtt                >  2:354929/1‑149 (MQ=255)
                                                                                                                                   gTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCACCAGGCGACCACGCGTGTCTTCCGCGTTTTGGTAtt                <  2:1033735/149‑1 (MQ=255)
                                                                                                                                     gCGGGTGGGTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCTTCGGCAATGGAAATATAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCACCAGGCTACCACGCGTGTCTTCCGCGTTTTGGTATTGc              <  1:634312/149‑1 (MQ=255)
                                                                                                                                        ggtggaAGGTCATCAACCACTCGCCCCCGCTTTCAGGGAACTGTTTCGCTTCTTCGTCGTCAATGCAAATAAAGGTGCCATCGCAAGCCGGATTCAAACAGGACATCACCAGGCGACCACGCGTGTCTTCCGCGTTTTGGTAATGCTTc           >  2:685875/1‑149 (MQ=255)
                                                                                                                                         gtggtTGATGACCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGATATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCACCAGGCTACCACGCGTGTATTACGCGTTTTGGTATTGCTTcc          >  2:677094/1‑148 (MQ=255)
                                                                                                                                          tggttggtCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCACCAGGCGACCACGCGTGTCTTCCGCGTTTTGGTATTGCTTCAc         <  2:624734/149‑1 (MQ=255)
                                                                                                                                           ggttggtCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCACCAGGCGACCACGCGTGTCTTCCGCGTTTTGGTATTGCTTCAcc        <  2:45090/149‑1 (MQ=255)
                                                                                                                                               ggtCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCACCAGGCGACCACGCGTGTCTTCCGCGTTTTGGTATTGCTTCACCGcac    >  2:709161/1‑149 (MQ=255)
                                                                                                                                                 tcatcaACCACTCGACCACGCATTCACGCAAATCTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGAGCCAACGGAAGCCGGTGTCATACAGGACATCACCAGGCGACCACGCGTGACTTCCGCGTTTTGGAATTCCTTCACCGcaccc  >  2:641101/1‑147 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GAGATGAATGGCCCCAAATCCTGATTACGGTGGGTACGTTTGGTGAAACGGTAGCGACGGAAGCTGTGTCCGGTCATCACAGTCATATCCTGAAGATGGCAGTTACCGCCCTCTTCACATACCGGACAGTCGTGCGGGTGGTTGGTCATCAACCACTCGACCACGCTTTCACGGAACTGTTTCGCTTCTTCGTCGTCAATGGAAATAAAGGTGCCATCGGAAGCCGGTGTCATACAGGACATCACCAGGCGACCACGCGTGTCTTCCGCGTTTTGGTATTGCTTCACCGCACAC  >  NC_000913/2399732‑2400025

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