Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,398,548 G→A 6.2% S520S (TCG→TCA hsdS → type I restriction modification system specificity protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,398,5480GA 6.2% 183.8 / 2.1 65S520S (TCG→TCAhsdStype I restriction modification system specificity protein
Reads supporting (aligned to +/- strand):  ref base G (29/31);  new base A (2/2);  total (31/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.99e-01

TAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAACACCGCATTGTCACAAAGGTCGATCAGCTTATGAAACTATGCGATACGCTAAAGCAACAGATAAATGTAGCTCGCAGCAAGCAAACAGAGCT  >  NC_002947/5398400‑5398691
                                                                                                                                                     |                                                                                                                                                
tAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCg                                                                                                                                                  >  7:134488/1‑149 (MQ=255)
 aaTGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCg                                                                                                                                                  >  3:195519/1‑148 (MQ=255)
 aaTGATTTACGACAAGATCATCAGTTTTCATCCTGTTGATCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGTTCAATATCTCg                                                                                                                                                  <  8:88230/148‑1 (MQ=255)
  aTGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGATCAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCa                                                                                                                                                <  3:128894/149‑1 (MQ=255)
   tGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCaa                                                                                                                                               <  4:51651/149‑1 (MQ=255)
        tCCGCCACGATCATCAGCTCTCCTCCGGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTTCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGTCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGaaa                                                                                                                                           <  4:274293/148‑1 (MQ=255)
          cGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGc                                                                                                                                        <  6:326118/149‑1 (MQ=255)
           gACAAGATCATCAGGTTTCATCATGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCt                                                                                                                                       >  8:2130/1‑149 (MQ=255)
                     tcaGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCt                                                                                                                             >  5:299905/1‑149 (MQ=255)
                           ttCATCCTGTTGAGCTGGGTGTTTATGGTAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATc                                                                                                                       >  1:195070/1‑149 (MQ=255)
                                 cTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGtt                                                                                                                  <  2:195070/148‑1 (MQ=255)
                                     tGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGct                                                                                                             <  6:299905/148‑2 (MQ=255)
                                              tgtTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGc                                                                                                                              >  7:273129/1‑123 (MQ=255)
                                              tgtTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGc                                                                                                                              <  8:273129/123‑1 (MQ=255)
                                                tttATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGa                                                                                                  <  4:244032/148‑1 (MQ=255)
                                                   aTGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGtt                                                                                                                  >  4:304907/1‑130 (MQ=255)
                                                   aTGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGtt                                                                                                                  <  3:304907/130‑1 (MQ=255)
                                                   aTGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCaa                                                                                              >  3:171766/1‑149 (MQ=255)
                                                           tttGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCACAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAACACCg                                                                                         <  4:300060/146‑1 (MQ=255)
                                                           tttGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCACAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAACACCg                                                                                         >  3:300060/1‑146 (MQ=255)
                                                                   tctTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTgg                                                                                                                               <  2:184903/101‑1 (MQ=255)
                                                                   tctTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTgg                                                                                                                               >  1:184903/1‑101 (MQ=255)
                                                                     tttGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAACACCGCATTGTCACAAAg                                                                            >  6:246002/1‑149 (MQ=255)
                                                                        gCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCa                                                                                                             >  4:256812/1‑113 (MQ=255)
                                                                        gCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCa                                                                                                             <  3:256812/113‑1 (MQ=255)
                                                                        gCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAACACCGCATTGTCACAAAGGTc                                                                         >  4:183091/1‑149 (MQ=255)
                                                                        gCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAACACAGCATTGTCACAAAgg                                                                           >  4:195708/1‑147 (MQ=255)
                                                                          ttGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTcc                                                                                                      <  2:153491/118‑1 (MQ=255)
                                                                          ttGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTcc                                                                                                      >  1:153491/1‑118 (MQ=255)
                                                                          ttGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAACACCGCATTGTCACAAAGGTCGa                                                                       >  1:315943/1‑149 (MQ=255)
                                                                          ttGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCAGTTAGCACCGCTCCGCGAGCAACACCGCATTGTCACAAAGGTCGa                                                                       >  7:307890/1‑149 (MQ=255)
                                                                          ttGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCACAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAACACCGCATTGTCACAAAGGTCGa                                                                       >  3:224113/1‑149 (MQ=255)
                                                                          ttGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGCCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAACACCGCATTGTCACAAAGGTCg                                                                        >  3:253903/1‑148 (MQ=255)
                                                                             aaCGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAACACCGCATTGTCACAAAGGTCGAt                                                                      >  6:42123/1‑147 (MQ=255)
                                                                                gCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAACACCGCATTGTCACAAAGGTCGATCAGCt                                                                 >  2:303389/1‑149 (MQ=255)
                                                                                  tGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAACACCGCATTGTCACAAAGGTCGATCAGCTTa                                                               >  7:305604/1‑149 (MQ=255)
                                                                                  tGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAACACCGCATTGTCACAAAGGTCGATCAGCTTa                                                               >  6:349718/1‑149 (MQ=255)
                                                                                  tGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAACACCGCATTGTCACAAAGGTCGATCAGCTTa                                                               >  8:320568/1‑149 (MQ=255)
                                                                                    gagaAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAACACCGCATTGTCACAAAGGTCGATCAGCTTAt                                                              <  3:195708/148‑1 (MQ=255)
                                                                                     agaAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAACACCGCATTGTCACAAAGGTCGATCAGCTTATGa                                                            <  3:183091/149‑1 (MQ=255)
                                                                                      caaaCTGCCAAAAATCTCGTAAACGCCAACTACGGCAAGGCTGGAACTCATGTCAAAATCTCGCAAGAAAAGCTGACTTTGGCTCCGCTCCCGTTCAGCACCGCTCCGCGAGCAACACCGCATTGTCACAAAGGTCGATCAGCTTATGa                                                            <  1:263711/148‑1 (MQ=255)
                                                                                         aCTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAACACCGCATTGTCACAAAGGTCGATCAGCTTATGAAAc                                                         >  2:171594/1‑148 (MQ=255)
                                                                                              cAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAACACCGCATTGTCACAAAGGTCGATCAGCTTATGAAACTATGc                                                    >  5:157397/1‑148 (MQ=255)
                                                                                                aaTATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAACACCGCATTGTCACAAAGGTCGATCAGCTTATGAAACTa                                                       <  3:36746/143‑1 (MQ=255)
                                                                                                aaTATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATACCCTTAGCACCTCTCCGCGAGCAACACCCCATAGTCACAAAGGTCGATCAGCTAATCAAACTa                                                       >  4:36746/1‑143 (MQ=255)
                                                                                                 atatCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAACACCGCATTGTCACAAAGGTCGATCAGCTTATGAAACTATGCGATa                                                <  8:305604/149‑1 (MQ=255)
                                                                                                    tctcGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAACACCGCATTGTCACAAAGGTCGATCAGCTTATGAAACTATGCGATACg                                              <  2:315943/148‑1 (MQ=255)
                                                                                                        gAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAACACCGCATTGTCACAAAGGTCGATCAGCTTATGAAACTATGCGATACGc                                             <  4:171766/145‑1 (MQ=255)
                                                                                                          aCAGGCCAAGTCCGGCATGGCTGAAAGTCAGGTCAATATCTCTCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAACACCGCATTGTCACAAAGGTCGATCAGCTTATGAAACTATGCGATACGCTAAAGc                                       <  1:303389/149‑1 (MQ=255)
                                                                                                              gCCAAGTCCGGCATGGCTGCAAGTCAGGTGCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAACACCGCATTGTCACAAAGGTCGATCAGCTTATGAAACTATGCGATACGCTAAAGCaa                                     >  1:26236/1‑148 (MQ=255)
                                                                                                               ccAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAACACCGCATTGTCACAAAGGTCGATCAGCTTATGAAACTATGCGATACGCTAAAGCaa                                     >  5:346668/1‑146 (MQ=255)
                                                                                                                cAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAACACCGCATTGTCACAAAGGTCGATCAGCTTATGAAACTATGCGATACGCTAAAGCaa                                     >  5:112205/1‑145 (MQ=255)
                                                                                                                    tCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAACACCGCATTGTCACAAAGGTCGATCAGCTTATGAAACTATGCGATACGCTAAAGCAACAGATaaa                             <  1:281626/149‑1 (MQ=255)
                                                                                                                    tCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAACACCGCATTGTCACAAAGGTCGATCAGCTTATGAAACTATGCGATACGCTAAAGCAACAGATaaa                             <  7:320568/149‑1 (MQ=255)
                                                                                                                     ccGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAACACCGCATTGTCACAAAGGTCGATCAGCTTATGAAACTATGCGATACGCTAAAGCAACAGATAAAt                            >  2:230200/1‑149 (MQ=255)
                                                                                                                      cGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAACACCGCATTGTCACAAAGGTCGATCAGCTTATGAAACTATGCGATACGCTAAAGCAACAGATAAAt                            <  8:134488/148‑1 (MQ=255)
                                                                                                                      cGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAACACCGCATTGTCACAAAGGTCGATCAGCTTATGAAACTATGCGATACGCTAAAGCAACAGATAAAt                            <  7:2130/148‑1 (MQ=255)
                                                                                                                       ggCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAACACCGCATTGTCACAAAGGTCGATCAGCTTATGAAACTATGCGATACGCTAAAGCAACAGATAAATGt                          <  5:349718/149‑1 (MQ=255)
                                                                                                                        gCATGGCTGCAAGTCAGGTCAATATCTCACAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAACACCGCATTGTCACAAAGGTCGATCAGCTTATGAAACTATGCGATACGCTAAAGCAACAGATAAATGTa                         <  4:224113/149‑1 (MQ=255)
                                                                                                                                 cAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAACACCGCATTGTCACAAAGGTCGATCAGCTTATGAAACTATGCGATACGCTAAAGCAACAGATAAATGTAGCTCgca                  <  6:112205/147‑1 (MQ=255)
                                                                                                                                     tCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAACACCGCATTGTCACAAAGGTCGATCAGCTTATGAAACTATGCGATACGCTAAAGCAACAGATAAATGTAGCTc                     <  3:146105/140‑1 (MQ=255)
                                                                                                                                     tCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAACACCGCATTGTCACAAAGGTCGATCAGCTTATGAAACTATGCGATACGCTAAAGCAACAGATAAATGTAGCTc                     >  4:146105/1‑140 (MQ=255)
                                                                                                                                        ggTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAACACCGCATTGTCACAAAGGTCGATCAGCTTATGAAACTATGCGATACGCTAAAGCAACAGATAAATGTAGCTCGCAGCAAGCaaa         <  5:42123/149‑1 (MQ=255)
                                                                                                                                          tCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAACACCGCATTGTCACAAAGGTCGATCAGCTTATGAAACTATGCGATACGCTAAAGCAACAGATAAATGTAGCTCGCAGCAAGCAAACa       <  1:171594/149‑1 (MQ=255)
                                                                                                                                             aaTATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAACACCGCATTGTCACAAAGGTCGATCAGCTTATGAAACTATGCGATACGCTAAAGCAAc                                    >  7:274973/1‑118 (MQ=255)
                                                                                                                                             aaTATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAACACCGCATTGTCACAAAGGTCGATCAGCTTATGAAACTATGCGATACGCTAAAGCAAc                                    <  8:274973/118‑1 (MQ=255)
                                                                                                                                             aaTATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAACACCGCATTGTCACAAAGGTCGATCAGCTTATGAAACTATGCGATACGCTAAAGCAACAGATAAATGTAGCTCGCAGCAAGCAAACaga     <  6:157397/149‑1 (MQ=255)
                                                                                                                                                 tctcGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAACACCGCATTGTCACAAAGGTCGATCAGCTTATGAAACTATGCGATACGCTAAAGCAACAGATAAATGTAGCTCGCAGCAAGCAAACAGAGCt  <  6:346668/148‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
TAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAACACCGCATTGTCACAAAGGTCGATCAGCTTATGAAACTATGCGATACGCTAAAGCAACAGATAAATGTAGCTCGCAGCAAGCAAACAGAGCT  >  NC_002947/5398400‑5398691

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