Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,353,170 G→A 8.1% intergenic (+89/+30) ghoT → / ← lysU toxin of the GhoTS toxin‑antitoxin system/lysine‑‑tRNA ligase/Ap4A synthetase/Ap3A synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,353,1700GA 8.1% 411.7 / 17.6 138intergenic (+89/+30)ghoT/lysUtoxin of the GhoTS toxin‑antitoxin system/lysine‑‑tRNA ligase/Ap4A synthetase/Ap3A synthetase
Reads supporting (aligned to +/- strand):  ref base G (63/61);  new base A (5/6);  total (68/70)
Fisher's exact test for biased strand distribution p-value = 7.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.31e-01

Only 100 of 140 total aligned reads displayed.

CCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTGGCGGCAGACCATA  >  NC_000913/4353021‑4353319
                                                                                                                                                     |                                                                                                                                                     
ccTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCg                                                                                                                                                       >  2:108812/1‑150 (MQ=255)
       ggAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACt                                                                                                                                                >  1:104435/1‑150 (MQ=255)
       ggAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACt                                                                                                                                                <  1:1147668/150‑1 (MQ=255)
        gAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTg                                                                                                                                               >  1:363532/1‑150 (MQ=255)
        gAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTg                                                                                                                                               <  2:1171212/150‑1 (MQ=255)
        gAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTCACTg                                                                                                                                               >  2:715362/1‑150 (MQ=255)
            aaCCTGGCCTATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTc                                                                                                                                           <  2:106011/150‑1 (MQ=255)
             aCCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCAGGCGACTGCTTCg                                                                                                                                          <  2:870038/150‑1 (MQ=255)
              ccTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGt                                                                                                                                         <  2:642583/150‑1 (MQ=255)
               cTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGtt                                                                                                                                        >  1:89253/1‑150 (MQ=255)
               cTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCAGGCGACTGCTTCGtt                                                                                                                                        <  1:92330/150‑1 (MQ=255)
                tGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTc                                                                                                                                       <  1:333084/150‑1 (MQ=255)
                tGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTc                                                                                                                                       <  1:1106811/150‑1 (MQ=255)
                tGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCAGGCGACTGCTTCGTTc                                                                                                                                       <  2:1016689/150‑1 (MQ=255)
                 ggCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCa                                                                                                                                      >  1:453404/1‑150 (MQ=255)
                   ccAATGAGTCTGCCTGGGGCGTTACTTTTTTCTCTCGGTTAGGAGCTGTAGTTGAGTGGTCCCGGTGGTTAGCGTGAAACGAAGCGGCGCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTGAtt                                                                                                                                    >  2:661785/1‑150 (MQ=255)
                    cAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTa                                                                                                                                   <  2:1601826/150‑1 (MQ=255)
                      aTGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTaaa                                                                                                                                 >  1:1000676/1‑150 (MQ=255)
                        gAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGt                                                                                                                               <  2:57923/150‑1 (MQ=255)
                         aGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTg                                                                                                                              <  2:629700/150‑1 (MQ=255)
                            cTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGaaa                                                                                                                           >  2:1556117/1‑150 (MQ=255)
                            cTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGCATTTGCGGAACGATGAAAGAACCAGTCGCCTGACTGCTTCGTTCATTAAAGTGaaa                                                                                                                           >  2:1560329/1‑150 (MQ=255)
                              gCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAatt                                                                                                                         >  1:29564/1‑150 (MQ=255)
                              gCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAatt                                                                                                                         >  1:609435/1‑150 (MQ=255)
                                 tgtgGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAAAGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAatttat                                                                                                                      >  2:1537453/1‑150 (MQ=255)
                                  gtgGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAatttatt                                                                                                                     <  2:1476987/150‑1 (MQ=255)
                                      cATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCtg                                                                                                                 <  1:98360/150‑1 (MQ=255)
                                       aTTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCtgt                                                                                                                <  2:488374/150‑1 (MQ=255)
                                            ttttttCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGgcg                                                                                                           <  1:550442/150‑1 (MQ=255)
                                             tttttCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGgcgc                                                                                                          >  2:48034/1‑150 (MQ=255)
                                              ttttCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCa                                                                                                         <  1:707355/150‑1 (MQ=255)
                                                  ctctctTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGc                                                                                                     >  1:508961/1‑150 (MQ=255)
                                                     tctTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCgg                                                                                                  >  2:500465/1‑150 (MQ=255)
                                                     cctTTAAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAATACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCgg                                                                                                  <  2:873938/149‑1 (MQ=255)
                                                      ctTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCggg                                                                                                 >  2:350869/1‑150 (MQ=255)
                                                      ctTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCggg                                                                                                 >  1:70133/1‑150 (MQ=255)
                                                      ctTTTAGGAGCTGTAGTTGAGTCGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAATCAGTCGCATGACTGCTTCGTTAATTTAAGTGACATTTATTTCTGTGGGAGCATCGCCggg                                                                                                 >  1:1526741/1‑150 (MQ=255)
                                                         ttAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAag                                                                                              <  1:611432/150‑1 (MQ=255)
                                                         ttAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAag                                                                                              <  1:246836/150‑1 (MQ=255)
                                                          tAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAaga                                                                                             >  2:320708/1‑150 (MQ=255)
                                                          tAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAaga                                                                                             >  1:783461/1‑150 (MQ=255)
                                                           aGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAagag                                                                                            >  2:422600/1‑150 (MQ=255)
                                                           aGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAagag                                                                                            >  1:95452/1‑150 (MQ=255)
                                                            ggAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCTCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGAATCGCCGGGAAgaga                                                                                           <  2:794297/150‑1 (MQ=255)
                                                            ggAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAgaga                                                                                           <  2:111026/150‑1 (MQ=255)
                                                             gAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGaa                                                                                          <  2:1318227/150‑1 (MQ=255)
                                                              aGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGACTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAAt                                                                                         <  1:406100/150‑1 (MQ=255)
                                                               gCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATa                                                                                        >  2:745918/1‑150 (MQ=255)
                                                                 tGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTg                                                                                                                                                   <  2:1254391/89‑1 (MQ=255)
                                                                 tGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTg                                                                                                                                                   >  1:1254391/1‑89 (MQ=255)
                                                                  gTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACg                                                                                     <  2:1317790/150‑1 (MQ=255)
                                                                      ttGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCTCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTcgc                                                                                 <  1:565260/150‑1 (MQ=255)
                                                                      ttGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTcgc                                                                                 <  1:1260090/150‑1 (MQ=255)
                                                                       tGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTcgcg                                                                                <  2:159478/150‑1 (MQ=255)
                                                                       tGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTcgcg                                                                                >  2:252989/1‑150 (MQ=255)
                                                                       tGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTcgcg                                                                                >  1:381173/1‑150 (MQ=255)
                                                                          gTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAAt                                                                             <  1:350869/150‑1 (MQ=255)
                                                                             gACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGt                                                                          >  1:1487571/1‑150 (MQ=255)
                                                                             gACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGt                                                                          >  2:1118634/1‑150 (MQ=255)
                                                                              aCCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTa                                                                         >  1:964125/1‑150 (MQ=255)
                                                                                 cGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATgg                                                                      >  1:846051/1‑150 (MQ=255)
                                                                                        tATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAg                                                               >  1:1349034/1‑150 (MQ=255)
                                                                                          tCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTa                                                             <  2:527455/150‑1 (MQ=255)
                                                                                          tCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTa                                                             <  1:206168/150‑1 (MQ=255)
                                                                                          tCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTa                                                             <  1:1169315/150‑1 (MQ=255)
                                                                                           cGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTaa                                                            >  1:782994/1‑150 (MQ=255)
                                                                                           cGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTaa                                                            >  2:923914/1‑150 (MQ=255)
                                                                                           cGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTaa                                                            >  2:21203/1‑150 (MQ=255)
                                                                                                 aCGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTCACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCa                                                      <  1:715362/150‑1 (MQ=255)
                                                                                                  cGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCAt                                                     <  2:70133/150‑1 (MQ=255)
                                                                                                        gTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCAt                                               >  2:161213/1‑150 (MQ=255)
                                                                                                          ctctcCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTc                                             <  1:1312099/150‑1 (MQ=255)
                                                                                                           tctcCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCg                                            <  2:962156/150‑1 (MQ=255)
                                                                                                           tctcCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCg                                            >  2:1483839/1‑150 (MQ=255)
                                                                                                           tctcCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCg                                            >  1:1330389/1‑150 (MQ=255)
                                                                                                              cctcAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTc                                         <  2:871568/150‑1 (MQ=255)
                                                                                                              cctcAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGACGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTc                                         >  1:848535/1‑150 (MQ=255)
                                                                                                               ctcAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCg                                        <  1:1537453/150‑1 (MQ=255)
                                                                                                                 cAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGAt                                      >  1:1284948/1‑150 (MQ=255)
                                                                                                                 cAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGAt                                      <  1:1556117/150‑1 (MQ=255)
                                                                                                                   aaaCTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATAc                                    <  2:1349034/150‑1 (MQ=255)
                                                                                                                    aaCTGGATTTGCGGAACGATGAACGAAGCAGTCTCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATAcc                                   <  2:453404/150‑1 (MQ=255)
                                                                                                                               cGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGAc                        >  2:1473428/1‑150 (MQ=255)
                                                                                                                                ggAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGAcc                       >  1:762479/1‑150 (MQ=255)
                                                                                                                                    cGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCg                   <  1:1032289/150‑1 (MQ=255)
                                                                                                                                      aTGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGt                 <  1:1596505/150‑1 (MQ=255)
                                                                                                                                         aaCGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTTGGTCGATACCAATACCCAGACCAGCGGTTgg              >  1:747454/1‑150 (MQ=255)
                                                                                                                                         aaCGAAGCAGTCAGGCGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTgg              >  1:331190/1‑150 (MQ=255)
                                                                                                                                          aCGAAGCAGTCAGGCGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTggc             >  1:1383899/1‑150 (MQ=255)
                                                                                                                                          aCGAAGCAGTCAGGCGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTggc             >  1:962703/1‑150 (MQ=255)
                                                                                                                                           cGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTggcg            <  1:1044172/150‑1 (MQ=255)
                                                                                                                                           cGAAGCAGTCAGGCGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTggcg            >  2:824039/1‑150 (MQ=255)
                                                                                                                                            gAAGCAGTCGCCTGACTGCTTCGTTCATTTAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTggcgg           <  2:845120/150‑1 (MQ=255)
                                                                                                                                            gAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTggcgg           >  1:769427/1‑150 (MQ=255)
                                                                                                                                                 aGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTGGCGGCAGAc      <  2:1039793/150‑1 (MQ=255)
                                                                                                                                                 aGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTGGCGGCAGAc      <  2:104435/150‑1 (MQ=255)
                                                                                                                                                 aGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTGGCGGCAGAc      <  2:375161/150‑1 (MQ=255)
                                                                                                                                                    cGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTGGCGGCAGACCat   <  2:184876/150‑1 (MQ=255)
                                                                                                                                                    cGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTGGCGGCAGACCat   <  2:811508/150‑1 (MQ=255)
                                                                                                                                                     gCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTGGCGGCAGACCata  >  1:10317/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTGGCGGCAGACCATA  >  NC_000913/4353021‑4353319

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