Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,070,887 C→T 7.2% W23* (TGG→TGA rutE ← putative malonic semialdehyde reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,070,8870CT 7.2% 252.8 / 5.8 83W23* (TGG→TGArutEputative malonic semialdehyde reductase
Reads supporting (aligned to +/- strand):  ref base C (40/36);  new base T (3/3);  total (43/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTTTCG  >  NC_000913/1070740‑1071034
                                                                                                                                                   |                                                                                                                                                   
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:647929/149‑1 (MQ=255)
         tGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTtgt                                                                                                                                           <  2:653440/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:698952/149‑1 (MQ=255)
                                                                                                     gggCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGt                                               >  2:455769/1‑149 (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                                                        >  2:1510234/1‑103 (MQ=255)
                                                                                                                                          tgtCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCa                                                        <  1:1016592/103‑1 (MQ=255)
                                                                                                                                            tctcCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAAc        >  1:339508/1‑149 (MQ=255)
                                                                                                                                                ccGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGTGTGCTAAGCGCACCGGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGttt    <  2:1047424/149‑1 (MQ=255)
                                                                                                                                                ccGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGttt    >  2:559980/1‑149 (MQ=255)
                                                                                                                                                  gCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTTTCg  >  1:1178135/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTTTCG  >  NC_000913/1070740‑1071034

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: