Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 573,037 C→T 5.1% A39A (GCC→GCT ninE → DLP12 prophage; conserved protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913573,0370CT 5.1% 254.9 / 2.2 79A39A (GCC→GCTninEDLP12 prophage; conserved protein
Reads supporting (aligned to +/- strand):  ref base C (33/41);  new base T (2/2);  total (35/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 = 7.32e-01

CTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCGTATGTAATGGCAACCCTGAAACGTCTGTACTGGCACATATCC  >  NC_000913/572898‑573183
                                                                                                                                           |                                                                                                                                                  
ctgcGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTgg                                                                                                                                             <  7:299116/147‑1 (MQ=255)
ctgcGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGt                                                                                                                                            >  2:239389/1‑148 (MQ=255)
   cgcgACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTg                                                                                                                                          <  6:2186/147‑1 (MQ=255)
           gCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGGCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAAt                                                                                                                                 <  5:261555/148‑1 (MQ=255)
             aggaggCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATgg                                                                                                                               <  1:239389/148‑1 (MQ=255)
                  gCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCg                                                                                                                          <  7:30816/148‑1 (MQ=255)
                  gCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGGCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCg                                                                                                                          <  8:299622/148‑1 (MQ=255)
                   cGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGa                                                                                                                         <  7:143974/148‑1 (MQ=255)
                     tAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACt                                                                                                                       >  1:115009/1‑148 (MQ=255)
                     tAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACt                                                                                                                       >  2:107812/1‑148 (MQ=255)
                      aaTGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTg                                                                                                                      >  8:155274/1‑148 (MQ=255)
                               acCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAagg                                                                                                             <  5:46117/148‑1 (MQ=255)
                                ccGCTTATTCGTGGCATGAACGAACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAagga                                                                                                            <  1:148980/148‑1 (MQ=255)
                                 cGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACt                                                                                                                       >  7:246969/1‑136 (MQ=255)
                                 cGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACt                                                                                                                       <  8:246969/136‑1 (MQ=255)
                                       ttCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCa                                                                                                     <  1:107812/148‑1 (MQ=255)
                                            gtCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGgc                                                                                                                    <  8:61362/128‑1 (MQ=255)
                                            gtCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGgc                                                                                                                    >  7:61362/1‑128 (MQ=255)
                                                tGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGAtt                                                                                            >  1:196820/1‑148 (MQ=255)
                                                   aCGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGAttt                                                                                           >  8:409011/1‑146 (MQ=255)
                                                   aCGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTg                                                                                          <  7:409011/147‑1 (MQ=255)
                                                        cacaTCTACCGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGaaaa                                                                                    <  2:210915/148‑1 (MQ=255)
                                                           aTCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGAAGCAAGGAGGAGTCATGGCTGATTTGAGAAAagca                                                                                 <  7:155274/148‑1 (MQ=255)
                                                               aCAGAGTACCAAATCGTCGCAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTg                                                                                               >  3:161198/1‑130 (MQ=255)
                                                               aCAGAGTACCAAATCGTCGCAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTg                                                                                               <  4:161198/130‑1 (MQ=255)
                                                                cAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGaaa                                                                                     <  7:39821/139‑1 (MQ=255)
                                                                cAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGaaa                                                                                     >  8:39821/1‑139 (MQ=255)
                                                                 agagTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGt                                                                           <  2:196820/148‑1 (MQ=255)
                                                                   agTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGgcagca                                                                                                                <  8:251125/109‑1 (MQ=255)
                                                                   agTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGgcagca                                                                                                                >  7:251125/1‑109 (MQ=255)
                                                                    gTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGt                                                                        >  5:309903/1‑148 (MQ=255)
                                                                       ccAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCg                                                                      <  3:39786/147‑1 (MQ=255)
                                                                        cAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCggg                                                                    <  4:166582/148‑1 (MQ=255)
                                                                        cAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCgg                                                                     >  3:201426/1‑147 (MQ=255)
                                                                         aaaTCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGa                                                                   <  6:309903/148‑1 (MQ=255)
                                                                            tcgtcgTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATg                                                                >  1:281079/1‑148 (MQ=255)
                                                                            tcgtcgTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATg                                                                >  8:161412/1‑148 (MQ=255)
                                                                               tcgtAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCa                                                             <  7:195761/148‑1 (MQ=255)
                                                                                  tAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAgg                                                           >  4:90700/1‑147 (MQ=255)
                                                                                    aGCGTAAACGTGAGCTGAAGCCATCCGAAATACCAACACTGCTCTGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGTCTGGGAGCAAGGAGGAGTCATGGCTGATTTGAGATAAGCAGCGCGTTGTCGGGAATGCCAGGTa                                                         <  3:304083/147‑1 (MQ=255)
                                                                                    aGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTa                                                         <  3:259454/147‑1 (MQ=255)
                                                                                    aGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGAT‑‑AAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGa                                                      <  7:28527/148‑1 (MQ=255)
                                                                                     gCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAg                                                       <  4:201426/148‑1 (MQ=255)
                                                                                         aaaCCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATc                                                   >  4:12226/1‑148 (MQ=255)
                                                                                         aaaCCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATc                                                   >  7:21769/1‑148 (MQ=255)
                                                                                          aaCCTGAGCTGAAGCCATCCGAAGTACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATc                                                   <  5:98062/147‑1 (MQ=255)
                                                                                            ccTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCt                                                <  3:90700/148‑1 (MQ=255)
                                                                                             cTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTg                                               >  1:285974/1‑148 (MQ=255)
                                                                                               gAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGc                                             >  4:191826/1‑148 (MQ=255)
                                                                                                 gCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCGt                                           <  2:281079/148‑1 (MQ=255)
                                                                                                      aGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCGTATGTa                                      <  8:21769/148‑1 (MQ=255)
                                                                                                      aGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGCGAATGCCAGGTAAGAATCCCTGGCGTa                                          >  8:94579/1‑144 (MQ=255)
                                                                                                      aGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGCGAATGCCAGGTAAGAATCCCTGGCGTa                                          <  7:94579/144‑1 (MQ=255)
                                                                                                       gCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCGTATGTaa                                     >  4:368485/1‑148 (MQ=255)
                                                                                                       gCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCGTATGTaa                                     <  2:115009/148‑1 (MQ=255)
                                                                                                          aTCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCGTATGTAATgg                                  >  1:246318/1‑148 (MQ=255)
                                                                                                            ccGAAATACCAACACTGCTCGGATATACCGCTAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCGTATGTAATg                                   <  1:121860/145‑1 (MQ=255)
                                                                                                            ccGAAATACCAACACTGCTCGGATATACCGCTAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCGTATGTAATg                                   <  1:127259/145‑1 (MQ=255)
                                                                                                            ccGAAATACCAACACTGCTCGGATATACCGCTAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCGTATGTAATg                                   >  2:121860/1‑145 (MQ=255)
                                                                                                            ccGAAATACCAACACTGCTCGGATATACCGCTAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCGTATGTAATg                                   >  2:127259/1‑145 (MQ=255)
                                                                                                            ccGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCa                                                             <  7:97544/119‑1 (MQ=255)
                                                                                                            ccGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCa                                                             >  8:97544/1‑119 (MQ=255)
                                                                                                            ccGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTa                                                         <  6:312963/123‑1 (MQ=255)
                                                                                                            ccGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTa                                                         >  5:312963/1‑123 (MQ=255)
                                                                                                               aaaTCCCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCGGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCGTATGTAATGGCAAcc                             <  5:91094/148‑1 (MQ=255)
                                                                                                               aaaTACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCGTATGTAATGGCAAcc                             <  3:191826/148‑1 (MQ=255)
                                                                                                               aaaTACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCGTATGTAATGGCAAcc                             >  1:156644/1‑148 (MQ=255)
                                                                                                                            gCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGa                                                                   <  4:39921/97‑1 (MQ=255)
                                                                                                                            gCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGa                                                                   >  3:39921/1‑97 (MQ=255)
                                                                                                                             cTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCGTATGTAATGGCAACCCTGAAACGTCTGTa               <  2:156644/148‑1 (MQ=255)
                                                                                                                               cGGATATACATCAAGCTTGGTTTATAAAAAATGGTTGCGACTTGAAGCAAGGATGAGTCATGGCTCATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGCCGTATGTAATGGCa                                <  1:302077/129‑1 (MQ=255)
                                                                                                                                ggATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAGTGCCAGGTAAGAATCCCTGGCGTATGTAATGGCAACCCTGAAACGTCTGTACTg            <  8:320737/148‑1 (MQ=255)
                                                                                                                                ggATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCGTATGTAATGGCAACCCTGAAACGTCTGTACTg            >  1:66732/1‑148 (MQ=255)
                                                                                                                                    ataCCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGgcagca                                                                                                                <  3:188266/44‑1 (MQ=255)
                                                                                                                                    ataCCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGgcagca                                                                                                                >  4:188266/1‑44 (MQ=255)
                                                                                                                                    ataCCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGTAGGAGTCAGGGATGCTTAGAGCAAAGCAGCGCGTCGTCGGGAATGCCAGGTAAGACTCCCTGGCGTATGTAATGGCAACCCTGAAACGTCTGGACTGGca         >  7:289150/1‑147 (MQ=255)
                                                                                                                                    ataCCGCCAGCTTGGTCGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCGTATGTAATGGCAACCCTGAAACGTCTGTACTGGca         >  5:305392/1‑147 (MQ=255)
                                                                                                                                       ccGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCGTATGTAATGGCAACCCTGAAACGTCTGTACTGGCACata     >  6:227389/1‑148 (MQ=255)
                                                                                                                                       ccGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCGTATGTAATGGCAACCCTGAAACGTCTGTACTGGCACata     >  6:167192/1‑148 (MQ=255)
                                                                                                                                          ccAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCGTATGTAATGGCAACCCTGAAACGTCTGTACTGGCACATATcc  >  1:337070/1‑148 (MQ=255)
                                                                                                                                           |                                                                                                                                                  
CTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCGTATGTAATGGCAACCCTGAAACGTCTGTACTGGCACATATCC  >  NC_000913/572898‑573183

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