Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,070,873 A→G 7.7% V28A (GTC→GCC)  rutE ← putative malonic semialdehyde reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,070,8730AG 7.7% 217.5 / 8.2 80V28A (GTC→GCC) rutEputative malonic semialdehyde reductase
Reads supporting (aligned to +/- strand):  ref base A (36/33);  new base G (3/3);  total (43/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.52e-01

GCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGA  >  NC_000913/1070730‑1071019
                                                                                                                                               |                                                                                                                                                  
gCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGgt                                                                                                                                               >  1:6836/1‑149 (MQ=255)
gCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGgt                                                                                                                                               <  1:1184520/149‑1 (MQ=255)
 cAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGgtg                                                                                                                                              >  2:184630/1‑149 (MQ=255)
        ccGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTcccc                                                                                                                                       >  1:1590361/1‑148 (MQ=255)
        ccGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCg                                                                                                                                       >  2:713302/1‑149 (MQ=255)
          gCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGcc                                                                                                                                     >  2:1178566/1‑149 (MQ=255)
            tGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGGCGGGAGACACCCGTCAg                                                                                                                                   <  1:727692/149‑1 (MQ=255)
            tGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGGCGGGAGACACCCGTCAg                                                                                                                                   <  1:1490658/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:698952/1‑149 (MQ=255)
                                 gTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCg                                                                                                              >  1:232549/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                                                                                                    <  2:660739/149‑1 (MQ=255)
                                           ttCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAggg                                                                                                    <  1:1161489/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:1086334/149‑1 (MQ=255)
                                                                                                            tcaGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGa                                   <  1:1488216/149‑1 (MQ=255)
                                                                                                            tcaGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCAGATGTGAGTGAGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCAAGGTGa                                   >  1:916877/1‑149 (MQ=255)
                                                                                                               gggCATAAATCTCCCGTAACGTCTCATCGCTGCCGGGAGGCTCCCGCCCGCCGTTGTGAGTGCGGGCAGCGGTGACCAGGGTGCTAACCGCACGTGGGCTAACGGCTTCGTTCATAGAATATCCTTTCAGGGAGGCTTCACGGtgttgt                                >  2:270975/1‑149 (MQ=255)
                                                                                                               gggCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGt                                <  2:698952/149‑1 (MQ=255)
                                                                                                               gggCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGt                                <  2:6836/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         >  2:603537/1‑149 (MQ=255)
                                                                                                                                      tcATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAg         >  1:1206318/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)
                                                                                                                                               |                                                                                                                                                  
GCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCATCAGGGCATAAATCTCCCGTAACGTCTCATCGCTGACGGGTGTCTCCCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGA  >  NC_000913/1070730‑1071019

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: