Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,509,750 A→G pseudogene (272/315 nt) insO → KpLE2 phage‑like element; IS911B regulator fragment

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,509,7500AG98.6% 244.8 / NA 70pseudogene (272/315 nt)insOKpLE2 phage‑like element; IS911B regulator fragment
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (45/23);  total (46/24)
Fisher's exact test for biased strand distribution p-value = 3.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAATGGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAACCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCACGATGACTC  >  NC_000913/4509616‑4509882
                                                                                                                                      |                                                                                                                                    
cGAATGGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAGCCGAt                                                                                                                                 >  2:11169/1‑140 (MQ=255)
  aaTGGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTAAGACCGCACGAATATAATGGTGGGTTACCACGAAACGAATCAGAAAGCCGATAc                                                                                                                               >  2:203962/1‑140 (MQ=255)
     gggTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTAAGACCGCACGAATATAATGGTGGGTTACCACAAATCGAATCCGAAAGCCGATACTgg                                                                                                                            >  1:252027/1‑140 (MQ=255)
      ggTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTAAGACCGCACGAATATAATAGTGGGTTTCCACCAAACGAATTAGAAAGCCGATACTGGa                                                                                                                           >  2:291542/1‑140 (MQ=255)
           cAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAGCCGATACTGGAAAAAc                                                                                                                      >  1:447297/1‑140 (MQ=255)
                 cGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAGCCGATACTGGAAAAACTCTAAc                                                                                                                >  2:508279/1‑140 (MQ=255)
                       aTGTAAGCTTCAGCGGTGCAGCGCACGCAGTAACGGAGTATATCGTTGGATATTACAGCGCACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAg                                                                                                          <  1:291542/140‑1 (MQ=255)
                        tGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAgg                                                                                                         <  2:573974/140‑1 (MQ=255)
                        tGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAgg                                                                                                         <  1:508279/140‑1 (MQ=255)
                        tGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAgg                                                                                                         <  2:31622/140‑1 (MQ=255)
                          tAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGcc                                                                                                       >  2:29629/1‑140 (MQ=255)
                              cTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCATAAACCCGATACTGGAAAAACTTTAACGCGGAGGCCAGtt                                                                                                   >  2:261900/1‑140 (MQ=255)
                               ttCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGttt                                                                                                  >  2:657432/1‑140 (MQ=255)
                                  aGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAGCGGATACTGTAAAAACACTAACGGGGAGGCCAGTTTTAg                                                                                               >  1:209708/1‑140 (MQ=255)
                                  aGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGTGGAGGCCAGTTTTAg                                                                                               >  1:124755/1‑140 (MQ=255)
                                  aGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAg                                                                                               >  1:58879/1‑140 (MQ=255)
                                  aGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAg                                                                                               >  1:409557/1‑140 (MQ=255)
                                  aGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAGCCGATACTGAAAAAACTCTAACGCGGAGGCCAGTTTTAg                                                                                               >  1:299574/1‑140 (MQ=255)
                                    cGATGCAGCTCACGCAACAACGTACTATATCCTTGGATAGTACAGCGCTCTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGtt                                                                                             <  2:252027/140‑1 (MQ=255)
                                       tGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGAc                                                                                          <  1:8438/140‑1 (MQ=255)
                                       tGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGAc                                                                                          <  1:115857/140‑1 (MQ=255)
                                                             tataTCGTTGGATATTACAGCGCACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGa                                                                    >  1:100288/1‑140 (MQ=255)
                                                             tataTCGTTGGATATTACAGCGCACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGa                                                                    >  2:518846/1‑140 (MQ=255)
                                                               tatCGTTGGATATTACAGCGCACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGc                                                                  >  1:160752/1‑140 (MQ=255)
                                                                  cGTTGGATATTACAGCGCACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCcgc                                                               <  1:165324/140‑1 (MQ=255)
                                                                     tGGATATTACAGCGCACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGt                                                            >  2:399130/1‑140 (MQ=255)
                                                                     tGGATATTACAGCGCACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGt                                                            >  1:500705/1‑140 (MQ=255)
                                                                       gATATTACAGCGCACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAATTGACCAGAACCGATTACCTGAAGGGAGCCGCGGTCg                                                          >  1:275554/1‑140 (MQ=255)
                                                                         tatTACAGCGCACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGcc                                                        >  1:575989/1‑140 (MQ=255)
                                                                         tatTACAGCGCACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGcc                                                        >  1:583715/1‑140 (MQ=255)
                                                                           ttACAGCGCACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTg                                                      <  2:100288/140‑1 (MQ=255)
                                                                                 cgcACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGtt                                                >  1:707771/1‑140 (MQ=255)
                                                                                 cgcACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGtt                                                >  1:426425/1‑140 (MQ=255)
                                                                                 cgcACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGtt                                                >  2:405501/1‑140 (MQ=255)
                                                                                 cgcACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGtt                                                >  2:412372/1‑140 (MQ=255)
                                                                                   cACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGc                                              >  2:706517/1‑140 (MQ=255)
                                                                                   cACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGc                                              >  1:498316/1‑140 (MQ=255)
                                                                                   cACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGc                                              >  2:104364/1‑140 (MQ=255)
                                                                                      tAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGt                                           >  2:615917/1‑140 (MQ=255)
                                                                                        aGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAg                                         <  2:500705/140‑1 (MQ=255)
                                                                                              cACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGgag                                   >  1:429864/1‑140 (MQ=255)
                                                                                                cGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGat                                 >  2:580095/1‑140 (MQ=255)
                                                                                                     ataATGGTGGGTTACCACCAAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATcc                            <  1:518846/140‑1 (MQ=255)
                                                                                                        aTGGTGGGTTACCACCAAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGTGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGc                         <  2:124755/140‑1 (MQ=255)
                                                                                                          ggtggGTTACCACCAAAAGAATAAGAAAGCCTATACTGGTAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCg                       <  2:300788/140‑1 (MQ=255)
                                                                                                           gtggGTTACCACCAAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGc                      >  1:41859/1‑140 (MQ=255)
                                                                                                             gggTTACCACCAAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGATATCATCCAGCCGCAg                    >  1:320938/1‑140 (MQ=255)
                                                                                                             gggTTACCACCAAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATcc                            >  2:388107/1‑132 (MQ=255)
                                                                                                             gggTTACCACCAAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATcc                            <  1:388107/132‑1 (MQ=255)
                                                                                                             gggTTACCACCAAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAg                    >  1:577892/1‑140 (MQ=255)
                                                                                                              ggTTACCACCAAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCTCCTGGCAATTGCAGTAACAGGTGATCATCCAGCCGCAGt                   >  2:224055/1‑140 (MQ=255)
                                                                                                              ggTTACCACCAAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCGTGGCAGTGGCAGTAGCAGGAGCTCATCCAGCCGCAGt                   >  1:300788/1‑140 (MQ=255)
                                                                                                              ggTTACCACCAAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGt                   >  2:713389/1‑140 (MQ=255)
                                                                                                              ggTTACCACCAAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGt                   >  2:139410/1‑140 (MQ=255)
                                                                                                                ttACCTCCAAACGAATCCGAAAACCGATACTAGAATAACTCTAACGCGTAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTAc                 <  2:209708/140‑1 (MQ=255)
                                                                                                                ttACCACCAAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTAc                 >  2:177618/1‑140 (MQ=255)
                                                                                                                  accaccAAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAg               <  1:139410/140‑1 (MQ=255)
                                                                                                                   ccaccaAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAgg              >  1:86986/1‑140 (MQ=255)
                                                                                                                    caccaAACGAATCAGAACGCCGATACTGGAATAACCCTATCGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGt             <  2:299574/140‑1 (MQ=255)
                                                                                                                      ccaAACGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCa           >  2:152028/1‑140 (MQ=255)
                                                                                                                        aaaCGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCAca         >  2:671048/1‑139 (MQ=255)
                                                                                                                        aaaCGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCACg         >  1:77819/1‑140 (MQ=255)
                                                                                                                         aaCGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCACGa        <  2:577892/140‑1 (MQ=255)
                                                                                                                         aaCGAATCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGAGCATCCAGCCACAGTACAGGTCACGa        >  1:224303/1‑140 (MQ=255)
                                                                                                                              aTCAGAAAGCCGATACTGGAAATACCCTATCGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCACGATGAct   <  1:203962/140‑1 (MQ=255)
                                                                                                                              aTCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCACGATGAct   <  2:160752/140‑1 (MQ=255)
                                                                                                                              aTCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCACGATGAct   <  1:706517/140‑1 (MQ=255)
                                                                                                                              aTCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCACGATGAct   <  1:657432/140‑1 (MQ=255)
                                                                                                                              aTCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCACGATGAct   <  1:29629/140‑1 (MQ=255)
                                                                                                                               tCAGAAAGCCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCACGATGActc  <  2:447297/140‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                    
CGAATGGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTAAGACCGCACGAATATAATGGTGGGTTACCACCAAACGAATCAGAAAACCGATACTGGAAAAACTCTAACGCGGAGGCCAGTTTTAGTTGACCACAACAGACTACCTGAAGGGAGCCGCGGTCGCCTGGCAGTTGCAGTAGCAGGAGATCATCCAGCCGCAGTACAGGTCACGATGACTC  >  NC_000913/4509616‑4509882

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