Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,426,178 A→T 6.0% intergenic (+198/‑276) ydaY → / → ynaA Rac prophage; putative uncharacterized protein YdaY/Rac prophage; putative prophage tail length tape measure domain‑containing protein YnaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,426,1780AT 6.0% 242.9 / 6.8 95intergenic (+198/‑276)ydaY/ynaARac prophage; putative uncharacterized protein YdaY/Rac prophage; putative prophage tail length tape measure domain‑containing protein YnaA
Reads supporting (aligned to +/- strand):  ref base A (45/44);  new base T (4/2);  total (49/46)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.27e-01

AAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAG  >  NC_000913/1426037‑1426317
                                                                                                                                             |                                                                                                                                           
aaGCGGGAAAAACTCAGGTGGCCCTCGAGAGCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGAc                                                                                                                                            <  2:626304/143‑1 (MQ=255)
 aGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGGGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGGGCCTCTGCCTTTATCTCTGACa                                                                                                                                           >  1:668586/1‑143 (MQ=255)
 aGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACa                                                                                                                                           >  1:278261/1‑143 (MQ=255)
   cGAGAAAAACTCAGGTTGACACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACaga                                                                                                                                         <  1:290571/143‑1 (MQ=255)
    gagaAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACagag                                                                                                                                        >  1:767583/1‑143 (MQ=255)
      gaAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCGTAGGGACCTGTTGTACGCCTAAACGGCACGATGCCGTGCCCGACGATATCCTGGAATGGCTAGGGAGCCTTTGCCTTTATCTCTGACAGAGAt                                                                                                                                      >  2:293862/1‑143 (MQ=255)
       aaaaaCTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGAGATTAAGGAGTTGTTGTACGCCTATTCAGTAATATGCCGCGCCCGACGTTATGCTGGAATGACTGCAGTGCCATTGCCTACATCTCTGACagatgga                                                                                                                                     >  2:279311/1‑139 (MQ=255)
          accTCAGGTTGCCACCGATACCGGAACCAGAACTAATCCCTGTTCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAAt                                                                                                                                  <  2:593285/141‑1 (MQ=255)
           aCTCAGGTTGCCACCGATACCGGAACCAGAAATAGCCCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGTCAGAGATaaagg                                                                                                                                 <  1:557735/143‑3 (MQ=255)
             tCAGGTTGCCACCGATACCGGAACCAGAAATAGACCGGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGCTAGATGaa                                                                                                                               >  1:331830/1‑143 (MQ=255)
                 gTTGCCACCGATACCGGAACCAGAAAGAGACCCGGTGCTTAAGGAGTTGTTGTACGGCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATATCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATAtt                                                                                                                           <  2:256834/143‑1 (MQ=255)
                     ccACCGATACCGGAACCAGAAATAGACCCGGTGCCCAAGGAGTTGTTGTACGCCTACTCGGCAATATCCCGTGCCCGACGTGATGCCGGGATTGCTCGGGTGCCTTCACCTTTATCTCTGACAGAGATATATGAAGAtttcga                                                                                                                       >  2:352790/1‑140 (MQ=255)
                      cACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTCATATCCCGTGCCCGCCGTTGTGCTGGAATGGCTGGGGTGCCTGTGCCTAAATCGCTGACAGGGCTAACCGCACATTTAGcc                                                                                                                      >  2:148140/1‑143 (MQ=255)
                          gATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCACGTGCCCGACGTTATGCTGTACTGGCTGGACAGCCTTTGCCTTTTTTTCTGACAGAACTAAATGAATATCTAGCCACta                                                                                                                  >  2:357326/1‑142 (MQ=255)
                            tACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACACAt                                                                                                                >  2:659165/1‑143 (MQ=255)
                               cGGAACCAGAAATAGACCCGGTGCTTGAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTTTCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCg                                                                                                             >  2:806496/1‑143 (MQ=255)
                               cGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTGTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCg                                                                                                             >  2:34783/1‑143 (MQ=255)
                                  aaCCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTa                                                                                                          >  2:661834/1‑143 (MQ=255)
                                      aGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAAGCGTAAGGCACTCATCCGGTattga                                                                                                      >  2:546686/1‑143 (MQ=255)
                                       ggaaTAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTattgat                                                                                                     <  2:425779/141‑1 (MQ=255)
                                       gAAATCGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTGATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTattgat                                                                                                     <  1:132656/143‑1 (MQ=255)
                                            aGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTTGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAGTGAATAGTTAGCCACTCATCCGGTATTGATTGAga                                                                                                >  1:315006/1‑142 (MQ=255)
                                            aGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGAGAGAGATCAAGGAAGATTTAGCCACTCATCCGGTATTGATTGAgc                                                                                                >  2:502242/1‑143 (MQ=255)
                                               ccgggTGCTTAGGGAGTTGTTGTACGCCTATGCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGcgcg                                                                                             <  2:835271/140‑1 (MQ=255)
                                                ccGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGATATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGAATTGATTGAGCTCGa                                                                                            >  2:140672/1‑143 (MQ=255)
                                                        tAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTAGCCTTTATCTCTGACAGAGATAAATGAAGATTTAGCCACTCATCCGGTATTGCTTCAGCGGGAAGAATTTg                                                                                    >  2:543113/1‑143 (MQ=255)
                                                          gggAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGaa                                                                                  <  1:122368/142‑1 (MQ=255)
                                                          aGGAGTTGTTGTACGCCTATTCGGTACTATCCCGTGCCCGACGTGATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGaa                                                                                  <  1:594924/143‑1 (MQ=255)
                                                          aGGAGTTGTGGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGaa                                                                                  <  1:497908/143‑1 (MQ=255)
                                                           ggcgatcttgtACGCCTCTTCGGTACTATCCCGTGCCCCACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAg                                                                                 <  1:54015/136‑1 (MQ=255)
                                                           ggAGTTGTTGTTCGCCTATTGGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAg                                                                                 <  2:719923/143‑1 (MQ=255)
                                                           ggAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAg                                                                                 >  2:3417/1‑143 (MQ=255)
                                                           ggAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAg                                                                                 <  1:514508/143‑1 (MQ=255)
                                                              gctgttgtACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAg                                                                              <  1:218654/141‑1 (MQ=255)
                                                                 gttgtACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGa                                                                           >  1:762325/1‑143 (MQ=255)
                                                                    gtACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAAGGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCt                                                                        >  2:299717/1‑143 (MQ=255)
                                                                     tACGCCTATTCGGTAATACCCCGTGCCCGACGTCATGCTGGAACGGCTGGGGTGCCTTTGACTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCtt                                                                       <  1:216244/143‑1 (MQ=255)
                                                                      aCGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGCATATTTAGCCACTCATCCGGTCTTTATTGAGCGCGATGAATTTGAAGCAGTGCTCttc                                                                      >  1:552654/1‑142 (MQ=255)
                                                                      aCGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCttt                                                                      <  1:200091/143‑1 (MQ=255)
                                                                         ccTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCa                                                                   >  1:264646/1‑143 (MQ=255)
                                                                          cTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGACGCCGCGATCTTTGCAc                                                                  >  2:321334/1‑143 (MQ=255)
                                                                           tATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACt                                                                 <  1:596252/143‑1 (MQ=255)
                                                                           tATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACt                                                                 <  2:316019/143‑1 (MQ=255)
                                                                                 ggattatCCCGTGCCCGGCGTTATGCTGACATGGCTGGGGTGCCTTGGCCGTTATCTCTGACAGAGATAAATGAATAGTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGa                                                           <  2:427262/139‑1 (MQ=255)
                                                                                 gTAATATCCCGTGCTCGACGTTATGCTGGCATGGCTGGGGTGCCTTTGCCTTTATCTCTGTCAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGa                                                           <  2:398293/143‑1 (MQ=255)
                                                                                 gTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAAGGAATATTTAGCAACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCCTTGCACTGGATGa                                                           <  1:710790/143‑1 (MQ=255)
                                                                                      atCCCGTGCCCGACGTTATTCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGt                                                      <  1:291690/143‑1 (MQ=255)
                                                                                       tCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACGAGTa                                                     >  1:677288/1‑143 (MQ=255)
                                                                                         ccGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTAtt                                                   <  2:120043/143‑1 (MQ=255)
                                                                                          cGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGAAGAATTTGAAGCCGCGATCTTTGCACTGGATGACCAGTAttt                                                  >  2:310464/1‑143 (MQ=255)
                                                                                          cGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCCCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTAttt                                                  <  1:434373/143‑1 (MQ=255)
                                                                                            tgcgcttcGTTATGCTGCAGTGCCTGGGTCTCCTTTGCCTTTATCTCTGAAGGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTc                                                <  2:488920/136‑1 (MQ=255)
                                                                                            tGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTc                                                <  2:161910/143‑1 (MQ=255)
                                                                                             gCCCGACGTTATGCTGGAATGGCTGGGGGGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCa                                               >  2:847608/1‑143 (MQ=255)
                                                                                              ccgggagTTACGCTGGAATGGCTGGGTTGCTTTTGCCTGTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAg                                              <  2:78869/137‑1 (MQ=255)
                                                                                              cccGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTGTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAg                                              <  1:374528/143‑1 (MQ=255)
                                                                                                cGACGTGATGCTGTAGTGGCGGGGGTGCCTTTGCCCTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCAGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGa                                            <  2:835915/143‑1 (MQ=255)
                                                                                                 gACGTTATGCTGGAATGGCTGGGGAGCCTTTGCCTTTATCTCTGACAGAGATAAATGAGTATTTAGCCAGTCATCCGGTATTGATTGAGCGCGATGAAATTGAAGCAGTGATCTTGGCACTGGATGACCAGTATTTGCAGGAg                                           >  2:90667/1‑143 (MQ=255)
                                                                                                   cGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCa                                         >  1:98823/1‑143 (MQ=255)
                                                                                                   cGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCa                                         <  2:278261/143‑1 (MQ=255)
                                                                                                       ggggtGGGATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATGGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAgtgt                                     <  1:218050/139‑1 (MQ=255)
                                                                                                        tGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAgtgtg                                    >  1:605218/1‑143 (MQ=255)
                                                                                                         gCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGCCCAGTATTTTCAGCACCAgtgtgt                                   >  2:152435/1‑143 (MQ=255)
                                                                                                          cTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGCGCAgtgtgtg                                  >  1:105213/1‑143 (MQ=255)
                                                                                                           tGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTCTCAGGCGCActctgtgt                                 >  2:768879/1‑143 (MQ=255)
                                                                                                            ggtaTGGGTGGTGTGGCTTTGCCTTTATCGCTGACGGGGATAAATGAATATTTGGCCACTCGTCCGGTATTGATTGAGCGCGATGAAGTTGAAGCAGTGATCTGTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTa                                <  2:287481/140‑1 (MQ=255)
                                                                                                            ggAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTa                                >  1:408923/1‑143 (MQ=255)
                                                                                                              aaTGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAgt                              >  1:66096/1‑143 (MQ=255)
                                                                                                               aTGGCGGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAgtt                             <  1:379234/143‑1 (MQ=255)
                                                                                                                 ggCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAgttgt                           >  1:552977/1‑143 (MQ=255)
                                                                                                                    tGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTaa                        <  1:396899/143‑1 (MQ=255)
                                                                                                                      gggTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGACGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAAtt                      <  1:158505/143‑1 (MQ=255)
                                                                                                                      gggTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCCATGAATTTCAAGCAGTCATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTCGTTGTTAAtt                      >  2:688670/1‑143 (MQ=255)
                                                                                                                        gTGCCATTGCCTTTATCGGGGACAGAGATAAATGGAGATCGAGCCACTCAGCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTAc                    <  1:127705/143‑1 (MQ=255)
                                                                                                                         tgcctttgccttTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACg                   >  1:696563/1‑143 (MQ=255)
                                                                                                                          gcctttgcctttATCTCTTACAGAGATATATTAATATTTAGCTACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGt                  <  1:680622/143‑1 (MQ=255)
                                                                                                                          gcctttgcctttATCTCTGATAAACATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGt                  <  2:743006/143‑1 (MQ=255)
                                                                                                                           cctttgcctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTa                 <  2:264646/143‑1 (MQ=255)
                                                                                                                             tttgcctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGCATTTCCAGGAGCAGTGTGTGTAGTTGTTATTTACGTaca               >  2:630181/1‑143 (MQ=255)
                                                                                                                               tgcctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTGCAct             >  2:71022/1‑143 (MQ=255)
                                                                                                                                 catttATCTCTGTCAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACActct           >  1:721451/3‑143 (MQ=255)
                                                                                                                                  atttATCTCTGTCAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCTGTGTGTGTAGTTGTTAATTACGTACACTCTg          >  2:793538/2‑143 (MQ=255)
                                                                                                                                  atttATCTCTGTCAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACtgtg          >  2:304836/2‑143 (MQ=255)
                                                                                                                                  atttATCTCTGTCAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGGGTGTGTAGTTGTTAATTACGTACACTCTg          >  2:405607/2‑143 (MQ=255)
                                                                                                                                   tttATGTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGt         <  1:802018/143‑1 (MQ=255)
                                                                                                                                   tttATCTCTGACAGAGATAAAAGCCTGTTTAGCCACTCATCCGGTATTGCTTGCGTGATACGGAGGCGATGCTTCCTTCTTTTCACTGTGTGACCAGTATTTTCAGGCGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGt         >  2:598974/1‑143 (MQ=255)
                                                                                                                                   tttATCTCCGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGt         >  2:463159/1‑143 (MQ=255)
                                                                                                                                   cgtATCTTCGACAGACATAAATCAATATTTAGCCAATCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTACTCGCACTAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGt         <  2:784256/141‑1 (MQ=255)
                                                                                                                                    ctATCACTGACAGATATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGtt        <  1:473334/142‑1 (MQ=255)
                                                                                                                                       tctcTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACa     <  1:539450/143‑1 (MQ=255)
                                                                                                                                       tctcTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACa     <  2:745340/143‑1 (MQ=255)
                                                                                                                                        ctctGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACaa    >  1:65086/1‑142 (MQ=255)
                                                                                                                                         tctGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACaga   >  1:457010/1‑143 (MQ=255)
                                                                                                                                         tctAACAGATATAAATGAAGATTTAGCCACTCATCCGGTAGTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTGCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACaga   <  2:363201/143‑1 (MQ=255)
                                                                                                                                          ctGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACagag  >  1:687662/1‑143 (MQ=255)
                                                                                                                                             |                                                                                                                                           
AAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAG  >  NC_000913/1426037‑1426317

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