Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,897,538 A→G 6.0% I111I (ATT→ATC ygcS ← putative MFS sugar transporter; membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,897,5380AG 6.0% 209.4 / 2.6 67I111I (ATT→ATCygcSputative MFS sugar transporter; membrane protein
Reads supporting (aligned to +/- strand):  ref base A (29/34);  new base G (2/2);  total (31/36)
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.80e-01

ACTTGCCAGCACATAGCCTACGGTCCACACCACGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTGCCAAGGAACAAACCAAGGAGCGCCG  >  NC_000913/2897391‑2897686
                                                                                                                                                   |                                                                                                                                                    
aCTTGCCAGCACATAGCCTACGGTCCACACCACGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCaa                                                                                                                                                     >  2:334784/1‑149 (MQ=255)
 cTTGCCAGCACATAGCCTACGGTCCACACCACGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAAt                                                                                                                                                    <  7:315500/149‑1 (MQ=255)
          acaTAGCCTACGGTCCACACCACGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGctc                                                                                                                                           >  4:196269/1‑149 (MQ=255)
           caTAGCCTACGGTCCACACCACGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCGATAAGATGctct                                                                                                                                          >  1:59242/1‑149 (MQ=255)
           caTAGCCTACGGTCCACACCACGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCGATAAGATGctct                                                                                                                                          >  1:53846/1‑149 (MQ=255)
                 cTACGGTCCACACCACGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCgtg                                                                                                                                    >  6:531/1‑149 (MQ=255)
                     ggTCCACACCACGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGtggtgg                                                                                                                                <  1:43790/149‑1 (MQ=255)
                       tCCACACCACGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATACTCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGc                                                                                                                               >  3:280810/1‑148 (MQ=255)
                         cacaCCACGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCaa                                                                                                                             >  5:120233/1‑148 (MQ=255)
                              cacGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAAtt                                                                                                                       <  5:531/149‑1 (MQ=255)
                              cacGCTGAATGCGCCCAACAAAATACCGCGATGGCGGCGCGGGGAAAATTCAACCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCTTGGTGGCAAAAAAtt                                                                                                                       <  5:150859/149‑1 (MQ=255)
                               acGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTg                                                                                                                      >  4:146200/1‑149 (MQ=255)
                                cGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGt                                                                                                                     <  1:334784/149‑1 (MQ=255)
                                cGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGt                                                                                                                     <  1:41309/149‑1 (MQ=255)
                                cGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGt                                                                                                                     <  1:12818/149‑1 (MQ=255)
                                 gCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGCCCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGt                                                                                                                     >  8:326946/1‑148 (MQ=255)
                                     aaTGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTaaaaa                                                                                                                >  8:245417/1‑149 (MQ=255)
                                     aaTGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGACAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGACAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTaaaaa                                                                                                                >  4:243871/1‑149 (MQ=255)
                                      aTGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAAc                                                                                                               >  5:265448/1‑149 (MQ=255)
                                         cgcCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGaa                                                                                                            <  1:291212/149‑1 (MQ=255)
                                                cAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCg                                                                                                     <  2:215742/149‑1 (MQ=255)
                                                      aCCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAAt                                                                                                 >  4:349858/1‑147 (MQ=255)
                                                       ccGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAg                                                                                              >  4:2372/1‑149 (MQ=255)
                                                         gcgATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCa                                                                                            <  4:280810/149‑1 (MQ=255)
                                                                cggcgCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGACAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCt                                                                                     <  3:196269/149‑1 (MQ=255)
                                                                   cgcgGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCGATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGaa                                                                                  <  2:53846/149‑1 (MQ=255)
                                                                   cgcgGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCGATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGaa                                                                                  <  2:59242/149‑1 (MQ=255)
                                                                    gcgGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAg                                                                                 <  7:245417/149‑1 (MQ=255)
                                                                                 agccagcAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGAtttttt                                                                    >  1:302879/1‑149 (MQ=255)
                                                                                 agccagcAAGGTGTGACCTACTGAATAATCGCCTCCCAGAACAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACTAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGAtttttt                                                                    <  3:291577/149‑1 (MQ=255)
                                                                                       cAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAACTGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGAc                                                              <  6:265448/149‑1 (MQ=255)
                                                                                          ggTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGt                                                                                                                     >  2:203109/1‑91 (MQ=255)
                                                                                          ggTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGt                                                                                                                     <  1:203109/91‑1 (MQ=255)
                                                                                           gtgtGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAat                                                          <  3:146200/149‑1 (MQ=255)
                                                                                               gACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATgg                                                      >  5:329843/1‑149 (MQ=255)
                                                                                                   tACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTATGGAGTGGTGGCAAAAAATTGTaaaaa                                                                                                                >  2:288277/1‑87 (MQ=255)
                                                                                                     cTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGga                                                 <  7:326946/148‑1 (MQ=255)
                                                                                                     cTGAATAATCCCCTCCCAGACCAAAGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCgg                                                  <  3:91312/147‑1 (MQ=255)
                                                                                                      tGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCg                                                   >  2:96440/1‑145 (MQ=255)
                                                                                                       gAATAATCGCCTCTCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCg                                                   >  1:274232/1‑144 (MQ=255)
                                                                                                        aataatCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCaaaaaa                                                                                                                         <  6:198214/73‑1 (MQ=255)
                                                                                                        aataatCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCaaaaaa                                                                                                                         >  5:198214/1‑73 (MQ=255)
                                                                                                        aataatCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGgaga                                               <  3:349858/147‑1 (MQ=255)
                                                                                                             tCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAACAAGCTGCAGGTGAAGATTTTTTGCCGACCAATAGGGTCGGAGATCCAcc                                         >  1:98616/1‑148 (MQ=255)
                                                                                                             tCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTc                                                    <  7:69437/137‑1 (MQ=255)
                                                                                                             tCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTc                                                    >  8:69437/1‑137 (MQ=255)
                                                                                                                    ccAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAg                                                                                              <  6:300994/88‑1 (MQ=255)
                                                                                                                    ccAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGaa                                   <  6:329843/147‑1 (MQ=255)
                                                                                                                    ccAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAAcc                                 >  5:91554/1‑149 (MQ=255)
                                                                                                                     cAGACCACTGCCAATCAAACTGCGCAGTCCAATAAGATGCTCTGGCGGGGTGGCAAAAATTTGTAAATACGAAGCAAGCGTAATCAGCAAAACGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCa                                <  2:98616/149‑1 (MQ=255)
                                                                                                                          cAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTg                           >  6:294896/1‑149 (MQ=255)
                                                                                                                           aaTGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTGc                          <  2:302879/149‑1 (MQ=255)
                                                                                                                            aTGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTGc                          >  5:182426/1‑148 (MQ=255)
                                                                                                                               ccAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTGCCaa                       <  3:2372/148‑1 (MQ=255)
                                                                                                                                     aaaaTGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTGCCAAGGAACaa                <  6:91554/149‑1 (MQ=255)
                                                                                                                                     aaaaTGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTGCCAAGGAACa                 >  3:237231/1‑148 (MQ=255)
                                                                                                                                        aTGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTGCCAAGGAACAAAcc             <  3:243871/149‑1 (MQ=255)
                                                                                                                                           cgcAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTGCCAAGGAACAAACCAAg          <  3:171417/149‑1 (MQ=255)
                                                                                                                                            gcAGTCCAATAAGATGCTCTGGCGTTGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTGCCAAGGAACAAACCAAgg         <  6:182426/149‑1 (MQ=255)
                                                                                                                                              aGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCaa                                      >  8:225392/1‑118 (MQ=255)
                                                                                                                                              aGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCaa                                      <  7:225392/118‑1 (MQ=255)
                                                                                                                                              aGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTGCCAAGGAACAAACCAAGGAg       >  2:354633/1‑149 (MQ=255)
                                                                                                                                              aGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTGCCAAGGAACAAACCAAGGAg       >  1:254479/1‑149 (MQ=255)
                                                                                                                                                tCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTc                                                    >  5:268579/1‑102 (MQ=255)
                                                                                                                                                tCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTc                                                    <  6:268579/102‑1 (MQ=255)
                                                                                                                                                  cAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTGCCAAGGAAc                  >  7:108120/1‑134 (MQ=255)
                                                                                                                                                  cAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTGCCAAGGAAc                  <  8:108120/134‑1 (MQ=255)
                                                                                                                                                  cAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTGCCAAGGAACAAACCAAGGAGCGcc   <  2:254479/149‑1 (MQ=255)
                                                                                                                                                   aaTAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTGCCAAGGAACaaa               >  2:189009/1‑136 (MQ=255)
                                                                                                                                                   aaTAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTGCCAAGGAACaaa               <  1:189009/136‑1 (MQ=255)
                                                                                                                                                   aaTAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTGCCAAGGAACAAACCAAGGAGCGcc   >  5:244169/1‑148 (MQ=255)
                                                                                                                                                   aaTAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTGCCAAGGAACAAACCAAGGAGCGCCg  >  6:4626/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ACTTGCCAGCACATAGCCTACGGTCCACACCACGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAAGCTGAAGGTGAAGATTTTTTGCCGACCAATATGGTCGGAGATCCACCCAAGAACCAGGCTGCCAAGGAACAAACCAAGGAGCGCCG  >  NC_000913/2897391‑2897686

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