Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,007,324 C→T 5.7% A487A (GCC→GCT recQ → ATP‑dependent DNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,007,3240CT 5.7% 227.8 / 3.4 70A487A (GCC→GCTrecQATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand):  ref base C (34/32);  new base T (2/2);  total (36/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 = 3.69e-01

ATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCG  >  NC_000913/4007178‑4007472
                                                                                                                                                  |                                                                                                                                                    
aTTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCa                                                                                                                                                    >  3:452666/1‑149 (MQ=255)
aTTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCa                                                                                                                                                    <  2:489146/149‑1 (MQ=255)
          cTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGcc                                                                                                                                          <  2:271575/149‑1 (MQ=255)
          cTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGcc                                                                                                                                          >  3:107084/1‑149 (MQ=255)
           taataaCCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGccc                                                                                                                                         >  2:251280/1‑149 (MQ=255)
             ataaCCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTa                                                                                                                                       <  3:296256/149‑1 (MQ=255)
              taaCCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTAc                                                                                                                                      >  1:608232/1‑149 (MQ=255)
                 ccAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAAc                                                                                                                                   >  1:276188/1‑149 (MQ=255)
                 ccAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAAc                                                                                                                                   >  1:450473/1‑149 (MQ=255)
                    gCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGa                                                                                                                                <  2:362597/149‑1 (MQ=255)
                    gCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGa                                                                                                                                >  2:376972/1‑149 (MQ=255)
                     cGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGAc                                                                                                                               <  2:582600/149‑1 (MQ=255)
                        aTCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACaga                                                                                                                            <  2:608232/149‑1 (MQ=255)
                            gcgACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCTCAGCATTCTGCCCTACAACTGACAGAGGcc                                                                                                                        >  1:402865/1‑149 (MQ=255)
                            gcgACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGGGACGCAAAATATTGCTCAGCATTCTGCCCTACAACTGACAGAGGcc                                                                                                                        >  1:226352/1‑149 (MQ=255)
                                cTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCcgcgc                                                                                                                    <  4:452666/149‑1 (MQ=255)
                                cTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCcgcgc                                                                                                                    <  4:74061/149‑1 (MQ=255)
                                      tCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCgg                                                                                                              >  3:617902/1‑149 (MQ=255)
                                      tCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCgg                                                                                                              >  3:77514/1‑149 (MQ=255)
                                        aTGACAAACTGAAAGTCTATGGCGTGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGtg                                                                                                            >  4:47900/1‑149 (MQ=255)
                                         tGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGtgc                                                                                                           <  4:169485/149‑1 (MQ=255)
                                         tGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGtgc                                                                                                           >  4:431279/1‑149 (MQ=255)
                                            cAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGtgctgc                                                                                                        <  1:460721/149‑1 (MQ=255)
                                                cTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCgcgg                                                                                                    <  1:251280/149‑1 (MQ=255)
                                                  gAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCgcggcg                                                                                                  >  1:4476/1‑149 (MQ=255)
                                                  gAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCgcggcg                                                                                                  >  4:751791/1‑149 (MQ=255)
                                                   aaaGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCTCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGa                                                                                                 <  2:402865/149‑1 (MQ=255)
                                                   aaaGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCTCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGa                                                                                                 <  2:226352/149‑1 (MQ=255)
                                                     aGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAAt                                                                                               <  3:289175/149‑1 (MQ=255)
                                                      gTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGaa                                                                                                >  3:566833/1‑147 (MQ=255)
                                                      gTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGaa                                                                                                <  4:566833/147‑1 (MQ=255)
                                                      gTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATc                                                                                              <  4:77514/149‑1 (MQ=255)
                                                      gTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATc                                                                                              <  1:763539/149‑1 (MQ=255)
                                                         tATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCctc                                                                                           >  1:433881/1‑149 (MQ=255)
                                                            ggCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCtt                                                                                         <  2:220018/148‑1 (MQ=255)
                                                            ggCATCGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCtt                                                                                         >  1:220018/1‑148 (MQ=255)
                                                             gCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCctc                                                                                           <  1:181803/145‑1 (MQ=255)
                                                             gCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCctc                                                                                           >  2:181803/1‑145 (MQ=255)
                                                             gCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTg                                                                                       >  1:394176/1‑149 (MQ=255)
                                                                   gtcGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACtt                                                                                 >  3:588641/3‑149 (MQ=255)
                                                                       tGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgccg                                                                             >  2:67525/1‑149 (MQ=255)
                                                                          tAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgccgtgc                                                                          >  4:216381/1‑149 (MQ=255)
                                                                          tAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAACATTCTGCCCTACAACTGACAGAGGCCGCTCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgccgtgc                                                                          >  2:44224/1‑149 (MQ=255)
                                                                                  aTGAACATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATc                                                                  <  3:216381/149‑1 (MQ=255)
                                                                                     aaCATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTg                                                               <  4:107084/149‑1 (MQ=255)
                                                                                     aaCATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTg                                                               >  1:54851/1‑149 (MQ=255)
                                                                                     aaCATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTg                                                               >  4:488424/1‑149 (MQ=255)
                                                                                        aTTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGgcg                                                            >  2:123703/1‑149 (MQ=255)
                                                                                         ttGGGTGAGCGTGATCCGCCAGCTGATTCACCGCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGgcgc                                                           <  4:414873/149‑1 (MQ=255)
                                                                                                   gTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGaa                                                 >  4:94576/1‑149 (MQ=255)
                                                                                                    tGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGaaa                                                <  2:564957/149‑1 (MQ=255)
                                                                                                    tGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGaaa                                                <  2:37139/149‑1 (MQ=255)
                                                                                                    tGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGaaa                                                <  2:127527/149‑1 (MQ=255)
                                                                                                      aTCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGc                                              >  4:160068/1‑149 (MQ=255)
                                                                                                         cGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGAt                                           <  4:377077/149‑1 (MQ=255)
                                                                                                          gCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATg                                          <  1:474589/149‑1 (MQ=255)
                                                                                                             aGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAg                                       <  4:659252/149‑1 (MQ=255)
                                                                                                                      aCCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTc                              <  3:464001/149‑1 (MQ=255)
                                                                                                                      aCCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTc                              >  1:37981/1‑149 (MQ=255)
                                                                                                                      aCCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTc                              >  1:609949/1‑149 (MQ=255)
                                                                                                                       ccTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCaa                                                                                    <  2:735595/94‑1 (MQ=255)
                                                                                                                       ccTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCaa                                                                                    >  1:735595/1‑94 (MQ=255)
                                                                                                                       ccTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcg                             >  3:576175/1‑149 (MQ=255)
                                                                                                                        cTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcgg                            >  4:578018/1‑149 (MQ=255)
                                                                                                                        cTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcgg                            >  3:288052/1‑149 (MQ=255)
                                                                                                                         tCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcggc                           <  1:234017/149‑1 (MQ=255)
                                                                                                                              cTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCaa                      <  2:394176/149‑1 (MQ=255)
                                                                                                                                 gTGACGCAAAATATTGCCCAGCATTCTGCCCTCCAACTGAAAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCACTTTGCAACTTGACGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAGTCGTTCCGAGGCAACTa                   >  3:8901/1‑149 (MQ=255)
                                                                                                                                      gCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATc              <  1:703627/149‑1 (MQ=255)
                                                                                                                                           atatTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCaaa         <  1:123703/149‑1 (MQ=255)
                                                                                                                                                 cccAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTc   >  4:583981/1‑149 (MQ=255)
                                                                                                                                                 cccAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTc   <  3:484595/149‑1 (MQ=255)
                                                                                                                                                  ccAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCg  >  2:387364/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
ATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCG  >  NC_000913/4007178‑4007472

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