Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 578,831 T→C 100% T21T (ACA→ACG)  borD ← DLP12 prophage; putative lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913578,8310TC98.7% 10622.3 / 31.2 3704T21T (ACA→ACG) borDDLP12 prophage; putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (27/24);  new base C (1849/1789);  total (1882/1822)
Fisher's exact test for biased strand distribution p-value = 7.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.59e-03
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 3705 total aligned reads displayed.

TGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTC  >  NC_000913/578685‑578977
                                                                                                                                                  |                                                                                                                                                  
 gggTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTc                                                                                                                                                 >  5:154607/1‑149 (MQ=255)
  ggTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCt                                                                                                                                                <  7:122029/149‑1 (MQ=255)
      tCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGGGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGtt                                                                                                                                             <  1:100289/148‑1 (MQ=255)
        tGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAg                                                                                                                                          >  8:315237/1‑149 (MQ=255)
          ttttAACAACATTTTCTGCGCCGCCTCAATTTTTAGCTGCATCGAGAGTTTTCTTCTGCCCAATTCCAGATACTAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGGTTGTTTTGAACAGTAAACGTCTGTTGa                                                                                                                                           <  3:204291/146‑1 (MQ=255)
          ttttAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGc                                                                                                                                         <  5:207814/148‑1 (MQ=255)
          ttttAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGca                                                                                                                                        >  1:335573/1‑149 (MQ=255)
          ttttAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGc                                                                                                                                         <  5:56553/148‑1 (MQ=255)
              aacaacATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATc                                                                                                                                    >  6:280865/1‑149 (MQ=255)
              aacaacATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATAATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATc                                                                                                                                    >  8:226840/1‑149 (MQ=255)
                caacaTTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCt                                                                                                                                  <  1:124142/149‑1 (MQ=255)
                caacaTTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCt                                                                                                                                  <  8:125276/149‑1 (MQ=255)
                caacaTTATCTGCGCCGCCACAAATTTTAGCTGCATCCACAGTTTTCTTCTGCCCAATTCCATAAAAGAAGAAATGATGGCTGATGGTTTCCTATGCTGCTACTGAAGCCGGTATGTTTTGAACAGTAAAAGACTGTTGAGCACATc                                                                                                                                    >  4:111193/1‑147 (MQ=255)
                 aacaTTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTg                                                                                                                                 >  3:34083/1‑149 (MQ=255)
                 aacaTTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTg                                                                                                                                 >  7:292348/1‑149 (MQ=255)
                 aacaTATTCTGAGCCGCCACAAAATTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCACAAACGAAGAAAAGATGGGTCATGCTTTCCTTTGCTCCTACTGCTCCCGGTTTCTTTTGAACAGAAAACGTCTGTTCAGCACATCCt                                                                                                                                  >  2:178236/1‑148 (MQ=255)
                   caTTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGt                                                                                                                                >  7:96694/1‑148 (MQ=255)
                    aTTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGtaa                                                                                                                              <  7:3680/149‑1 (MQ=255)
                       ttCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTTGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGtaata                                                                                                                            <  1:200964/148‑1 (MQ=255)
                       ttCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGtaata                                                                                                                            <  6:95457/148‑1 (MQ=255)
                            cgccgcCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAgg                                                                                                                      <  1:312279/149‑1 (MQ=255)
                            cgccgcCACAAATTTTAGCGGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAgg                                                                                                                      <  4:341679/149‑1 (MQ=255)
                              ccgccACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAggg                                                                                                                     <  3:312966/148‑1 (MQ=255)
                                 ccACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGCAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAg                                                                                                                 >  6:11125/1‑149 (MQ=255)
                                   acaAATTTTAGCTGCATCGACAGTTTTCTTCTGCGCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAgcg                                                                                                               >  8:179713/1‑149 (MQ=255)
                                     aaaTTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAgcg                                                                                                               >  2:345297/1‑147 (MQ=255)
                                     aaaTTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAgcg                                                                                                               >  6:250900/1‑147 (MQ=255)
                                             gCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAGGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAg                                                                                                     <  7:174124/149‑1 (MQ=255)
                                             gCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACTTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGa                                                                                                      <  4:225106/148‑1 (MQ=255)
                                             gCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGc                                                                                                        >  1:14826/1‑146 (MQ=255)
                                             gCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGa                                                                                                      <  2:28845/148‑1 (MQ=255)
                                             gCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGa                                                                                                      <  8:291027/148‑1 (MQ=255)
                                                gCATCGACAGTTATCTTATGCCCAATTCCATAAAATAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTTCTGCAGGTTTGTTATGAACAGTAAAAGTCTGTTGAGCACATACTGTAATAAGCATGGCCAGAGAATTATCGAGTAg                                                                                                  >  5:348143/1‑149 (MQ=255)
                                                 cATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTa                                                                                                   <  6:138743/147‑1 (MQ=255)
                                                      aCAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCAtttt                                                                                            <  6:349172/149‑1 (MQ=255)
                                                       cAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCa                                                                                                >  1:74576/1‑144 (MQ=255)
                                                        aGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCAtttttt                                                                                          <  8:34884/149‑1 (MQ=255)
                                                         gTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCAttttttt                                                                                         >  1:175045/1‑149 (MQ=255)
                                                         gTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCAttttttt                                                                                         <  8:42761/149‑1 (MQ=255)
                                                           tttCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCa                                                                                       >  4:28118/1‑149 (MQ=255)
                                                            ttcttcTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCa                                                                                       <  2:374890/148‑1 (MQ=255)
                                                              cttctGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATg                                                                                     <  3:254080/148‑1 (MQ=255)
                                                               ttctGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAgg                                                                                                                      <  6:219900/114‑1 (MQ=255)
                                                                 ctGCACAATTCCTGAAACGAAGAAATGCTGGGTGCTGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGTCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGgtg                                                                                  <  3:85087/148‑1 (MQ=255)
                                                                    cccAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGc                                                                                                                    <  8:255463/111‑1 (MQ=255)
                                                                    cccAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTAt                                                                              >  3:183556/1‑149 (MQ=255)
                                                                     ccAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTAtt                                                                             >  6:19058/1‑149 (MQ=255)
                                                                     ccAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTAtt                                                                             >  4:65539/1‑149 (MQ=255)
                                                                       aaTTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGgtgt                                                                                 <  3:15231/143‑1 (MQ=255)
                                                                       aaTTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTcc                                                                           >  1:28596/1‑149 (MQ=255)
                                                                            cAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATg                                                                      >  7:19993/1‑149 (MQ=255)
                                                                            cAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATg                                                                      >  5:37771/1‑149 (MQ=255)
                                                                                  cGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGCTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCttttt                                                                <  5:289715/149‑1 (MQ=255)
                                                                                   gaagaaATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCttttt                                                                >  2:320313/1‑148 (MQ=255)
                                                                                   gaagaaATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTg                                                               <  5:188627/149‑1 (MQ=255)
                                                                                    aagaaATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGTTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTCGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTg                                                               >  5:84147/1‑148 (MQ=255)
                                                                                     agaaATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGaa                                                             <  4:310129/149‑1 (MQ=255)
                                                                                     agaaATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGaa                                                             <  8:376037/149‑1 (MQ=255)
                                                                                      gaaATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTAATGGTGTTATTCCCGATGCTTTTTGAAg                                                            >  1:40969/1‑149 (MQ=255)
                                                                                                gTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGAGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAAt                                                  >  6:314044/1‑149 (MQ=255)
                                                                                                 tGATGGTTTCCTTTTGTTCTACTGCTGCCTGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCTGGGCCAGCGCAGTAGCGATTAGCATTTTTTTCATGGTGTTATTCCCGATTCTTTTTGAAGTTCGCTGAATc                                                 <  4:198536/149‑1 (MQ=255)
                                                                                                 tGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATc                                                                                                                                    <  5:262109/66‑1 (MQ=255)
                                                                                                 tGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATc                                                 >  4:143486/1‑149 (MQ=255)
                                                                                                  gATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCg                                                >  3:127470/1‑149 (MQ=255)
                                                                                                    tGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTa                                              <  1:358961/149‑1 (MQ=255)
                                                                                                      gTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGaa                                                             <  4:97835/132‑1 (MQ=255)
                                                                                                      gTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTAt                                             <  6:37771/148‑1 (MQ=255)
                                                                                                      gTTTCCTTTGGTGCTACTGCGGCGGGTTTGTTTTGAACAGTATACGTCTTTTGAGCACATCCTGTAATAAGCAGTGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTGCCGATTCTTTTTGAAGTTCGCAGAATCGTAt                                             <  1:228516/148‑1 (MQ=255)
                                                                                                       tttCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTAtgt                                           <  7:393658/149‑1 (MQ=255)
                                                                                                        ttCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTAtgtg                                          >  4:366411/1‑149 (MQ=255)
                                                                                                            tttGGTGCTACTGCTGCCGGTTTGATTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAg                                       >  2:47033/1‑148 (MQ=255)
                                                                                                                gTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGaaaa                                   >  2:261477/1‑148 (MQ=255)
                                                                                                                gTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGaaaa                                   <  7:68630/148‑1 (MQ=255)
                                                                                                                 tGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAAtt                                 <  8:202312/149‑1 (MQ=255)
                                                                                                                  gCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAAtt                                 >  3:558/1‑148 (MQ=255)
                                                                                                                   cTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTa                                <  2:237209/148‑1 (MQ=255)
                                                                                                                       tgctgcCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTaaaca                            >  2:70067/1‑148 (MQ=255)
                                                                                                                       tgctgcCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTaaaca                            >  7:345046/1‑148 (MQ=255)
                                                                                                                       tgctgcCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTaaaca                            >  7:318144/1‑148 (MQ=255)
                                                                                                                          tgccggtttgtttTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAAcc                        >  4:170501/6‑149 (MQ=255)
                                                                                                                            ccggtttgtttTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAAccc                       >  5:233145/4‑148 (MQ=255)
                                                                                                                             cggtttgtttTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCttt                                                                  >  8:73059/3‑104 (MQ=255)
                                                                                                                               gtttgtttTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTaaa                   <  7:259739/149‑1 (MQ=255)
                                                                                                                               gttggttgTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGGGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAACCCCTaaa                   <  2:215101/149‑1 (MQ=255)
                                                                                                                                tttgtttTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTaaa                   >  1:258767/1‑148 (MQ=255)
                                                                                                                                tttgtttTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTaaa                   >  3:361292/1‑148 (MQ=255)
                                                                                                                                gttgtttTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCATGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAGCAAACCCTaaa                   >  6:16300/2‑148 (MQ=255)
                                                                                                                                 ttgtttTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTTGAAAATTAAACAAACCCTAAAc                  >  2:122634/1‑148 (MQ=255)
                                                                                                                                   gtttTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGGAGAAAATTAAACAAACCCTAAACAAt               >  3:285711/1‑149 (MQ=255)
                                                                                                                                   gtttTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCAGTTTTTTCATGGTGTTATTGCCGATGCTTTTTGAAGTTCGCAGACTCGTATGTGTAGAAAATTAAACAAACCCTAAACAAt               >  6:63426/1‑149 (MQ=255)
                                                                                                                                    ttttGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAAt               >  7:146683/1‑148 (MQ=255)
                                                                                                                                    ttttGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAAt               >  4:251486/1‑148 (MQ=255)
                                                                                                                                      ttGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGa             <  3:234226/148‑1 (MQ=255)
                                                                                                                                        gAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGtt          >  2:151366/1‑149 (MQ=255)
                                                                                                                                         aaCAGTAAGCGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCATTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGtt          >  3:103059/1‑148 (MQ=255)
                                                                                                                                         aaCAGTAAACGTCTGTTGAGCACATCCTGTATTAAGCAGGGACAGCGCAGTAGCGAGTAGCAATTTTTTCAAGGTGTGATTCACGAAGCTTTTTGAAGTACGCAGAATAGTATGTGTACAAAATAAAACAAACACTAAACAATTAGtt          >  2:100082/1‑148 (MQ=255)
                                                                                                                                         aaCAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGtt          >  5:312996/1‑148 (MQ=255)
                                                                                                                                          aCAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTg         >  7:234805/1‑148 (MQ=255)
                                                                                                                                             gTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGATCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAAt     <  8:149677/149‑1 (MQ=255)
                                                                                                                                              tAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGt           <  5:12619/142‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTC  >  NC_000913/578685‑578977

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