Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,309,875 T→G 6.0% K369T (AAA→ACA) ‡ alsA ← D‑allose ABC transporter ATP binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,309,8750TG 6.0% 355.3 / 12.9 116K369T (AAA→ACA) ‡alsAD‑allose ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base T (55/54);  new base G (3/4);  total (58/58)
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.01e-01

Only 100 of 120 total aligned reads displayed.

GGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGATATCTTTGCCATTAAGACGGATT  >  NC_000913/4309726‑4310024
                                                                                                                                                     |                                                                                                                                                     
ggACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGATGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCtt                                                                                                                                                       >  1:1285749/1‑150 (MQ=255)
  aCTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTc                                                                                                                                                     >  2:1250047/1‑150 (MQ=255)
  aCTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTc                                                                                                                                                     <  1:1045706/150‑1 (MQ=255)
  aCTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGACATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTc                                                                                                                                                     <  1:1324338/150‑1 (MQ=255)
   cTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCa                                                                                                                                                    >  1:1023087/1‑150 (MQ=255)
   cTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCa                                                                                                                                                    <  2:362804/150‑1 (MQ=255)
       ctgctgATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACt                                                                                                                                                >  1:70957/1‑150 (MQ=255)
       ctgctgATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACAGAATGACATTTCAGCGAAAGCAGTTAGCGTTGATTTTAAGCGGTACGTTGCTCTTCAAATTCATGAAACAACCCCATCGCGCCTTTATAGACGCCGTCTTTAAGACt                                                                                                                                                >  2:591049/1‑150 (MQ=255)
            gATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGc                                                                                                                                           <  1:583623/150‑1 (MQ=255)
            gATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCAGTCTGAAAGACggcggc                                                                                                                                           <  2:613389/150‑1 (MQ=255)
             aTTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGTGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCAGTCTGAAAGACGGCGGCt                                                                                                                                          <  1:205217/150‑1 (MQ=255)
                tCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGAt                                                                                                                                       >  1:306561/1‑150 (MQ=255)
                tCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCAGTCTGAAAGACGGCGGCtata                                                                                                                                       <  1:1162075/150‑4 (MQ=255)
                      ggAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCAt                                                                                                                                 <  1:917013/150‑1 (MQ=255)
                        agagTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGt                                                                                                                               <  1:11067/150‑1 (MQ=255)
                          agTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTc                                                                                                                             <  1:1070101/150‑1 (MQ=255)
                             tCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGa                                                                                                                          >  1:981222/1‑150 (MQ=255)
                              cGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAg                                                                                                                         >  2:110087/1‑150 (MQ=255)
                                gTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCg                                                                                                                       >  1:1018161/1‑150 (MQ=255)
                                gTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCg                                                                                                                       >  2:1217902/1‑150 (MQ=255)
                                    tatTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCa                                                                                                                                                    <  1:370954/117‑1 (MQ=255)
                                    tatTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCa                                                                                                                                                    >  2:370954/1‑117 (MQ=255)
                                     atTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGa                                                                                                                  >  2:277908/1‑150 (MQ=255)
                                     atTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGa                                                                                                                  <  1:1369623/150‑1 (MQ=255)
                                      ttCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGaa                                                                                                                 >  2:1318868/1‑150 (MQ=255)
                                         tGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAg                                                                                                              <  1:125454/150‑1 (MQ=255)
                                            tttACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCTCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTg                                                                                                           <  1:320375/150‑1 (MQ=255)
                                              tACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTggg                                                                                                         >  1:439601/1‑150 (MQ=255)
                                              tACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAACCCCATCGAGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTggg                                                                                                         >  2:1448808/1‑150 (MQ=255)
                                                ccGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGa                                                                                                       >  2:1027973/1‑150 (MQ=255)
                                                ccGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGa                                                                                                       >  1:461615/1‑150 (MQ=255)
                                                  gAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGCGCGATGGAAAAGTTGGGGaaa                                                                                                     >  2:438393/1‑150 (MQ=255)
                                                      gACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAg                                                                                                                         <  1:1036544/126‑1 (MQ=255)
                                                      gACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAg                                                                                                                         >  2:1036544/1‑126 (MQ=255)
                                                      gACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAAc                                                                                                 >  2:589185/1‑150 (MQ=255)
                                                      gACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAAc                                                                                                 >  2:1073169/1‑150 (MQ=255)
                                                        cATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCg                                                                                               >  1:1299685/1‑150 (MQ=255)
                                                        cATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCg                                                                                               <  1:468953/150‑1 (MQ=255)
                                                           ttCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTa                                                                                            <  1:613954/150‑1 (MQ=255)
                                                            tCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTAt                                                                                           >  1:279122/1‑150 (MQ=255)
                                                             cAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATc                                                                                          >  1:800091/1‑150 (MQ=255)
                                                               gcgcCAGCATTTCGCGTTGTTTTTCAGCGGTACGTTGGTCGTCAACTTCATGAAACAAGCCCATCGCGGCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATccc                                                                                        <  1:1508309/150‑1 (MQ=255)
                                                                 gcCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCcgg                                                                                      >  1:606542/1‑150 (MQ=255)
                                                                   cagcagTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCcggcg                                                                                    >  1:284401/1‑150 (MQ=255)
                                                                    agcagTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCcggcgg                                                                                   <  2:1299685/150‑1 (MQ=255)
                                                                      cagTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCt                                                                                 <  1:346951/150‑1 (MQ=255)
                                                                      cagTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCt                                                                                 <  1:1073169/150‑1 (MQ=255)
                                                                             cgTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTg                                                                          <  2:1285749/150‑1 (MQ=255)
                                                                              gTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGa                                                                         >  1:1317971/1‑150 (MQ=255)
                                                                              gTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGa                                                                         <  2:279122/150‑1 (MQ=255)
                                                                               ttGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACTGTTATCCCGGCGGCTTTCAGTGAt                                                                        <  2:1023087/150‑1 (MQ=255)
                                                                               ttGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGAt                                                                        >  2:1003718/1‑150 (MQ=255)
                                                                                tGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATTCCGGCGGCTTTCAGTGATg                                                                       <  2:1121500/150‑1 (MQ=255)
                                                                                    tttCAGCGGTACGTTGCTCGTCAACTTCATGAAACTAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTTATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAg                                                                   <  1:591049/150‑1 (MQ=255)
                                                                                         gCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGa                                                                         >  2:131398/1‑139 (MQ=255)
                                                                                         gCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGa                                                                         >  2:125488/1‑139 (MQ=255)
                                                                                         gCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGa                                                                         <  1:125488/139‑1 (MQ=255)
                                                                                           ggTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATccc                                                            <  2:1317971/150‑1 (MQ=255)
                                                                                             tACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCtt                                                          >  2:673882/1‑150 (MQ=255)
                                                                                                gTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCttttt                                                       >  2:865124/1‑150 (MQ=255)
                                                                                                    cTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCAc                                                   <  1:891618/150‑1 (MQ=255)
                                                                                                    cTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCAc                                                   >  2:339976/1‑150 (MQ=255)
                                                                                                    cTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCAc                                                   >  2:340899/1‑150 (MQ=255)
                                                                                                     tCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACg                                                  <  1:1003718/150‑1 (MQ=255)
                                                                                                     tCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCCGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACg                                                  <  1:110087/150‑1 (MQ=255)
                                                                                                      cGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGCCTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACgg                                                 >  1:1332513/1‑150 (MQ=255)
                                                                                                      cGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCACGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACgg                                                 >  1:1225319/1‑150 (MQ=255)
                                                                                                           aCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATcc                                            >  2:266841/1‑150 (MQ=255)
                                                                                                           aCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATcc                                            <  2:357785/150‑1 (MQ=255)
                                                                                                            cTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCa                                           >  1:1178174/1‑150 (MQ=255)
                                                                                                               cATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAggg                                        >  2:73090/1‑150 (MQ=255)
                                                                                                                aTGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAgggg                                       <  2:1568936/150‑1 (MQ=255)
                                                                                                                aTGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAgggg                                       <  1:865124/150‑1 (MQ=255)
                                                                                                                 tGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGt                                      >  1:1384015/1‑150 (MQ=255)
                                                                                                                         gCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTgg                              <  2:298421/150‑1 (MQ=255)
                                                                                                                         gCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTgg                              >  1:527982/1‑150 (MQ=255)
                                                                                                                          cccATCGCGCCTTTATAGCCGCCTTGTTTCAGACTGCGGCTGATCGGGATGTTATGAGCGATGGAAAAGTTGGGGAAAAAAGCGTTATCCCGGCGGCTTTCAGTGATTTAAGCCATCCCTTTTTTCACGGCATCCAGGGGGGAACGTGga                             <  2:800091/150‑1 (MQ=255)
                                                                                                                          cccATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGga                             >  1:1160125/1‑150 (MQ=255)
                                                                                                                          cccATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGga                             >  1:621562/1‑150 (MQ=255)
                                                                                                                           ccATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGgag                            <  2:1114145/150‑1 (MQ=255)
                                                                                                                           ccATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGgag                            >  1:1218015/1‑150 (MQ=255)
                                                                                                                             aTCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGgagag                          <  2:1178174/150‑1 (MQ=255)
                                                                                                                             aTCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGgagag                          <  1:1348287/150‑1 (MQ=255)
                                                                                                                                gcgcCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTTGTGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACTGCATCCAGTTGTTAAAGTGGAGAtata                       <  2:306561/150‑1 (MQ=255)
                                                                                                                                gcgcCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGata                       >  2:452482/1‑150 (MQ=255)
                                                                                                                                  gcCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGATATc                     <  1:1318868/150‑1 (MQ=255)
                                                                                                                                   ccTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACAGCATCCAGGGGTGAAc                                  >  1:763243/1‑136 (MQ=255)
                                                                                                                                   ccTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACAGCATCCAGGGGTGAAc                                  <  2:763243/136‑1 (MQ=255)
                                                                                                                                      ttATAGCCGCAGTATTTCAGACTGGGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGATATCTTTg                 <  2:1585029/150‑1 (MQ=255)
                                                                                                                                       tataGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGGGATGTAAGCCATCCCTTTTTTCACGGCATCCAGTGGTGAACGTGGAGAGATATCTTTGc                <  2:802885/150‑1 (MQ=255)
                                                                                                                                           gccgccGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGATATCTTTGCCAtt            >  1:142366/1‑150 (MQ=255)
                                                                                                                                           gccgcaGTCTGAAAGACGGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGATATCTTTGCCAtt            >  2:324042/1‑150 (MQ=255)
                                                                                                                                           gccgcaGTCTGAAAGACGGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGATATCTTTGCCAtt            >  2:1418507/1‑150 (MQ=255)
                                                                                                                                             cgccgTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTTAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGATATCTTTGCCATTat          <  1:1309295/150‑2 (MQ=255)
                                                                                                                                              gccgTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTg                                                                          <  2:1242672/85‑1 (MQ=255)
                                                                                                                                                 gTCTTTCAGACTGCGGCTGATCGCGATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGATATCTTTGCCATTAAGACg      <  2:95201/150‑1 (MQ=255)
                                                                                                                                                 gTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGATATCTTTGCCATTAAGACg      <  1:381737/150‑1 (MQ=255)
                                                                                                                                                 gTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGATATCTTTGCCATTAAGACg      <  1:1095944/150‑1 (MQ=255)
                                                                                                                                                    tttCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGATATCTTTGCCATTAAGACGGAt   >  1:604890/1‑150 (MQ=255)
                                                                                                                                                    tttCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGATATCTTTGCCATTAAGACGGAt   >  2:1101369/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGATATCTTTGCCATTAAGACGGATT  >  NC_000913/4309726‑4310024

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