Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,426,178 A→T 10.3% intergenic (+96/‑276) ydaY → / → tmpR pseudogene, Rac prophage;Phage or Prophage Related/Rac prophage; pseudogene, tail protein family;Phage or Prophage Related; putative alpha helix protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,426,1780AT10.3% 280.4 / 20.7 89intergenic (+96/‑276)ydaY/tmpRpseudogene, Rac prophage;Phage or Prophage Related/Rac prophage; pseudogene, tail protein family;Phage or Prophage Related; putative alpha helix protein
Reads supporting (aligned to +/- strand):  ref base A (39/39);  new base T (5/4);  total (45/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.07e-01

TGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAG  >  NC_000913/1426034‑1426317
                                                                                                                                                |                                                                                                                                           
tGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACag                                                                                                                                          >  4:209254/1‑148 (MQ=255)
      cGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATaa                                                                                                                                    <  3:1062493/148‑1 (MQ=255)
          aaaaaCTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACagag                                                                                                                                        <  2:878481/140‑1 (MQ=255)
          aaaaaCTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACagag                                                                                                                                        >  1:878481/1‑140 (MQ=255)
          aaaaaCTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGa                                                                                                                                >  4:690991/1‑148 (MQ=255)
                tCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATAttt                                                                                                                          >  4:295574/1‑148 (MQ=255)
                           ccGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGTCAgaga                                                                                                                                       >  2:735387/1‑124 (MQ=255)
                           ccGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGTCAgaga                                                                                                                                       <  1:735387/124‑1 (MQ=255)
                                 ccGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATaaa                                                                                                                                   >  4:932020/1‑122 (MQ=255)
                                 ccGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATaaa                                                                                                                                   <  3:932020/122‑1 (MQ=255)
                                  cGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTatt                                                                                                        >  3:531493/1‑148 (MQ=255)
                                    gAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTattgn                                                                                                      <  3:209254/148‑2 (MQ=255)
                                          gAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACag                                                                                                                                          <  4:712618/106‑1 (MQ=255)
                                          gAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACag                                                                                                                                          >  3:712618/1‑106 (MQ=255)
                                              tAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATaaa                                                                                                                                   <  1:633219/109‑1 (MQ=255)
                                              tAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATaaa                                                                                                                                   >  2:633219/1‑109 (MQ=255)
                                                  cccGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGaa                                                                                        >  4:190824/1‑148 (MQ=255)
                                                    cGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGGCAGAGATAAATGAATATTTAGCCACt                                                                                                                   <  1:847426/119‑1 (MQ=255)
                                                    cGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGGCAGAGATAAATGAATATTTAGCCACt                                                                                                                   >  2:847426/1‑119 (MQ=255)
                                                    cGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAAtt                                                                                      >  1:202240/1‑148 (MQ=255)
                                                           tAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCa                                                                               <  3:295574/148‑1 (MQ=255)
                                                              ggAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTg                                                                            <  3:690991/148‑1 (MQ=255)
                                                                   tgttgtACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCtt                                                                       <  4:531493/148‑1 (MQ=255)
                                                                    gttgtACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCttt                                                                      <  3:190824/148‑1 (MQ=255)
                                                                      tgtACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGc                                                                    >  3:970231/1‑148 (MQ=255)
                                                                        tACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAgaga                                                                                                                                       <  4:106458/79‑1 (MQ=255)
                                                                        tACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAgaga                                                                                                                                       >  3:106458/1‑79 (MQ=255)
                                                                           gCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTgg                                                               <  2:202240/148‑1 (MQ=255)
                                                                            ccTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTgg                                                               >  1:461482/1‑147 (MQ=255)
                                                                            ccTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTgg                                                               <  2:461482/147‑1 (MQ=255)
                                                                               aTTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCa                                                                   <  4:716012/140‑1 (MQ=255)
                                                                               aTTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCa                                                                   >  3:716012/1‑140 (MQ=255)
                                                                                 tCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGa                                                                                   <  1:9048/122‑1 (MQ=255)
                                                                                 tCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGa                                                                                   >  2:9048/1‑122 (MQ=255)
                                                                                     tAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGTCAGAGATaaa                                                                                                                                   <  1:752257/70‑1 (MQ=255)
                                                                                     tAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGTCAGAGATaaa                                                                                                                                   >  2:752257/1‑70 (MQ=255)
                                                                                           cccGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTattgattg                                                                                                   >  2:661086/1‑96 (MQ=255)
                                                                                           cccGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTattgattg                                                                                                   <  1:661086/96‑1 (MQ=255)
                                                                                           cccGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTc                                                >  1:311900/1‑147 (MQ=255)
                                                                                           cccGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTc                                                <  2:311900/147‑1 (MQ=255)
                                                                                                     aCGTTATGCTGGAATGGGTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCa                                                                   >  4:916688/1‑118 (MQ=255)
                                                                                                     aCGTTATGCTGGAATGGGTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCa                                                                   <  3:916688/118‑1 (MQ=255)
                                                                                                     aCGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTg                                                                <  3:682860/121‑1 (MQ=255)
                                                                                                     aCGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTg                                                                >  4:682860/1‑121 (MQ=255)
                                                                                                      cGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAgtgtg                                    <  4:970231/148‑1 (MQ=255)
                                                                                                             cTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGc                                                                                <  3:35577/97‑1 (MQ=255)
                                                                                                             cTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGc                                                                                >  4:35577/1‑97 (MQ=255)
                                                                                                               ggAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGaa                                                                                  >  3:327731/1‑93 (MQ=255)
                                                                                                               ggAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGaa                                                                                  <  4:327731/93‑1 (MQ=255)
                                                                                                                gAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATaa                                                                                                                                    >  4:961538/1‑42 (MQ=255)
                                                                                                                gAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATaa                                                                                                                                    <  3:961538/42‑1 (MQ=255)
                                                                                                                gAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACt                                                                 >  4:75388/1‑109 (MQ=255)
                                                                                                                gAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACt                                                                 <  3:75388/109‑1 (MQ=255)
                                                                                                                gAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAgttgt                           >  4:997947/1‑147 (MQ=255)
                                                                                                                gAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAgttgt                           <  3:997947/147‑1 (MQ=255)
                                                                                                                      cTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTattgat                                                                                                     >  4:642718/1‑67 (MQ=255)
                                                                                                                      cTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTattgat                                                                                                     <  3:642718/67‑1 (MQ=255)
                                                                                                                           gTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATGTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTaa                        <  2:334038/139‑1 (MQ=255)
                                                                                                                           gTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTaa                        >  1:334038/1‑139 (MQ=255)
                                                                                                                                    cctttATCTCTGACAGAGATAAATGAATGTTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTAc      >  3:727437/1‑148 (MQ=255)
                                                                                                                                    cctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGt         >  1:751707/1‑145 (MQ=255)
                                                                                                                                    cctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGt         <  2:751707/145‑1 (MQ=255)
                                                                                                                                      tttATCTCTGTCAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCa                                                                               >  2:5150/1‑73 (MQ=255)
                                                                                                                                      tttATCTCTGTCAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCa                                                                               <  1:5150/73‑1 (MQ=255)
                                                                                                                                      tttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGaa                                                                                        >  3:700487/1‑64 (MQ=255)
                                                                                                                                      tttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCa                                                                               <  2:395529/73‑1 (MQ=255)
                                                                                                                                      tttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCa                                                                               >  1:395529/1‑73 (MQ=255)
                                                                                                                                      tttATCTCTGACAGAGATAAATGAATATTTAGCCACCCATCCGGTATTGATTGAGCGCGATGaa                                                                                        <  4:700487/64‑1 (MQ=255)
                                                                                                                                       ttATCTCTGTCAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGAcc                                                         >  4:832639/1‑94 (MQ=255)
                                                                                                                                       ttATCTCTGTCAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGAcc                                                         <  3:832639/94‑1 (MQ=255)
                                                                                                                                       ttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAgttgt                           >  4:73985/1‑124 (MQ=255)
                                                                                                                                       ttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGCATTTTCAGGAGCAGTGTGTGTAgttgt                           <  3:73985/124‑1 (MQ=255)
                                                                                                                                        tATCTCTGTCAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACagag  >  2:590701/1‑148 (MQ=255)
                                                                                                                                        tATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCtt                                                                       <  4:178515/79‑1 (MQ=255)
                                                                                                                                        tATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCtt                                                                       >  3:178515/1‑79 (MQ=255)
                                                                                                                                        tATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGt                                                      <  3:574226/96‑1 (MQ=255)
                                                                                                                                        tATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGt                                                      >  4:574226/1‑96 (MQ=255)
                                                                                                                                            tctGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTac                <  1:185225/130‑1 (MQ=255)
                                                                                                                                            tctGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTac                >  2:185225/1‑130 (MQ=255)
                                                                                                                                             ctGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGc                                                                                >  4:554828/1‑65 (MQ=255)
                                                                                                                                             ctGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGc                                                                                <  3:554828/65‑1 (MQ=255)
                                                                                                                                             ctGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAg                                                                              >  4:16927/1‑67 (MQ=255)
                                                                                                                                             ctGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAg                                                                              <  3:16927/67‑1 (MQ=255)
                                                                                                                                             ctGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAgg                                             <  1:900083/100‑1 (MQ=255)
                                                                                                                                             ctGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAgg                                             >  2:900083/1‑100 (MQ=255)
                                                                                                                                               gACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAgtgtg                                    >  1:759355/1‑107 (MQ=255)
                                                                                                                                               gACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAgtgtg                                    <  2:759355/107‑1 (MQ=255)
                                                                                                                                               gACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAgtgt                                     <  1:913687/106‑1 (MQ=255)
                                                                                                                                               gACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAgtgt                                     >  2:913687/1‑106 (MQ=255)
                                                                                                                                                |                                                                                                                                           
TGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAG  >  NC_000913/1426034‑1426317

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