Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,146,922 A→G 5.0% pseudogene (351/417 nt) yghX ← pseudogene, dienlactone hydrolase family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,146,9220AG 5.0% 250.3 / 2.7 79pseudogene (351/417 nt)yghXpseudogene, dienlactone hydrolase family
Reads supporting (aligned to +/- strand):  ref base A (33/42);  new base G (2/2);  total (35/44)
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

GATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAACCCCCGGATAGATATACGCCTCATAAACCTTATTATTGGCTTTCAACGCCGCCTCGTAAGCAGGCCAGCCCTCGTTGATTCGGGTGTCCAGTTCCGCGAAGTGGAGTAGT  >  NC_000913/3146774‑3147069
                                                                                                                                                    |                                                                                                                                                   
gATATTTTAGGCAAATAACGAAATTTGCGCACGCTATTTTTTCATGGTTGTTTAAGACGTAATCCCGTAATAGATAAACCTATCAGGAGAGATATATATCGAACCATTTAAGTGTCATTTGCCAGGCAAGATACGCGCCAGATTTGTCa                                                                                                                                                     >  4:3698/1‑149 (MQ=255)
gATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCa                                                                                                                                                     <  3:13126/149‑1 (MQ=255)
 atatTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCAt                                                                                                                                                    <  3:199404/149‑1 (MQ=255)
    ttttAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATaa                                                                                                                                                  <  5:2378/148‑1 (MQ=255)
            aaaGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTg                                                                                                                                         >  6:10245/1‑149 (MQ=255)
            aaaGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTg                                                                                                                                         >  8:157309/1‑149 (MQ=255)
            aaaGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTg                                                                                                                                         <  5:10245/149‑1 (MQ=255)
            aaaGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTg                                                                                                                                         <  8:229664/149‑1 (MQ=255)
            aaaGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTg                                                                                                                                         <  8:40029/149‑1 (MQ=255)
              aGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGa                                                                                                                                       <  1:147775/149‑1 (MQ=255)
                aaCGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAAt                                                                                                                                     >  2:145701/1‑149 (MQ=255)
                 aCGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCGTAACGGGGCGTGGAAt                                                                                                                                     >  4:141866/1‑148 (MQ=255)
                 aCGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCGTAACGGGGCGTGGAAt                                                                                                                                     >  4:138823/1‑148 (MQ=255)
                       cttGCGCACGCTATCTTACCATGCCCGCATACGACGTACTCCCGTAAGAGATAAACCTATCAGGAGAGATACTTATCGAACCATTTCAGTGGCCTTTGCCAGGCAAGATCCGCGGCAGATTTGGCATAACGGGGCGTGGAAGCATTATg                                                                                                                              <  4:204527/148‑1 (MQ=255)
                        ttGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTACTCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGAGCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTAt                                                                                                                               <  3:190426/147‑1 (MQ=255)
                        ttGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATg                                                                                                                              <  2:8979/148‑1 (MQ=255)
                        ttGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATg                                                                                                                              <  7:157309/148‑1 (MQ=255)
                           cgcACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGaa                                                                                                                           <  3:232225/148‑1 (MQ=255)
                              aCGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCGTAACGGGGCGTGGAATCATTATGGAATCCg                                                                                                                       <  3:138823/149‑1 (MQ=255)
                              aCGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCGTAACGGGGCGTGGAATCATTATGGAATCCg                                                                                                                       <  3:141866/149‑1 (MQ=255)
                              aCGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCg                                                                                                                       >  1:11141/1‑149 (MQ=255)
                                gCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTACTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAAAGGGGCGTGGAATCATTATGGAATCCGt                                                                                                                      >  4:11024/1‑148 (MQ=255)
                                 cTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGa                                                                                                                    <  7:212736/149‑1 (MQ=255)
                                  gaTTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTGGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGAt                                                                                                                   >  6:163468/2‑149 (MQ=255)
                                      ttttCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAAc                                                                                                               >  7:240539/1‑149 (MQ=255)
                                       tttCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAAc                                                                                                               <  8:141855/148‑1 (MQ=255)
                                           aTGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACgg                                                                                                                                               >  3:124281/1‑112 (MQ=255)
                                           aTGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACgg                                                                                                                                               <  4:124281/112‑1 (MQ=255)
                                           aTGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAAccccc                                                                                                           <  6:14812/148‑1 (MQ=255)
                                              cTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAACCCCCGga                                                                                                        <  3:101800/148‑1 (MQ=255)
                                                tGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGt                                                                                                                      >  6:134935/1‑132 (MQ=255)
                                                tGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGt                                                                                                                      <  5:134935/132‑1 (MQ=255)
                                                 gTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAACCCCCGgatag                                                                                                     <  2:11141/148‑1 (MQ=255)
                                                  tttAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAACCCCCGgatagat                                                                                                   <  8:240539/149‑1 (MQ=255)
                                                  tttAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAACCCCCGgataga                                                                                                    >  8:209515/1‑148 (MQ=255)
                                                    tAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAAcccc                                                                                                            <  6:117312/138‑1 (MQ=255)
                                                    tAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAAcccc                                                                                                            >  5:117312/1‑138 (MQ=255)
                                                     aaGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAACCCCCGGATAGAtata                                                                                                >  8:152170/1‑149 (MQ=255)
                                                      aGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAACCCCCGGATAGAtata                                                                                                <  3:3698/148‑1 (MQ=255)
                                                      aGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAACCCCCGGATAGATATAc                                                                                               <  5:106764/149‑1 (MQ=255)
                                                       gccGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAACCCCCGGATAGATATACg                                                                                              >  2:161126/3‑149 (MQ=255)
                                                       gACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAACCCCCGGATAGATATACg                                                                                              >  8:154402/1‑149 (MQ=255)
                                                       gACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAACCCCCGGATAGATATACg                                                                                              >  6:73129/1‑149 (MQ=255)
                                                          gTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAACCCCCGGATAGATATACGCCt                                                                                           >  2:201105/1‑149 (MQ=255)
                                                          gTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAACCCCCGGATAGATATACGCCt                                                                                           >  7:172877/1‑149 (MQ=255)
                                                            aaTCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAACCCCCGGATAGATATAc                                                                                               >  3:69132/1‑143 (MQ=255)
                                                            aaTCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAACCCCCGGATAGATATAc                                                                                               <  4:69132/143‑1 (MQ=255)
                                                              tCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAACCCCCGGATAGATATACGCCTCAt                                                                                        >  8:165192/1‑148 (MQ=255)
                                                                 cGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAACCCCCGGATAGATATACGCCTCATAGAc                                                                                    >  3:26964/1‑149 (MQ=255)
                                                                 cGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAACCCCCGGATAGATATACGCCTCATAAAc                                                                                    <  4:200215/149‑1 (MQ=255)
                                                                  gTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAACCCCCGGATAGATATACGCCTCa                                                                                         <  4:188284/143‑1 (MQ=255)
                                                                  gTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAACCCCCGGATAGATATACGCCTCa                                                                                         >  3:188284/1‑143 (MQ=255)
                                                                     agagATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAACCCCCGGATAGATATACGCCTCATAAACCtta                                                                                <  6:98676/149‑1 (MQ=255)
                                                                      gagaTAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAACCCCCGGATAGATATACGCCTCATAAACCttat                                                                               >  3:220728/1‑149 (MQ=255)
                                                                            aaCCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAACCCCCGGATAGATATACGCCTCATAAACCTTATTATTg                                                                          <  7:154402/148‑1 (MQ=255)
                                                                             aCCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAACCCCCGGATAGATATACGCCTCATAAACCTTATTATTGGc                                                                        <  1:145701/149‑1 (MQ=255)
                                                                             aCCTATCAGGAGAGATATTTATCGAACCATTTCAGGTTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACTTGGCGTGGAATCATTATGGAATCCGTGATTAACCCCCGTATAGATATACGCCTCATAAACCTTATTATTGGc                                                                        <  1:155558/149‑1 (MQ=255)
                                                                               cTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTAGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAACCCCCGGATAGATATACGCCTCATAAACCTTATTATTGGCtt                                                                      >  3:128667/1‑149 (MQ=255)
                                                                                  tCAGGAGAGATATTTATCGAAGCATTTCAGTGTCCTTTGCCTGGCAAGATCGGCGGCAGATTTGTCATAACGGGTCGTGGAATCATTATGGAATCCGTGATTAACCCCCGGATAGATATACGCCTCATAAACCTTATTATTGGCTTTc                                                                    <  5:105728/148‑1 (MQ=255)
                                                                                  tCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAACCCCCGGATAGATATACGCCTCATAAACCTTATTATTGGCttt                                                                     <  5:96619/147‑1 (MQ=255)
                                                                                              tttATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAACCCCCGGATAGATATACGCCTCATAAACCTTATTATTGGCTTTCAACGCCGCCTCg                                                        >  3:226263/1‑148 (MQ=255)
                                                                                               ttATCGAACCATTTCAGTGTCCTTTGCCAGCCAAGATCGGCGTCAGATTTGTCATAACGGGGCGTGGAAACATTATCGAATCCGTGATTAACCCCCGGATAGATTTACGCCTCATAAACGTTTTTATTGGCTTTCAACGCCGCCTCGt                                                       >  6:115913/1‑148 (MQ=255)
                                                                                                 aTCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAACCCCCGGATAGATATACGCCTCATAAACCTTATTATTGGCTTTCAACGCCGCCTCGTaa                                                     <  1:201105/148‑1 (MQ=255)
                                                                                                 aTCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATAGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAACCCCCGGATAGATATACGCCTCATAAACCTTATTATTGGCTTTCAACGCCGCCTCGTaa                                                     <  4:75185/148‑1 (MQ=255)
                                                                                                    gAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAACCCCCGGATAGATATACGCCTCATAAACCTTATTATTGGCTTTCAACGCCGCCTCGTAAGCAg                                                 <  7:209515/149‑1 (MQ=255)
                                                                                                    gAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGGGGAATCATTATGGAATCCGTGATTAACCCCCGGATAGATATACGCCTCATAAACCTTATTATTGGCTTTCAACGCCGCCTCGTAAGCAg                                                 <  5:73129/149‑1 (MQ=255)
                                                                                                     aaCCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAACCCCCGGATAGATATACGCCTCATAAACCTTATTATTGGCTTTCAACGCCGCCTCGTAAGCAgg                                                >  3:231518/1‑149 (MQ=255)
                                                                                                      aCCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAACCCCCGGATAGATATACGCCTCATAAACCTTATTATTGGCTTTCAACGCCGCCTCGTAAGCAGGc                                               <  1:96227/149‑1 (MQ=255)
                                                                                                         aTTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGt                                                                                                                      >  1:208494/1‑75 (MQ=255)
                                                                                                         aTTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGt                                                                                                                      <  2:208494/75‑1 (MQ=255)
                                                                                                            tCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAACCCCCGGATAGATATACGCCTCATAAACCTTATTATTGGCTTTCAACGCCGCCTCGTAAGCAGGCCAGccc                                         >  3:4661/1‑149 (MQ=255)
                                                                                                                       tGCCAGGCAAGATCGGCGGCAGATTTGGCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAACCCCCGGATAGATATACGCCTCATAAACCTTATTATTGGCTTTCAACGCCGCCTCGTAAGCAGGCCAGCCCTCGTTGATTc                               <  8:90412/148‑1 (MQ=255)
                                                                                                                            ggCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAACCCCCGGATAGATATACGCCTCATAAACCTTATTATTGGCTTTCAACGCCGCCTCGTAAGCAGGCCAGCCCTCGTTGATTCGGgtgt                         >  3:116484/1‑149 (MQ=255)
                                                                                                                              cAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAACCCCCGGATAGATATACGCCTCATAAACCTTATTATTGGCTTTCAACGCCGCCTCGTAAGCAGGCCAGCCCTCGTTGATTCGGGTGTc                        >  3:105138/1‑148 (MQ=255)
                                                                                                                                 gATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAACCCCCGGATAGATATACGCCTCATAAACCTTATTATTGGCTTTCAACGCCGCCTCGTAAGCAGGCCAGCCCTCGTTGATTCGGGTGTCCAGt                    >  8:66524/1‑149 (MQ=255)
                                                                                                                                           aGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAACCCCCGGATAGATATACGCCTCATAAACCTTATTATTGGCTTTCAACGCCGCCTCGTAAGCAGGCCAGCCCTCGTTGATTCGGGTGTCCAGTTCCGCGAAGt          >  4:100412/1‑149 (MQ=255)
                                                                                                                                             aTTTGTCATAACGGGCGGGGGAATCATTATGGAATCCGTGCTAAACCCCCGGATAGATATACGCCTCATAAACCTTATTATTGGCTTTCAACGCCGCCTCGTAAGCAGGCCAGCCCTCGTAGATTCGGGCGTCCAGTTCCGCGAAGTgg        <  4:128667/149‑1 (MQ=255)
                                                                                                                                              tttGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAACCCCCGGATAGATATACGCCTCATAAACCTTATTATTGGCTTTCAACGCCGCCTCGTAAGCAGGCCAGCCCTCGTTGATTCGGGTGTCCAGTTCCGCGAAGTGGa       <  8:172877/149‑1 (MQ=255)
                                                                                                                                              tttGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAACCCCCGGATAGATATACGCCTCATAAACCTTATTATTGGCTTTCAACGCCGCCTCGTAAGCAGGCCAGCCCTCGTTGATTCGGGTGTCCAGTTCCGCGAAGTGGa       <  4:220728/149‑1 (MQ=255)
                                                                                                                                                tGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAACCCCCGGATAGATATACGCCTCATAAACCTTATTATTGGCTTTCAACGCCGCCTCGTAAGCAGGCCAGCCCTCGTTGATTCGGGTGTCCAGTTCCGCGAAGTGGag      <  4:105138/148‑1 (MQ=255)
                                                                                                                                                  tCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAACCCCCGGATAGATATACGCCTCATAAACCTTATTATTGGCTTTCAACGCCGCCTCGTAAGCAGGCCAGCCCTCGTTGATTCGGGTGTCCAGTTCCGCGAAGTGGagtag   >  8:221429/1‑149 (MQ=255)
                                                                                                                                                    aTAACGGGGCGTGGAATCATTATGGAATCCGTGATTAACCCCCGGATAGATATACGCCTCATAAACCTTATTATTGGCTTTCAACGCCGCCTCGTAAGCAGGCCAGCCCTCGTTGATTCGGGTGTCCAGTTCCGCGAAGTGGagtagt  >  2:31847/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGGCAGATTTGTCATAACGGGGCGTGGAATCATTATGGAATCCGTGATTAACCCCCGGATAGATATACGCCTCATAAACCTTATTATTGGCTTTCAACGCCGCCTCGTAAGCAGGCCAGCCCTCGTTGATTCGGGTGTCCAGTTCCGCGAAGTGGAGTAGT  >  NC_000913/3146774‑3147069

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