Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,491,466 G→A 5.3% A354V (GCC→GTC)  idnT ← L‑idonate and D‑gluconate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,491,4660GA 5.3% 216.4 / 2.7 78A354V (GCC→GTC) idnTL‑idonate and D‑gluconate transporter
Reads supporting (aligned to +/- strand):  ref base G (40/32);  new base A (2/2);  total (43/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.61e-01

CCTTTAAATAGCCAGAAGCCAGGGTCGTTTACGTGTGACGCGATCACGCTGCCCGCACCGGTTGCCAGTACCATTAAAGCGGGATCGGCATGGGTAACGTTGATAATCGGCAACACCACACCCGCGGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCAGCACATCAATAACGGCGAAAGTGTAGTTCCGGTCATTAAGTGTGAAATATAGTGCCCGACACCGCTATCTACTAATACCTGCTTAAACGCGCCGCCGCCAGCGAT  >  NC_000913/4491319‑4491611
                                                                                                                                                   |                                                                                                                                                 
ccTTTCAATACCCCGACGCCAGGGTCGGTTACATGAGACGAGACCACGCTGCCCGCACCGGGTGCCAATACCATTAAAGCGGGATCGGCATGGGTAACGTTGATAACCGGCAACACCACCCCCGCGGAGGTAATGGCCGCGACGGTGGc                                                                                                                                                  <  2:968541/149‑1 (MQ=255)
 cTTTAAATAGCCAGAAGCCAGGGTCGTTTACGTGTGACGCTATCACGCTGCCCGCACCGGTTGCCAGTACCATTAAAGCGGGATCGGCATGGGTAACGTTGATAATCGGCAACACCACACCCGCGGTGGTAATGGCCGCGCCGGTGGcc                                                                                                                                                 >  2:967262/1‑148 (MQ=255)
 cTTTAAATAGCCAGAAGCCAGGGTCGTTTACGTGTGACGCGATCACGCTGCCCGCACCGGTTGCCAGTACCATTAAAGCGGGATCGGCATGGGTAACGTTGATAATCGGCAACACCACACCCGCGGCGGTAATGGCCGCGCCGGTGGCa                                                                                                                                                 >  2:327327/1‑149 (MQ=255)
  tttAAATAGCCAGAAGCCAGGGTCGTTTACGTGTGACGCGATCACGCTGCCCGCACCGGTTGCCAGTACCATTAAAGCGGGATCGGCATGGGTAACGTTGATAATCGGCAACACCACACCCGCGGTGGTAATGGCCGCGACGGTGGcgg                                                                                                                                                >  2:1268554/1‑147 (MQ=255)
   ttAAATAGCCAGAAGCCAGGGTCGTTTACGTGTGACGCGATCACGCTGCCCGCACCGGTTGCCAGTACCATTAAAGCGGGATCGGCATGGGTAACGTTGATAATCGGCAACACCACACCCGCGGTGGTAATGGCCGCGACGGTGACaga                                                                                                                                               >  2:1138599/1‑149 (MQ=255)
   ttAAATAGCCAGAAGCCAGGGTCGTTTACGTGTGACGCGATCACGCTGCCCGCACCGGTTGCCAGTACCATTAAAGCGGGATCGGCATGGGTAACGTTGATAATCGGCAACACCACACCCGCGGTGGTAATGGCCGCGACGGGGACaga                                                                                                                                               >  2:1399577/1‑149 (MQ=255)
       aTAGCCAGAAGCCAGGGTCGTTTACGTGTGACGCGATCACTCTGCCCGCACCGGTTGCCAGTACCATTAAAGCGGGATCGGCATGGGAAACGTTGATAATCGGCAACACCACACCCGCGCTGGTAACGGCCGCGACGGGGCCAGAGccc                                                                                                                                           >  2:1086220/1‑149 (MQ=255)
         aGCCAGAAGCCAGGGTCGTTTACGTGTGACGCGATCACGCTGCCCGCACCGGTTGCCAGTACCATTAAAGCGGGATCGGCATGGGTAACGTTGATAATCGGCAACACCACACCCGCGGTGGTAATGGCCGCGACGGTGGgagacacaag                                                                                                                                         >  2:1187851/1‑143 (MQ=255)
           ccAGAAGCCAGGGTCGTTTACGTGTGACGCGATCACGCTGCCCGCACCGGTTGCCAGTACCATTAAAGCGGGATCGGCATGGGTAACGTTGATAATCGGCAACACCACACCCGCGGTGGTAATGGCCGCGACGGTGGCAGAGCCCagag                                                                                                                                       <  1:1206205/149‑1 (MQ=255)
             aGAAGCCAGGGTCGTTTACGTGTGACGCGAGCACGCTGCCCGCACCGGTGGCCAGTACCATTAAAGCGGGATCGGCATGGGTAACGTTGATAATCGGCAGCACCACACCCGCGGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCg                                                                                                                                     <  2:341637/149‑1 (MQ=255)
                  ccAGGGTCGTTTACGTGTCACGCGATCACGCTGCCCGCACCGGTTGCCAGTACCATAAAAGCGGGATCGGCATGGGTAACGTTGATAATCGGCAACACCACACCCGCGGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACg                                                                                                                                >  2:937782/1‑149 (MQ=255)
                    aGGGTCGTTTACGTGTGACGCGATCACGCTGCCCGCACCGGTTGCCAGTACCATTAAAGCGGGATCGGCATGGGTAACGTTGATAATCGGCAACACCACACCCGCGGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCa                                                                                                                              >  2:1022261/1‑149 (MQ=255)
                        tCGTTTACGTGTGACGCGATCACGCTGCCCGCACCGGTTGCCAGTACCATTAAAGCGGGGTCGGCATGGGTAACGTTGATAATCGGCAACACCACACCCGCGGTGGTAATGGGCGGGGGGGTGGCAGAGCCCAGGGGGATACGCAagag                                                                                                                          >  1:314515/1‑149 (MQ=255)
                           tttACGTGTGACGCGATCACGCTGCCCGCACCGGTTGCCAGTACCATTAAAGCGGGATCGGCATGGGTAACGTTGATAATCGGCAACACCACACCCGCGGTGGTAATGGCCGCGACGGTGGCAGGGCCCAAAGCGATACGCAACcgcga                                                                                                                       >  2:307749/1‑148 (MQ=255)
                           tttACGTGTGACGCGATCACGCTGCCCGCACCGGTTGCCAGTACCATTAAAGCGGGATCGGCATGGGTAACGTTGATAATCGGCAACACCACACCCGCGGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAgcgc                                                                                                                       >  1:824323/1‑149 (MQ=255)
                             tACTTGTGACGCTATCACGCTGCCCGCCCCTGTTGCCAGTACCATTAAAGCGGGATCGGCATGGGTAACGTTGATAATCGGCAACACCACACCCGCGGTGGTAATGGCCGCGACGGGGTCAGATCCCAGAGCGATACGCAACAGcgccg                                                                                                                     <  1:1074724/149‑1 (MQ=255)
                             tACGTGTGACGCGATCACGCTGCCCGCACCGGTTGCCAGTACCATTAAAGCGGGATCGGCATGGGTAACGTTGATAATCGGCAACACCACACCCGCGGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGcgccg                                                                                                                     <  1:85447/149‑1 (MQ=255)
                                gtgtGACGCGATCACGCTGCCCGCACCGGTTGCCAGTACCATTAAAGCGGGATCGGCATGGGTAACGTTGATAATCGGCAACACCACACCCGCGGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCac                                                                                                                  <  1:139491/149‑2 (MQ=255)
                                    gACCCGATCACGCTGCCCGCACCGGTTGCCAGGACCATTAAAGCGGGATCGGCATGGGTAACGTTGATAATCGGCAACACCACACCCGCGGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGt                                                                                                              <  2:559779/149‑1 (MQ=255)
                                       gcgATCACGCTGCCCGCACCGGTTGCCAGTACCATTAAAGCGGGATCGGCATGGGTAACGTTGATAATCGGCAACACCACACCCGCGGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCa                                                                                                           >  1:3597/1‑149 (MQ=255)
                                       ccgAACACCCTGCCCATCCCGGTTGCCAGGACCAGTAAAGCGGAATCGGTATGGGTACCTTTGATAATCGTCAACTCCACACTCGCGGTGGTAATGGCCGCGACGTTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCa                                                                                                           <  2:1129742/148‑1 (MQ=255)
                                        ccaTCACGCTGCCCGTCCCGGTGGCCAGTACCATTAAAGCGGGATCGGCATGGGTAACGTTGATAAGCGGCAACACCACACCCGCGGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCAg                                                                                                          <  1:378838/147‑1 (MQ=255)
                                          aTCACGCTGCCCGCACCGGTTGCCAGTACCATTAAAGCGGGATCGGCATGGGTAACGTTGATAATCGGCAACACCACACCCGCGGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCTCAACAGTCCAGca                                                                                                        <  1:993630/149‑1 (MQ=255)
                                           tCACGCTGCCCGCACCGGTTGCCAGTACCATTAAAGCGGGATCGGCATGGGTAACGTTGATAATCGGCAACACCACACCCGCGGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAGCAGCGCCGCAACAGTCCAGcac                                                                                                       >  1:54334/1‑149 (MQ=255)
                                                    ccGCCCAGGTTGCCAGTACCATTAAAGCGGGATCGGCATGGGTAACGTTGATAATCGGCAACACCACACCCGCGGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGGCCAGCACATCAATGAc                                                                                              >  2:1335875/1‑149 (MQ=255)
                                                       cACCGGTTGCCAGTACCATTAAAGCGGGATCGGCATGGGTAACGTTGATAATCGGCAACACCACACCCGCGGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCAGCACATCAATAACGGc                                                                                           <  1:238327/149‑1 (MQ=255)
                                                          cGGTTGCCAGTACCATTAAAGCGGGATCGGCATGGGTAACGTTGATAATCGGCAACACCACACCCGCGGTGGTACTGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCAGCACACCAATAACGGCGaa                                                                                        >  2:334994/1‑149 (MQ=255)
                                                          cGGTTGCCAGTACCATTAAAGCGGAATCGGCATGGGTGACGTTGATAATCGGCAACACCACCCCCGCCGTGGTAATGGCCGCGACGGTGGCGTAGCCCAGAGCCGTACGCACCACCGCAGCAACAGTCCCACACATCAAAACCGGGGaa                                                                                        >  2:49533/1‑149 (MQ=255)
                                                            gTTGCCAGTACCATTAAAGCGGGATCGGCATGGGTAACGTTGATAATCGGCAACACCACACCCGCGGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCAGCACATCAATAACGGCGAAAg                                                                                      >  1:1261973/1‑149 (MQ=255)
                                                             ttGCCAGTACCATTAAAGCGGGATCGGCATGGGTAACGTTGATAATCGGCAACACCACACCCGCGGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGGCCAGCACATCAATCACGGCGAAAga                                                                                     >  1:591340/1‑148 (MQ=255)
                                                                ccAGTACCATTAAAGCGGGATCGGCATGGGTAACGTTGATAATCGGCAACACCACACCCGCGGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACCGTCCAGCACATCAATAACGGCGAAAGTGTa                                                                                  <  2:97180/149‑1 (MQ=255)
                                                                ccAGTACCATTAAAGCGGGATCGGCATGGGTAACGTTGATAATCGGCAACACCACACCCGCGGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCAGCACATCAATAACGGCGAAAGTGTa                                                                                  <  1:838710/149‑1 (MQ=255)
                                                                ccAGTACCATTAAAGCGGGATCGGCATGGGTAACGTTGATAATCGGCAACACCACACCCGCGGTGGTAATGGCCGCGACGGTGACAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCAGCACATCAATAACGGCGAAAGTGTa                                                                                  <  1:1399577/149‑1 (MQ=255)
                                                                ccAGTACCATTAAAGCGGGATCGGCATGGGTAACGTTGATAATCGGCAACACCACACCCGCGGTGGTAATGGCCGCGACGGTGACAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCAGCACATCAATAACGGCGAAAGTGTa                                                                                  <  1:1138599/149‑1 (MQ=255)
                                                                   gTACCATTAAAGCGGGAGCGGCATGGGTAACGTTGATAATCGGCAACACCACACCCGCGGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCAGCACATCAATAACGGCGAAAGTGTAGtt                                                                               <  1:1311769/149‑1 (MQ=255)
                                                                    tACCATTAAAGCGGGATCGGCATGGGTAACGTTGATAATCGGCAACACCACACCCGCGGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCAGCACATCAATAACGGCGAAAGTGTAGTTc                                                                              >  1:580370/1‑149 (MQ=255)
                                                                    tACCATAAAAGCGGGATCCGCATGGGTAACGTTGATAATCCGCAACACCACACCCGCGGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCAGCACATCAATAACGGCGAAAGTGTAGTTc                                                                              <  1:144729/149‑1 (MQ=255)
                                                                     aCCATTAAAGCGGGATCGGCATGGGTAACGTTGATAATCGGCAACACCACACCCGCGGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCAGCACATCAATAACGGCCAAAGTGTAGTTcg                                                                             >  1:244234/1‑148 (MQ=255)
                                                                      ccATTAAAGCGGGATCGGCAGGGGTAACGTTGATAATCGGCAACACCACACCCGCGGTGGTAATGGCCGCGACGGTGGCAGCGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCAGCACATCAATAACGGCGAAAGTGTAGTTCCg                                                                            <  2:1239856/149‑1 (MQ=255)
                                                                       cATTAAAGCGGGATCGGCATGGGTAACGTTGATAATCGGCAACACCACACCCGCGGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCAGCACATCAATAACGGCGAAAGTGTAGTTCCgg                                                                           >  1:1182286/1‑149 (MQ=255)
                                                                         ttAAAGCGGGATCGGCATGGGTAACGTTGATAATCGGCAACACCACACCCGCGGTGGTACTGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCAGCACATCAATAACGGCGAAAGTGTAGTTCCGGTc                                                                         >  1:160108/1‑149 (MQ=255)
                                                                         ttAAAGCGGGATCGGCATGGGTAACGTTGATAATCGGCAACACCACACCCGCGGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCAGCACATCAATAACGGCGAAAGTGTAGTACCGGTc                                                                         >  2:515120/1‑149 (MQ=255)
                                                                               cGGGATCGGCATGGGTAACGTTGATAATCGGCAACACCACCCCCGCGGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCAGCACATCAATAACGGCGAAAGGGTAGTTCCGGTCAGTAAg                                                                   >  1:582841/1‑149 (MQ=255)
                                                                                  gaccGGCATGGGTAACGTTGATAACCGGCAACACCACCCCCGCGGGGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCAGCACATCAATAACGGCGAAAGTGTAGTTCCGGTCATTAAgtgt                                                                <  1:1448042/146‑1 (MQ=255)
                                                                                          tGGGTAACGTTGATAATCGGCAACACCACACCCGCGGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTGCAGCACATCAATAACAGCGAAAGTGTAGGTCCCGTCATTAAGTGGAAAAtata                                                        >  2:614333/1‑149 (MQ=255)
                                                                                                  gtggATAACCGGCACCCCCCCCCCCGCGGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCCACAGTCCAGCACATCAATAACGGCGAAAGTGTAGTTCCGGTCATTAAGTGTGAAATATAGTGCCCGa                                                <  2:791953/146‑1 (MQ=255)
                                                                                                  gTTGATACTCGGCAACACCACACCCGCGGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCAGCACATCAATAACGGCGAAAGTGTAGTTCCGGTCATTAAGTGTGAAATATAGTGCCCGa                                                <  1:606139/149‑1 (MQ=255)
                                                                                                  gTTGATAATCGGCAACACCACACCCGCGGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCAGCACATCAATAACGGCGAAAGTGTAGTTCCGGTCATTAAGTGTGAAATATAGTGCCCGa                                                >  1:1038291/1‑149 (MQ=255)
                                                                                                    tGATAATCGGCAACACCACACCCGCGGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACCGCGCCGCAACAGTCCAGCACATCAATAACGGCGAAAGTGTAGTTCCGGTCATTAAGTGTGAAATATAGGGCCCGtca                                              >  2:1294098/1‑149 (MQ=255)
                                                                                                       tAATCGGCACCACCACACCCGCGGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCAGCACATCAATAACGGCGAAAGTGTAGTTCCGGTCATTAAGTGTGAAATATAGTGCCCGACACCg                                           <  1:963246/149‑1 (MQ=255)
                                                                                                       tAATCGGCAACACCACACCCGCGGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCAGCACATCAATAACGGCGAAAGTGTAGTTCCGGTCATTAAGTGTGAAATATAGTGCCCGACACCg                                           >  2:1442423/1‑149 (MQ=255)
                                                                                                       tAATCGGCAACACCACACCCGCGGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCAGCACATCAATAACGGCGAAAGAGTAGTTCCGGTCCTCAAGTGGGCAATGTCGTGGCCCGCACCg                                           >  2:1276141/1‑149 (MQ=255)
                                                                                                       tAATCGGCAACACCACACCCGCGGGTGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCAGCACATCAATAACGGCGAAAGTGTAGTTCCGGTCATTAAGTGTGAAATATAGTGCCCGACACCg                                           <  1:887765/149‑1 (MQ=255)
                                                                                                        aaTCGGCAACACCACACCCGCGGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCAGCACATCAATAACGGCGAAAGTGTAGATCCGGCCATTAAGAGTGAAATATAGTGCCCGACACAGc                                          >  2:293843/1‑149 (MQ=255)
                                                                                                        aaTCGGCAACACCACACCCGCGGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCAGCACATCAATAACGGCGAAAGAGGAGTTCCGGCCATTAAGTGTGAAATATAGTGCCCGACACCGc                                          >  2:153765/1‑149 (MQ=255)
                                                                                                          tCGGCAACACCACACCCGCGGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCAGCACATCAATAACGGCGAAAGTGTAGTTCCGGTCATTAAGTGTGAAATATAGTGCCCGACCCCGCTa                                        >  1:502095/1‑149 (MQ=255)
                                                                                                           cGGCAACACCACCCCCGCGGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCAGCACATCAATAACGGCGAAAGTGCAGTTCCGGTCATTAAGTGTGAAATATAGTGCCCGACACCGCTAt                                       <  2:424571/149‑1 (MQ=255)
                                                                                                            ggCAACACCACACCCGCGGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCAGCACATCAATAACGGCGAAAGTGTAGTTCCGGTCATTAAGTGTGAAATATAGTGCCCGACACCGCTacc                                      >  2:459456/1‑147 (MQ=255)
                                                                                                              cAACACCACCCCCGCGGTGGTAAGGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCAGCACATCAATAACGGCGAAAGTGTAGCTCCGGCCATTAAGTGTGAAATATAGTGCCCGACACCGCTATcta                                    <  2:375391/149‑1 (MQ=255)
                                                                                                                acaccacaccCGCAGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCAGCACATCAATAACGGCGAAAGTGTAGTTCCGGTCATTAAGTGTGAAATATAGTGCCCGACACCGCTATctact                                  <  1:1187851/149‑1 (MQ=255)
                                                                                                                 caccacaccCGCGGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCAGCACATCAATAACGGCAAAAGTGTAGTTCCGGTCATTAAGTGTGAAATATAGTGCCCGACACCGCTATctacta                                 >  2:774625/1‑149 (MQ=255)
                                                                                                                  accacaccCGCGGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCAGCACATCAATAACGGCGAAAGTGTAGTTCCGGTCATTAAGTGTGAAATATAGTGCCCGACACCGCTATCTACTaa                                >  1:1103720/1‑149 (MQ=255)
                                                                                                                   ccacaccCGCGGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCAGCACATCAATAACGGCGAGAGTGGAGTTCCGATCATTAAGTGCGAAATATAGAGGCCGACACCGCTATCTACTAAt                               >  2:963175/1‑149 (MQ=255)
                                                                                                                   ccacaccCGCGGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCAGCACATCAATAACGGCGAAAGTGTAGTTCCGGTCATTAAGTGTGAAATATAGTGCCCGACACCGCTATCTACTAAt                               >  1:1445880/1‑149 (MQ=255)
                                                                                                                     ccaccGGCGGTGGTAAGGGCCGCGCGGGTGGCAGAGCCCAGAGCGCTACGCAACAGCGCCGCAACAGTCCAGCACATCAATAGCGGCGAAAGTGTAGTTCCGGTCATTAAGTGTGAAATATAGTGCCCGACACCGCTATCTACTAATAc                             <  1:98290/148‑1 (MQ=255)
                                                                                                                        cccGCGGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCAGCACATCAATAACGGCGAAAGTGTAGTTCCGGTCATTAAGTGTGAAATATAGTGCCCGACACCGCTATCTACTAATACCTg                          <  1:1007446/149‑1 (MQ=255)
                                                                                                                         ccGCGGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCAGCACATCAATAACGGCGGAAACCCAGATCCGGTCATTAAGTGTGAAATATAGTGCCCGACACCGCTATCTACTAATACCTGc                         >  1:1427095/1‑149 (MQ=255)
                                                                                                                         ccGCGGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCAGCACATCAATAACGGCGAAAGTGTAGTTCCGGTCATTAAGTGTGAAATATAGTGCCCGACACCGCTATCTACTAATACCTGc                         <  1:517066/149‑1 (MQ=255)
                                                                                                                         ccGCGGGGGGAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATCCGACACAGCGCTGCAACAGTCCAGCACATCAATAACGGCGAAAGTGTAGTTCCGGTCATTAAGTGTGAAATATAGTGCCCGACACCGCTATCTACTAATACCTGc                         <  1:334994/149‑1 (MQ=255)
                                                                                                                          cgcgGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCAGCACATCAATAACGGCGAAAGTGTAGTTCCGGTCATTAAGTGTGAAAAATAGTGCCCGACACCGCTATCTACTAATACCTGCt                        >  1:413860/1‑149 (MQ=255)
                                                                                                                             ggtggtAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCAGCACATCAATAACGGCGAAAGTGTAGTTCCGGTCATTAAGTGTGAAATATAGTGCCCGACACCGCTATCTACTAATACCTGCTTaa                     >  1:30698/1‑149 (MQ=255)
                                                                                                                                    aTGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCAGCACATCAATAACGGCGAAAGTGTAGTTCCGGTCATTAAGTGTGAAATATAGTGCCCGACACCGCTATCTACTAATACCTGCTTAAACGcgcc              <  2:1194645/149‑1 (MQ=255)
                                                                                                                                      ggCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCAGCACATCAATAACGGCGAAAGTGTAGTTCCGGTCATTAAGTGTGAAATATAGTGCCCGACACCGCTATCTACTAATACCTGCTTAAACGcgccgc            <  1:1317414/149‑1 (MQ=255)
                                                                                                                                          gcgACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCAGCACATCAATAACGGCGAAAGTGTAGTTCCGGTCATTAAGTGTGAAATATAGTGCCCGACACCGCTATCTACTAATACCTGCTTAAACGCgccgccgcc        >  1:482431/1‑149 (MQ=255)
                                                                                                                                           cgACGGTGGCAGAGCCCAAAGCCAAACCCAACAGCGCCGCAACAGTCCAGCACATCAACAACGGCGAAAGTGTAGTTCCGGTCATTAAGTGTGAAATATAGCGCCCGACAACCCTATCTACTAATACCTGCTTAAACGCGCCGCCGCCa       <  2:314515/149‑1 (MQ=255)
                                                                                                                                             aCGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCAGCACATCAATAACGGCGAAAGTGTAGTTCCGGTCATTAAGTGTGAAATATAGTGCCCGACACCGCTATCTACTAATACCTGCTTAAACGCGCCGCCGCCAGc     <  1:327327/149‑1 (MQ=255)
                                                                                                                                               ggTGGCCGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCAGCACATCAATAACGGCGAAAGTGTAGTTCCGGTCATTAAGTGTGAAATATAGTGCCCGACACCGCTATCTACTAATACCTGCTTAAACGCGCCGCCGCCAGCGa   <  1:19031/149‑1 (MQ=255)
                                                                                                                                                gTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCAGCACATCAATAACGGCGAAAGTGTAGTTCCGGTCATTAAGTGTGAAATATAGTGCCCGACACCGCTATCTACTAATACCTGCTTAAACGCGCCGCCGCCAGCGat  >  1:151732/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CCTTTAAATAGCCAGAAGCCAGGGTCGTTTACGTGTGACGCGATCACGCTGCCCGCACCGGTTGCCAGTACCATTAAAGCGGGATCGGCATGGGTAACGTTGATAATCGGCAACACCACACCCGCGGTGGTAATGGCCGCGACGGTGGCAGAGCCCAGAGCGATACGCAACAGCGCCGCAACAGTCCAGCACATCAATAACGGCGAAAGTGTAGTTCCGGTCATTAAGTGTGAAATATAGTGCCCGACACCGCTATCTACTAATACCTGCTTAAACGCGCCGCCGCCAGCGAT  >  NC_000913/4491319‑4491611

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