Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,448,750 C→T 11.1% Q49* (CAG→TAG)  smrB → putative DNA endonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,448,7500CT11.1% 196.4 / 13.3 72Q49* (CAG→TAG) smrBputative DNA endonuclease
Reads supporting (aligned to +/- strand):  ref base C (29/35);  new base T (4/4);  total (33/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

AAATGAAAAAGAAAACAACACTTAGCGAGGAGGACCAGGCTCTGTTTCGCCAGTTGATGGCGGGGACTCGCAAGATTAAGCAGGACACGATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAAACTGCGCCGTGGCGATTATTCGCCGGAGTTGTTTTT  >  NC_000913/2448604‑2448898
                                                                                                                                                  |                                                                                                                                                    
aaaTGAAAAAGAAAACAACACTTAGCGAGGAGGACCAGGCTCTGTTTCGCCAGTTGATGGCGGGGACTCGCAAGATTAAGCAGGACACGATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCTAg                                                                                                                                                    >  7:464045/1‑149 (MQ=255)
 aaTGAAAAAGAAAACAACACTTAGCGAGGAGGACCAGGCTCTGTTTCGCCAGTTGATGGCGGGGACTCGCAAGATTAAGCAGGACACGATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAg                                                                                                                                                    <  7:316461/148‑1 (MQ=255)
 aaTGAAAAAGAAAACAACACTTAGCGAGGAGGACCAGGCTCTGTTTCGCCAGTTGATGGCGGGGACTCGCAAGATTAAGCAGGACACGATTGTCCACCGACCGCAGCCTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCTAg                                                                                                                                                    <  5:2576/148‑1 (MQ=255)
   tGAAAAAGAAAACAACACTTAGCGAGGAGGACCAGGCTCTGTTTCGCCAGTTGATGGCGGGGACTCGCAAGATTAAGCAGGACACGATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAg                                                                                                                                                 >  1:280510/1‑149 (MQ=255)
   tGAAAAAGAAAACAACACTTAGCGAGGAGGACCAGGCTCTGTTTCGCCAGTTGATGGCGGGGACTCGCAAGATTAAGCAGGACACGATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAg                                                                                                                                                 >  7:79273/1‑149 (MQ=255)
     aaaaaGAAAACAACACTTAGCGAGGAGGACCAGGCTCTGTTTCGCCAGTTGATGGCGGGGACTCGCAAGATTAAGCAGGACACGATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAg                                                                                                                                                 >  1:360437/1‑147 (MQ=255)
      aaaaGAAAACAACACTTAGCGAGGAGGACCAGGCTCTGTTTCGCCAGTTGATGGCGGGGACTCGCAAGATTAAGCAGGACACGATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAg                                                                                                                                              <  6:29490/149‑1 (MQ=255)
         aGAAAACAACACTTAGCGAGGAGGACCAGGCTCTGTTTCGCCAGTTGATGGCGGGGACTCGCAAGATTAAGCAGGACACGATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCt                                                                                                                                           <  1:281403/149‑1 (MQ=255)
          gAAAACAACACTTAGCGAGGAGGACCAGGCTCTGTTTCGCCAGTTGATGGCGGGGACTCGCAAGATTAAGCAGGACACGATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCt                                                                                                                                           >  3:271037/1‑148 (MQ=255)
           aaaaCAACACTTAGCGAGGAGGACCAGGCTCTGTTTCGCCAGTTGATGGCGGGGACTCGCAAGATTAAGCAGGACACGATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCTGa                                                                                                                                         >  7:168717/1‑149 (MQ=255)
            aaaCAACACTTAGCGAGGAGGACCAGGCTCTGTTTCGCCAGTTGATGGCGGGGACTCGCAAGATTAAGCAGGACACGATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCTGAt                                                                                                                                        <  8:168717/149‑1 (MQ=255)
                aacaCTTAGCGAGGAGGACCAGGCTCTGTTTCGCCAGTTGATGGCGGGGACTCGCAAGATTAAGCAGGACACGATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCTGATgc                                                                                                                                      >  3:129468/1‑147 (MQ=255)
                aacaCTTAGCGAGGAGGACCAGGCTCTGTTTCGCCAGTTGATGGCGGGGACTCGCAAGATTAAGCAGGACACGATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCTGATgc                                                                                                                                      >  4:284705/1‑147 (MQ=255)
                 acacTTAGCGAGGAGGACCAGGCTCTGTTTCGCCAGTTGATGGCGGGGACTCGCAAGATTAAGCAGGACACGATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCTGATgccag                                                                                                                                   <  5:85238/149‑1 (MQ=255)
                    cTTAGCGAGGAGGACCAGGCTCTGTTTCGCCAGTTGATGGCGGGGACTCGCAAGATTAAGCAGGACACGATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCTGATgccagcca                                                                                                                                <  3:284705/149‑1 (MQ=255)
                     ttAGCGAGGAGGACCAGGCTCTGTTTCGCCAGTTGATGGCGGGGACTCGCAAGATTAAGAAGGACACGATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCTAGGAGCAGGCTGATGCCAGCCat                                                                                                                               >  3:136006/1‑149 (MQ=255)
                              aggaCCAGGCTCTGTTTCGCCAGTTGATGGCGGGGACTCGCAAGATTAAGCAGGACACGATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCTAGGAGCAGGCTGATGCCAGCCATTATTTCTcc                                                                                                                      <  5:349332/149‑1 (MQ=255)
                                   cAGGCTCTGTTTCGCCAGTTGATGGCGGGGACTCGCAAGATTAAGCAGGACACGATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAg                                                                                                                                              <  4:166124/120‑1 (MQ=255)
                                   cAGGCTCTGTTTCGCCAGTTGATGGCGGGGACTCGCAAGATTAAGCAGGACACGATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAg                                                                                                                                              >  3:166124/1‑120 (MQ=255)
                                   cAGGCTCTGTTTCGCCAGTTGATGGCGGGGACTCGCAAGATTAAGCAGGACACGATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGa                                                                                                                 <  4:205490/149‑1 (MQ=255)
                                        tctGTTTCGCCAGTTGATGGCGGGGACTCGCAAGATTAAGCAGGACACGATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGTTGAAACGCTTGATCCAGGa                                                                                                                                                  <  8:383039/111‑1 (MQ=255)
                                        tctGTTTCGCCAGTTGATGGCGGGGACTCGCAAGATTAAGCAGGACACGATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGttt                                                                                                             >  8:95837/1‑148 (MQ=255)
                                           gTTTCGCCAGTTGATGGCGGGGACTCGCAAGATTAAGCAGGACACGATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGc                                                                                                         >  2:285356/1‑149 (MQ=255)
                                                gCCAGTTGATGGCGGGGACTCGCAAGATTAAGCAGGACACGATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGtt                                                                                                     <  4:207628/148‑1 (MQ=255)
                                                 ccAGTTGATGGCGGGGACTCGCAAGATTAAGCAGGACACGATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCTAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGttat                                                                                                   <  8:464045/149‑1 (MQ=255)
                                                 ccAGTTGATGGCGGGGACTCGCAAGATTAAGCAGGACACGATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGttat                                                                                                   >  2:143935/1‑149 (MQ=255)
                                                  cAGTTGATGGCGGGGACTCGCAAGATTAAGCAGGACACGATAGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGttat                                                                                                   >  6:369137/1‑148 (MQ=255)
                                                     gtGATGGCGGGGACTCGCAAGATTAAGCAGGACACGATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTaaa                                                                                               >  7:330961/2‑149 (MQ=255)
                                                       gATGGCGGGGACTCGCAAGATTAAGCAGGACACGATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAAt                                                                                              <  1:285356/148‑1 (MQ=255)
                                                        aTGGCGGGGACTCGCAAGATTAAGCAGGACCCGATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATAc                                                                                            <  3:152351/149‑1 (MQ=255)
                                                        aTGGCGGGGACTCGCAAGATTAAGCAGGACACGATTGTCCACCGACCGCAGCGTAAAAAACTCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATa                                                                                             <  6:49791/148‑1 (MQ=255)
                                                        aTGGCGGGGACTCGCAAGATTAAGCAGGACACGATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATa                                                                                             <  7:346439/148‑1 (MQ=255)
                                                            cGGGGACTCGCAAGATTAAGCAGGACACGATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGaa                                                                                        <  8:79273/149‑1 (MQ=255)
                                                             ggggACTCGCAAGATTAAGCAGGACACGATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAg                                                                                       <  6:394503/149‑1 (MQ=255)
                                                               ggACTCGCAAGATTAAGCAGGACACGATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCTGATGCCATCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGt                                                                                     >  8:215912/1‑149 (MQ=255)
                                                               ggACTCGCAAGATTAAGCAGGACACGATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGt                                                                                     >  4:403504/1‑149 (MQ=255)
                                                                   tCGCAAGATTAAGCAGGACACGATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCTAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCgg                                                                                 >  2:238937/1‑149 (MQ=255)
                                                                   tCGCAAGATTAAGCAGGACACGATTGTACACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGCAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCgg                                                                                 >  8:295148/1‑149 (MQ=255)
                                                                     gCAAGATTAAGCAGGACACGATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTTAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGt                                                                                <  6:166496/148‑1 (MQ=255)
                                                                     gCAAGATTAAGCAGGACACGATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGGTATTAAATACCGAAGGTCCGGt                                                                                <  2:280510/148‑1 (MQ=255)
                                                                      cAAGATTAAGCAGGACACGATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTg                                                                               >  4:273465/1‑148 (MQ=255)
                                                                        aGATTAAGCAGGACACGATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGaaa                                                                            <  4:129468/149‑1 (MQ=255)
                                                                         gATTAAGCAGGACACGATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTc                                                                                    >  5:244199/1‑140 (MQ=255)
                                                                         gATTAAGCAGGACACGATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTc                                                                                    <  6:244199/140‑1 (MQ=255)
                                                                               gCAGGACACGATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCTAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTc                                                                     <  4:136006/149‑1 (MQ=255)
                                                                               gCAGGACACGATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTc                                                                     >  6:393836/1‑149 (MQ=255)
                                                                                     cacGATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCgg                                                               <  1:413315/149‑1 (MQ=255)
                                                                                      acGATTGTCCACCGACCGCAGCGTAAAACAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCTGATGCCCGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGa                                                              >  5:102333/1‑149 (MQ=255)
                                                                                      acGATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTc                                                                     <  2:345400/142‑1 (MQ=255)
                                                                                      acGATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTc                                                                     >  1:345400/1‑142 (MQ=255)
                                                                                        gATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGttatt                                                                                                  <  8:210835/111‑1 (MQ=255)
                                                                                        gATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGttatt                                                                                                  >  7:210835/1‑111 (MQ=255)
                                                                                            gTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGCTCCAGGAGCAGGCTGCTGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATCCCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCa                                                         <  5:337547/148‑1 (MQ=255)
                                                                                            gTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAg                                                        >  7:413977/1‑149 (MQ=255)
                                                                                                 ccgaccgCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCAtt                                                   <  3:384292/149‑1 (MQ=255)
                                                                                                    accgCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAATGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCAtttt                                                 >  6:401744/1‑148 (MQ=255)
                                                                                                     ccgCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGa                                               <  7:295148/149‑1 (MQ=255)
                                                                                                      cgCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAg                                              <  5:401744/149‑1 (MQ=255)
                                                                                                      cgCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAg                                              <  5:393883/149‑1 (MQ=255)
                                                                                                      cgCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAg                                              <  2:360437/149‑1 (MQ=255)
                                                                                                      cgCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAg                                              <  4:271037/149‑1 (MQ=255)
                                                                                                               aaaaaaTCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGaaa                                     >  3:410284/1‑149 (MQ=255)
                                                                                                                aaaaaTCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGa                                               >  7:416535/1‑138 (MQ=255)
                                                                                                                aaaaaTCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGa                                               <  8:416535/138‑1 (MQ=255)
                                                                                                                             gTGCCGGTGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAAACTGCGCCGTGGCGa                       <  8:330961/149‑1 (MQ=255)
                                                                                                                                ccGGTGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAAACTGCGCCGTGGCGatta                    <  5:393836/149‑1 (MQ=255)
                                                                                                                                   gTGAAACGCTTGATCTAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAAACTGCGCCGTGGCGATTATTc                 >  7:448706/1‑149 (MQ=255)
                                                                                                                                      aaaCGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAAACTGCGCCGTGGCGATTATTCGcc              >  2:442052/1‑149 (MQ=255)
                                                                                                                                      aaaCGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGGCAGCCATTTTGAGGCGAAGAAACTGCGCCGTGGCGATTATTCGc               <  5:369137/148‑1 (MQ=255)
                                                                                                                                           cTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATCTGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAAACTGCGCCGTGGCGATTATTCGCCGGAgt         <  1:143935/149‑1 (MQ=255)
                                                                                                                                                tCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAAACTGCGCCGTGGCGATTATTCGCCGGAGTTGttt    >  1:172039/1‑149 (MQ=255)
                                                                                                                                                tCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAAACTGCGCCGTGGCGATTATTCGCCGGAGTTGttt    >  1:165517/1‑149 (MQ=255)
                                                                                                                                                  cAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAAACTGCGCCGTGGCGATTATTCGCCGGAGTTGttttt  >  4:126540/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
AAATGAAAAAGAAAACAACACTTAGCGAGGAGGACCAGGCTCTGTTTCGCCAGTTGATGGCGGGGACTCGCAAGATTAAGCAGGACACGATTGTCCACCGACCGCAGCGTAAAAAAATCAGCGAAGTGCCGGTGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAAACTGCGCCGTGGCGATTATTCGCCGGAGTTGTTTTT  >  NC_000913/2448604‑2448898

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