Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,299,056 T→C 7.0% Q104R (CAG→CGG)  fdhF ← formate dehydrogenase‑H, selenopolypeptide subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,299,0560TC 7.0% 183.7 / 5.4 71Q104R (CAG→CGG) fdhFformate dehydrogenase‑H, selenopolypeptide subunit
Reads supporting (aligned to +/- strand):  ref base T (32/34);  new base C (2/3);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.40e-01

CGACCGATTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGCGCGGGGTCAGGATCTGG  >  NC_000913/4298910‑4299202
                                                                                                                                                  |                                                                                                                                                  
cgaccgaTTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTgg                                                                                                                                                  >  2:100384/1‑149 (MQ=255)
  accgaTTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTTTTGGTAGCAATAACGGCGCGCGCAAATTTTTGCATTGCAGAGTTGGTTTCGTTCCCCGTACCGAGCGAGGAGCCGGTCGTCTGGAt                                                                                                                                                >  1:576766/1‑149 (MQ=255)
   ccgaTTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATg                                                                                                                                               >  1:831209/1‑149 (MQ=255)
     gaTTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCCGGATgcc                                                                                                                                             <  1:341875/149‑3 (MQ=255)
       ttGGTGCAGACCTGCAAACTATGGGCCGTGTCAGACACGAGCGAAGCAGTTAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCAt                                                                                                                                           <  1:560110/149‑1 (MQ=255)
            gCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCCCGCGAGGAGCCGGTCGTCTGGATGGCATCCGGa                                                                                                                                      >  1:520034/1‑149 (MQ=255)
                aCCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGt                                                                                                                                  <  1:924602/149‑1 (MQ=255)
                aCCGGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCCGGATGCCATCCAGAcgac                                                                                                                                  <  2:365556/149‑4 (MQ=255)
                 ccTGCTTCCGATGGGCCGGGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTa                                                                                                                                 <  2:598765/149‑1 (MQ=255)
                    gCAACCGATGGGCCGTGTCAGCCACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTGACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGGATCCGGGCCGttcgt                                                                                                                              >  1:500475/1‑145 (MQ=255)
                        ccGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGAGGGCAGCCGGACCGAACTTcgct                                                                                                                          >  2:500573/1‑149 (MQ=255)
                           aTGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATATTTGGTTTCGTTACCCGTAACACGCGAGGCGCCGGTCGTCTGGATGGCATCCGGGCCGGACTTCTCATTg                                                                                                                       >  2:986768/1‑149 (MQ=255)
                            tGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGa                                                                                                                      <  2:315463/149‑1 (MQ=255)
                              ggCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATg                                                                                                                    <  1:379653/149‑1 (MQ=255)
                                ccGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGCCCGTACTTCTCTTTGATGgc                                                                                                                  >  1:166308/1‑149 (MQ=255)
                                    gtCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCCTCCGGACCGTACTTCTCTTTGATGGCGCTc                                                                                                              >  2:962530/1‑149 (MQ=255)
                                     tCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCa                                                                                                             >  1:326818/1‑149 (MQ=255)
                                     tCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAAGAACGGCGCGCGCAATTTTGTGCATTACATACTTGGTTTCGTTACCCGAACCACGCGAGG‑TCCGGTCGTCCGGAAGGCAGCCGGAGGGTACATGTGATTGCTGGCGCTCAg                                                                                                            >  2:226582/1‑149 (MQ=255)
                                     gcAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGCCCGTACTTCTCTTTGATGGCGCTCa                                                                                                             >  2:639595/2‑149 (MQ=255)
                                      cAGACAAGAGATCAGCATTCAAAGTTATTTGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATTGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAg                                                                                                            <  1:205275/149‑1 (MQ=255)
                                       aGACACGAGCGCAGCAGTCAACCTTATTGGTACCAAAAACGGCGCGCGCAAATTTTTCCATTACATAGTTGGTTGCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAgg                                                                                                           <  2:417274/149‑1 (MQ=255)
                                         acacGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGgcg                                                                                                         <  1:100384/149‑1 (MQ=255)
                                             gAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGAGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGAAATTCTCTTTGATTGCGCTCAGGCGCTCg                                                                                                     >  1:610763/1‑149 (MQ=255)
                                                    gcagTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTAGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCACCGt                                                                                              >  1:335935/1‑149 (MQ=255)
                                                      agTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTaa                                                                                            <  1:115478/149‑1 (MQ=255)
                                                         cAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGGTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTc                                                                                         <  1:32556/149‑1 (MQ=255)
                                                          aaCGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCa                                                                                        >  2:297982/1‑149 (MQ=255)
                                                           aCGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAg                                                                                       <  2:326818/149‑1 (MQ=255)
                                                              ttATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGc                                                                                    <  1:128263/149‑1 (MQ=255)
                                                                 ttGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTc                                                                                 <  2:1013910/149‑1 (MQ=255)
                                                                    gTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGCCCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGGAATTCAGTGCctcctc                                                                              >  2:912806/1‑149 (MQ=255)
                                                                       ccAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCa                                                                           <  1:622313/149‑1 (MQ=255)
                                                                        cAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAg                                                                          >  1:88903/1‑149 (MQ=255)
                                                                        cAATAACGGCGCGCGCAAATTTTTGCATAACCTAGTTGGTTTCGTTACCCGTTCCCCGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCCGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAg                                                                          >  2:12315/1‑149 (MQ=255)
                                                                         aaTAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGGGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCGCAGGCGCTCGGCAACGTAATTAAGTGCCTCATCCCAgg                                                                         >  1:173806/1‑149 (MQ=255)
                                                                          aTAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATCGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCCGGa                                                                        >  1:660162/1‑149 (MQ=255)
                                                                          aTAACGGCGCGCGCAAATTTTTGCATTACATAGGTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGa                                                                        <  2:316075/149‑1 (MQ=255)
                                                                             aCGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAAc                                                                     >  2:9691/1‑149 (MQ=255)
                                                                                     cgcAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGa                                                             <  1:620135/149‑1 (MQ=255)
                                                                                      gcAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGGCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAg                                                            <  2:226087/149‑1 (MQ=255)
                                                                                        aaaTTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCGTCATCCCCGGAAACAGGTTCGAGtt                                                          >  2:985492/1‑149 (MQ=255)
                                                                                             tttGCATTACAGAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTgccg                                                     >  1:371878/1‑149 (MQ=255)
                                                                                              ttGCATTACATAGTTGGTTTCGTTTCCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATGCGGACCGTACTTCTCTTTCATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTgccgc                                                    <  2:255767/149‑1 (MQ=255)
                                                                                               tGCATTACATAGTTGGTTTCGATACCCGTACCACCCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACCTAATTCCGTGCCTCATCCCACGAAACAGGTGCGAGTTTgccgcc                                                   <  2:576766/149‑1 (MQ=255)
                                                                                                     aCATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGATACAGGTTCGAGTTTGCCGCCACGCTg                                             >  2:289776/1‑149 (MQ=255)
                                                                                                      cATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTgg                                            <  2:512160/149‑1 (MQ=255)
                                                                                                           ttGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCACGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACg                                       >  2:371317/1‑149 (MQ=255)
                                                                                                             ggTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGa                                     <  1:923952/149‑1 (MQ=255)
                                                                                                             ggTATCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCGCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGGTCGAGTTTGCCGCCACGCTGGCGACGGa                                     <  1:1014969/149‑1 (MQ=255)
                                                                                                               tttCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATc                                   >  1:50561/1‑149 (MQ=255)
                                                                                                               tttCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGGTCGAGTTTGCCGCCACGCGGGCGACGGCTc                                   >  1:929738/1‑149 (MQ=255)
                                                                                                                 tcgctccccGTCCCACGCGAGGAGCCGGGCGGCTGGATGGCATCCGGCCCGTACTTCGCTTTGCGGGCGCTCCGGCGCTCGGCACCGTAATTCAGTGCCTCATCCCCGGAAACCGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCAt                                 <  1:711358/143‑1 (MQ=255)
                                                                                                                   gTTACCCGCACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGCCCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACTTAATTCAGTGCCTCATCACAGGAAACAGGTTAGACTTTGCCGCCCCCCTGGCGCCGGATCATgg                               >  2:230027/1‑149 (MQ=255)
                                                                                                                    ttACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCAGGACCGTACATCTCGTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACTCTGGCGACGGTTCATggg                              >  1:921552/1‑149 (MQ=255)
                                                                                                                        ccGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGtt                          >  2:69122/1‑149 (MQ=255)
                                                                                                                           tACCCCGCGACGAGCCGGTCGTCTGGCTGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTc                       <  2:1019584/149‑1 (MQ=255)
                                                                                                                            aCCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCa                      >  1:270895/1‑149 (MQ=255)
                                                                                                                              cACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAgg                    >  1:391618/1‑149 (MQ=255)
                                                                                                                              cACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAgg                    <  2:47504/149‑1 (MQ=255)
                                                                                                                                   gaggagCCGGTCGTCTGGATGGCATCCGGACCGGACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGCGCgg               <  2:138385/149‑1 (MQ=255)
                                                                                                                                   gaggagCCGGTCGGCTGGATGGCATCCGGACCGGAATTCTCTTTGATGGAGCTCAGGCGCTCGGCAACGTAAATCAGTGCATCTTCCCAGGAAACAGCGTCGAGTTAGCCGCCACGCTGGCGACGGATCATGGGGGGGTGCCGGCGCgg               <  2:212559/149‑1 (MQ=255)
                                                                                                                                       agCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGCTGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGCGCGGGGTc           <  1:534442/149‑1 (MQ=255)
                                                                                                                                         ccGGTCGTCTGGATGGCATCCGGACAGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCCACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGCGCGGGGTCAg         <  2:831209/149‑1 (MQ=255)
                                                                                                                                           ggTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCGAGGAAACAGGTTCGAGTTTGCCGCCCCGCTGGCGACGGATCATGGGGGTGTTCAGGCGCGGGGTCAGGa       >  1:253462/1‑149 (MQ=255)
                                                                                                                                            gtcgtcTGGATGGCATCCGGACCGGACTTCTCTTTGATGGCGCTCAGGCGCGCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGCGCGGGGTCAGGAt      <  2:166308/149‑1 (MQ=255)
                                                                                                                                              cgtcTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGCGCGGGGTCAGGATCt    <  2:387170/149‑1 (MQ=255)
                                                                                                                                               gtcCGGATGCCATCCAGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCa                      <  2:518813/130‑1 (MQ=255)
                                                                                                                                               gtcCGGATGCCATCCAGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCa                      >  1:518813/1‑130 (MQ=255)
                                                                                                                                                tcTGGTTGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCGGCCACGCTGGCGACGGATCATGGGGGGTTTCAGGCGCGGGGTCAGGATCTgg  <  2:610763/149‑1 (MQ=255)
                                                                                                                                                gcTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGCGCGGGGTCAGGATCTgg  <  2:8715/148‑1 (MQ=255)
                                                                                                                                                gcTGGATGGCAACCGGACCGTACTTCTCTTTGATGGCGCTCAGCCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGCGCGGGGTCAGGATCTgg  <  2:520034/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CGACCGATTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGCGCGGGGTCAGGATCTGG  >  NC_000913/4298910‑4299202

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: ‑