Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,141,799 G→C 5.3% F139L (TTC→TTG dcd ← deoxycytidine triphosphate deaminase; dCTP deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,141,7990GC 5.3% 413.2 / 12.0 113F139L (TTC→TTGdcddeoxycytidine triphosphate deaminase; dCTP deaminase
Reads supporting (aligned to +/- strand):  ref base G (54/53);  new base C (3/3);  total (57/56)
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.54e-02

Only 100 of 113 total aligned reads displayed.

CTGGCTACCGCGCCCTGCTGGTTGCGATATTTCGCATCTTCACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACC  >  NC_000913/2141653‑2141943
                                                                                                                                                  |                                                                                                                                                
ctggctACCGCGCCCTGCTGGTTGCGATATTTCGCATCTTCACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGaa                                                                                                                                                <  4:743108/149‑1 (MQ=255)
ctggctACCGCGCCCTGCTGGTTGCGATATTTCGCATCTTCACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGCACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGaa                                                                                                                                                <  4:1380170/149‑1 (MQ=255)
    ctACCGCGCCCTGCTGGTTGCGATATTTCGCATCTTCACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTcc                                                                                                                                            >  2:104865/1‑149 (MQ=255)
       ccGCGCCCTGCTGGTTGCGATATTTCGCATCTTCACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGc                                                                                                                                         >  1:1553266/1‑149 (MQ=255)
             ccTGCTGGTTGCGATATTTCGCATCTTCACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATg                                                                                                                                   >  1:1343757/1‑149 (MQ=255)
                 ctgGTTGCGATATTTCGCATCTTCACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAac                                                                                                                               <  2:1343757/149‑1 (MQ=255)
                 ctgGTTGCGATATTTCGCATCTTCACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAac                                                                                                                               <  2:248053/149‑1 (MQ=255)
                 ctgGTTGCGATATTTCGCATCTTCACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAac                                                                                                                               <  1:317759/149‑1 (MQ=255)
                        cGATATTTCGCATCTTCACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAaccagacca                                                                                                                        >  4:513570/1‑149 (MQ=255)
                         gATATTTCGCATCTTCACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAAt                                                                                                                                    <  2:386379/136‑1 (MQ=255)
                         gATATTTCGCATCTTCACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAaccagaccag                                                                                                                       <  3:388668/149‑1 (MQ=255)
                         gATATTTCGCATCTTCACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAaccagaccag                                                                                                                       >  2:990584/1‑149 (MQ=255)
                          atatTTCGCATCTTCACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCa                                                                                                                                           >  1:1204694/1‑128 (MQ=255)
                          atatTTCGCATCTTCACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCa                                                                                                                                           <  2:1204694/128‑1 (MQ=255)
                          atatTTCGCATCTTCACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGc                                                                                                                      >  3:388887/1‑149 (MQ=255)
                              ttCGCATCTTCACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCa                                                                                                                                           <  1:649286/124‑1 (MQ=255)
                              ttCGCATCTTCACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCa                                                                                                                                           >  2:649286/1‑124 (MQ=255)
                                    tCTTCACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGAt                                                                                                            >  3:717487/1‑149 (MQ=255)
                                      ttCACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGc                                                                                                          >  1:1601328/1‑149 (MQ=255)
                                       tCACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCg                                                                                                         <  2:1087203/149‑1 (MQ=255)
                                        cACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCgg                                                                                                        <  4:717487/149‑1 (MQ=255)
                                           ggcggTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTgc                                                                                                     >  1:1558560/1‑149 (MQ=255)
                                            gcggTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTgcg                                                                                                    >  2:329870/1‑149 (MQ=255)
                                                 tGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTg                                                                                               <  2:1558560/149‑1 (MQ=255)
                                                  gTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATg                                                                                                                                   >  3:1863357/1‑112 (MQ=255)
                                                  gTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATg                                                                                                                                   <  4:1863357/112‑1 (MQ=255)
                                                   tAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGa                                                                                                                 >  1:1445313/1‑129 (MQ=255)
                                                   tAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGa                                                                                                                 <  2:1445313/129‑1 (MQ=255)
                                                   tAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCgg                                                                                                        >  3:67884/1‑138 (MQ=255)
                                                   tAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCgg                                                                                                        <  4:67884/138‑1 (MQ=255)
                                                    aaGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACg                                                                                            <  4:388887/149‑1 (MQ=255)
                                                     aGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGt                                                                                           <  4:1695747/149‑1 (MQ=255)
                                                      ggTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCaa                                                                                                                                <  1:1259337/111‑1 (MQ=255)
                                                       gTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGc                                                                                         <  1:990584/149‑1 (MQ=255)
                                                        tCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGAt                                                                                                                <  1:1389240/125‑1 (MQ=255)
                                                        tCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGAt                                                                                                                >  2:1389240/1‑125 (MQ=255)
                                                        tCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCa                                                                                        <  3:513570/149‑1 (MQ=255)
                                                        tCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCa                                                                                        <  2:918304/149‑1 (MQ=255)
                                                           cACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATg                                                                                                           <  1:444637/127‑1 (MQ=255)
                                                           cACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATg                                                                                                           >  2:444637/1‑127 (MQ=255)
                                                              cgccgGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCa                                                                                  >  2:1542162/1‑149 (MQ=255)
                                                                      cggAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGa                                                                          >  4:629298/1‑149 (MQ=255)
                                                                      cggAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGa                                                                          <  1:329870/149‑1 (MQ=255)
                                                                      cggAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGa                                                                          >  2:261897/1‑149 (MQ=255)
                                                                       ggAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCg                                                                                                                   >  3:1892114/1‑107 (MQ=255)
                                                                       ggAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCg                                                                                                                   <  4:1892114/107‑1 (MQ=255)
                                                                       ggAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGAc                                                                         >  1:1667746/1‑149 (MQ=255)
                                                                       ggAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGAc                                                                         >  3:592756/1‑149 (MQ=255)
                                                                                cTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGAc                                                                         >  4:667206/1‑140 (MQ=255)
                                                                                cTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGAc                                                                         <  3:667206/140‑1 (MQ=255)
                                                                                cTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTg                                                                >  2:886323/1‑149 (MQ=255)
                                                                                 tCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGa                                                               >  2:570796/1‑149 (MQ=255)
                                                                                 tCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCACCCCCAGACGCGCCAGTGa                                                               >  4:1059115/1‑149 (MQ=255)
                                                                                   aaaGCTCAGCGTACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAaccagaccag                                                                                                                       >  3:1081659/1‑91 (MQ=255)
                                                                                   aaaGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAaccagaccag                                                                                                                       <  4:1081659/91‑1 (MQ=255)
                                                                                   aaaGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAgg                                                             >  2:518456/1‑149 (MQ=255)
                                                                                   aaaGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAgg                                                             <  2:1601328/149‑1 (MQ=255)
                                                                                     aGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGaa                                                           <  3:629298/149‑1 (MQ=255)
                                                                                          aGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAaccagacc                                                                                                                         >  2:1572609/1‑82 (MQ=255)
                                                                                          aGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGccc                                                      <  1:886323/149‑1 (MQ=255)
                                                                                            cgcACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGt                                                    <  1:570796/149‑1 (MQ=255)
                                                                                               aCCAATTAACATGCCCGGGCGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTcca                                                 >  2:1066681/1‑149 (MQ=255)
                                                                                               aCCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGaa                                                           >  3:1360102/1‑139 (MQ=255)
                                                                                               aCCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGaa                                                           <  4:1360102/139‑1 (MQ=255)
                                                                                                 cAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTccagc                                               <  2:1553266/149‑1 (MQ=255)
                                                                                                  aaTTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTccagcc                                              >  4:747795/1‑149 (MQ=255)
                                                                                                   aTTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAAACGCGCCAGTGAGGAACGCCCGTccagcca                                             <  2:1201711/149‑1 (MQ=255)
                                                                                                         aTGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCccac                                       >  2:373168/1‑149 (MQ=255)
                                                                                                                  cGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGc                              <  1:1542162/149‑1 (MQ=255)
                                                                                                                   gCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCa                                                                            >  2:730565/1‑102 (MQ=255)
                                                                                                                   gCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCa                                                                            <  1:730565/102‑1 (MQ=255)
                                                                                                                   gCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCt                             <  2:1667746/149‑1 (MQ=255)
                                                                                                                   gCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCt                             <  4:592756/149‑1 (MQ=255)
                                                                                                                   gCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCt                             <  1:261897/149‑1 (MQ=255)
                                                                                                                     aGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTgg                           >  4:429693/1‑149 (MQ=255)
                                                                                                                        gcCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCg                                                     <  3:566252/120‑1 (MQ=255)
                                                                                                                          cAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTccagccag                                            >  2:78740/1‑127 (MQ=255)
                                                                                                                          cAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTccagccag                                            <  3:1065711/127‑1 (MQ=255)
                                                                                                                          cAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTccagccag                                            >  4:1065711/1‑127 (MQ=255)
                                                                                                                          cAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTccagccag                                            <  1:78740/127‑1 (MQ=255)
                                                                                                                           aGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGt                     >  4:125683/1‑149 (MQ=255)
                                                                                                                                  gCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGt                                                                 <  3:1075977/98‑1 (MQ=255)
                                                                                                                                  gCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGt                                                                 >  4:1075977/1‑98 (MQ=255)
                                                                                                                                    ttACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCGATCTGGTGGGCTGGc                          >  3:1265727/1‑135 (MQ=255)
                                                                                                                                      aCCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCAc                  >  1:574935/1‑141 (MQ=255)
                                                                                                                                      aCCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCAc                  <  2:574935/141‑1 (MQ=255)
                                                                                                                                       ctggAGTTCTACAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTcca                                                 >  4:1754006/3‑109 (MQ=255)
                                                                                                                                       ctggAGTTCTACAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTcca                                                 <  3:1754006/107‑1 (MQ=255)
                                                                                                                                       ctggAGTTCTACAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTc                    <  3:1440648/136‑1 (MQ=255)
                                                                                                                                       ctggAGTTCTACAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAg         >  2:682578/3‑149 (MQ=255)
                                                                                                                                       ccGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGc                                                                                         >  3:926546/1‑69 (MQ=255)
                                                                                                                                        tggAGTTCTACAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGc        <  1:682578/148‑1 (MQ=255)
                                                                                                                                        cGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTccagccag                                            >  4:515324/1‑113 (MQ=255)
                                                                                                                                        cGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGc        <  3:747795/149‑1 (MQ=255)
                                                                                                                                        cGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGc        <  1:467770/149‑1 (MQ=255)
                                                                                                                                          gAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCGGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGt      >  3:1212851/1‑149 (MQ=255)
                                                                                                                                             ttGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCt                             >  1:1626544/1‑123 (MQ=255)
                                                                                                                                              tGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCAcc  >  1:1577899/1‑149 (MQ=255)
                                                                                                                                                tAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGa                                                               <  1:278181/86‑1 (MQ=255)
                                                                                                                                                tAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGa                                                               >  2:278181/1‑86 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
CTGGCTACCGCGCCCTGCTGGTTGCGATATTTCGCATCTTCACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACC  >  NC_000913/2141653‑2141943

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