Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,070,882 T→A 8.1% E25V (GAG→GTG)  rutE ← putative malonic semialdehyde reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,070,8820TA 8.1% 236.5 / 8.7 82E25V (GAG→GTG) rutEputative malonic semialdehyde reductase
Reads supporting (aligned to +/- strand):  ref base T (39/36);  new base A (4/3);  total (43/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

CCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAAC  >  NC_000913/1070738‑1071028
                                                                                                                                                |                                                                                                                                                  
ccGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTcccc                                                                                                                                                >  1:1590361/1‑148 (MQ=255)
ccGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCg                                                                                                                                                >  2:713302/1‑149 (MQ=255)
  gCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGcc                                                                                                                                              >  2:1178566/1‑149 (MQ=255)
    tGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGGCGGGAGACACCCGTCAg                                                                                                                                            <  1:1490658/149‑1 (MQ=255)
    tGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGGCGGGAGACACCCGTCAg                                                                                                                                            <  1:727692/149‑1 (MQ=255)
       aaaGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCCATTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGCCGGGTGTCACCCGCCAGCCg                                                                                                                                         >  1:103109/1‑149 (MQ=255)
         aGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGtt                                                                                                                                       >  1:321868/1‑149 (MQ=255)
           tGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTtgt                                                                                                                                     <  2:653440/149‑1 (MQ=255)
           tGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTtgt                                                                                                                                     <  2:647929/149‑1 (MQ=255)
            gCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGGCTCATCGCTGGCGGGAGACACCCGTCAGcgatgaga                                                                                                                                    <  1:30666/149‑8 (MQ=255)
             ccGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGa                                                                                                                                   <  1:416888/149‑1 (MQ=255)
                     aGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGc                                                                                                                           >  2:1086334/1‑149 (MQ=255)
                         gTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCg                                                                                                                       >  1:232549/1‑149 (MQ=255)
                         gTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCg                                                                                                                       >  1:698952/1‑149 (MQ=255)
                             ttttCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGa                                                                                                                   >  1:600510/1‑149 (MQ=255)
                              tttCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGaa                                                                                                                  >  2:1488216/1‑149 (MQ=255)
                               ttCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCAATTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAAc                                                                                                                 <  1:1324296/149‑1 (MQ=255)
                                   ttCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAggg                                                                                                             <  1:1161489/149‑1 (MQ=255)
                                   ttCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAggg                                                                                                             <  2:660739/149‑1 (MQ=255)
                                    tCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGt                                                                                                            <  1:124965/149‑1 (MQ=255)
                                            gCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAgcg                                                                                                    >  2:1045143/1‑149 (MQ=255)
                                                   aaCACGATCCGTGCCGGAGAACAGATAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATAGCTGACGGGTGTCTCCCGCCAGCCGTTGGGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTgg                                                                                             >  2:141977/1‑149 (MQ=255)
                                                    acacGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGATCACATTTCATCATGGCATAAATCTCACGTATCGTCTCATCGCTGACGGGTGTCTACCGCCAGCCGTTGTGTGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTggg                                                                                            <  2:125298/149‑1 (MQ=255)
                                                    acacGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGACCCCATTTCATCAGGGCATAAATCTCCAGTAACTTCTAATCGCTGCCGGGTGTCTCCCGCCCGCCGGTGTGAGTGCGGGCCTCGGTGAACAGGGTGCTAAGCGAACCTggg                                                                                            >  2:66635/1‑149 (MQ=255)
                                                     cacGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATATCATCAGGGCATAAATCTCCCGTAAAGTCACATCGCAGACGGGTGTCTCCCGCCAGCCGTTGTGAGAGCGGGCATCGGTGAACAGGGTGCTAAGCGCAGCTGGGc                                                                                           >  1:988069/1‑149 (MQ=255)
                                                                 cGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGAGTCTCCCGCCTGCCGTTGTGAGTGCGGGAATCGGTGCACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGt                                                                               >  2:29829/1‑149 (MQ=255)
                                                                   gagaACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTc                                                                             >  2:554523/1‑149 (MQ=255)
                                                                   gagaACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTc                                                                             <  1:184630/149‑1 (MQ=255)
                                                                    agaACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGGCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCa                                                                            <  2:232549/149‑1 (MQ=255)
                                                                      aaCAGTTAGCTGATGTCGGCCCCAATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGCCGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGCTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATa                                                                          >  1:268300/1‑149 (MQ=255)
                                                                      aaCAGTGAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATa                                                                          <  2:103109/149‑1 (MQ=255)
                                                                        cAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGc                                                                        >  2:885421/1‑149 (MQ=255)
                                                                             aGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTc                                                                   >  2:659014/1‑149 (MQ=255)
                                                                             aGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTc                                                                   >  2:203125/1‑149 (MQ=255)
                                                                              gCTGATGTCGGCCCCCATTTCATCTGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTcc                                                                  <  2:1517427/149‑1 (MQ=255)
                                                                              gCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCCCTGGCGGGAGACACCCGTCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTcc                                                                  >  2:745439/1‑149 (MQ=255)
                                                                              gCTGATGTAGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGACAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTcc                                                                  <  2:1394040/149‑1 (MQ=255)
                                                                               cTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCt                                                                 >  1:493829/1‑149 (MQ=255)
                                                                                tGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCCCCTGGGCTAACGGCTTCGTTCATAGCAATTCCtt                                                                >  1:608141/1‑149 (MQ=255)
                                                                                 gATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTa                                                               <  1:1462684/149‑1 (MQ=255)
                                                                                  aTGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTAc                                                              >  2:1135029/1‑149 (MQ=255)
                                                                                   tGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACa                                                             <  1:752578/149‑1 (MQ=255)
                                                                                        gCCCGCATTTCATCAGGGCATATATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGGTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCTTAGCAATTCCTTACAGggcg                                                        <  1:195526/149‑1 (MQ=255)
                                                                                        gCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGggcg                                                        <  1:1211875/149‑1 (MQ=255)
                                                                                         cccccATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGggcgg                                                       >  1:537294/1‑149 (MQ=255)
                                                                                           cccATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCt                                                     >  2:1240424/1‑149 (MQ=255)
                                                                                               tttCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCAc                                                 >  1:517786/1‑149 (MQ=255)
                                                                                                 tcatcaGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACgg                                               <  1:1178566/149‑1 (MQ=255)
                                                                                                   atcaGGGCAGAAATCTCCCGACACGTCTCATCTCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTATGCGCAACTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTg                                             <  1:66635/149‑1 (MQ=255)
                                                                                                    tgggggCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGa                                            <  2:478298/146‑1 (MQ=255)
                                                                                                    tcaGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGa                                            <  1:1488216/149‑1 (MQ=255)
                                                                                                    tcaGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGa                                            <  1:1086334/149‑1 (MQ=255)
                                                                                                    tcaGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCAGATGTGAGTGAGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCAAGGTGa                                            >  1:916877/1‑149 (MQ=255)
                                                                                                       gggCATAAATCTCCCGTAACGTCTCATCGCTGCCGGGAGGCTCCCGCCCGCCGTTGTGAGTGCGGGCAGCGGTGACCAGGGTGCTAACCGCACGTGGGCTAACGGCTTCGTTCATAGAATATCCTTTCAGGGAGGCTTCACGGtgttgt                                         >  2:270975/1‑149 (MQ=255)
                                                                                                       gggCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGt                                         <  2:6836/149‑1 (MQ=255)
                                                                                                       gggCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGt                                         >  2:455769/1‑149 (MQ=255)
                                                                                                       gggCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGt                                         <  2:698952/149‑1 (MQ=255)
                                                                                                         gCATAAATCTCCCGTAACGGCTCATCGCTGCCGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTaa                                       <  2:23729/149‑1 (MQ=255)
                                                                                                          cATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGGCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAAc                                      <  1:498809/149‑1 (MQ=255)
                                                                                                              aaTCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGc                                  <  1:1135029/149‑1 (MQ=255)
                                                                                                              aaTCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGc                                  >  1:1470367/1‑149 (MQ=255)
                                                                                                                tctcCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTg                                <  2:938950/149‑1 (MQ=255)
                                                                                                                tctcCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGGTGTGAGGGCGGGCATCGGGGAACAGGGCGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTg                                <  1:1401223/149‑1 (MQ=255)
                                                                                                                    ccGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCaa                            >  1:510091/1‑149 (MQ=255)
                                                                                                                       tAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGcc                         >  1:473684/1‑149 (MQ=255)
                                                                                                                          cGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGt                      >  2:354249/1‑149 (MQ=255)
                                                                                                                              tcATCGCTGACGTGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAg                  <  2:321868/149‑1 (MQ=255)
                                                                                                                              tcATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAg                  >  1:1206318/1‑149 (MQ=255)
                                                                                                                              tcATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAg                  >  2:603537/1‑149 (MQ=255)
                                                                                                                               cATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGTGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGt                 >  2:453259/1‑149 (MQ=255)
                                                                                                                               cATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGt                 <  1:713302/149‑1 (MQ=255)
                                                                                                                               cATCGCTGACGGGGGGCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGt                 <  1:37393/149‑1 (MQ=255)
                                                                                                                                aTCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTa                >  2:863947/1‑149 (MQ=255)
                                                                                                                                 tCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTaa               >  1:1437442/1‑149 (MQ=255)
                                                                                                                                 cggCTGGCGGGAGACACCCGTCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGAAAGCCCGTTGAGTaa               >  2:623895/3‑149 (MQ=255)
                                                                                                                                    cTGGCGGGAGACACCCGTCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTGCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACag            >  1:1004465/1‑149 (MQ=255)
                                                                                                                                     tGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACaga           <  2:600510/149‑1 (MQ=255)
                                                                                                                                         gggTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCTCGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCAt       <  2:1590361/149‑1 (MQ=255)
                                                                                                                                            tgtCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCa                                                  <  1:1510234/103‑1 (MQ=255)
                                                                                                                                            tgtCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCa                                                  >  2:1016592/1‑103 (MQ=255)
                                                                                                                                            tgtCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCa                                                  <  1:1016592/103‑1 (MQ=255)
                                                                                                                                            tgtCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCa                                                  >  2:1510234/1‑103 (MQ=255)
                                                                                                                                              tctcCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAAc  >  1:339508/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
CCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAAC  >  NC_000913/1070738‑1071028

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: ‑