Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,309,875 T→G 6.2% 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.2% 298.9 / 10.7 97K369T (AAA→ACA) ‡alsAD‑allose ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base T (46/45);  new base G (2/4);  total (48/49)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.44e-01

GGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGATATCTTTGCCATTAAGACGGATT  >  NC_000913/4309726‑4310024
                                                                                                                                                     |                                                                                                                                                     
ggACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCtt                                                                                                                                                       >  2:1628526/1‑150 (MQ=255)
ggACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCtt                                                                                                                                                       <  1:1371254/150‑1 (MQ=255)
  aCTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTc                                                                                                                                                     >  1:535687/1‑150 (MQ=255)
   cTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCa                                                                                                                                                    <  2:1664469/150‑1 (MQ=255)
      tCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGAc                                                                                                                                                 >  1:360842/1‑150 (MQ=255)
          gtgATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCg                                                                                                                                             <  2:1532980/149‑1 (MQ=255)
          ctgATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCg                                                                                                                                             >  2:802243/1‑150 (MQ=255)
           tgATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCgg                                                                                                                                            >  2:1090939/1‑150 (MQ=255)
              tttCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTg                                                                                                                                         >  1:78543/1‑150 (MQ=255)
              tttCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTg                                                                                                                                         <  1:1311218/150‑1 (MQ=255)
              tttCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCAGTCTGAAAGACggcggc                                                                                                                                           <  1:416101/148‑1 (MQ=255)
              tttCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCAGTCTGAAAGACggcggc                                                                                                                                           >  2:416101/1‑148 (MQ=255)
                tCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCAGTCTGAAAGACGGCGGCtata                                                                                                                                       <  2:127914/150‑4 (MQ=255)
                 cccccGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATc                                                                                                                                      >  1:951400/1‑150 (MQ=255)
                  ccccGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCg                                                                                                                                     >  2:363055/1‑150 (MQ=255)
                  ccccGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCg                                                                                                                                     >  1:712285/1‑150 (MQ=255)
                  ccccGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCAGTCTGAAAGACGGCGGCtataaa                                                                                                                                     <  2:762847/150‑6 (MQ=255)
                           gTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCt                                                                                                                            <  1:1303711/150‑1 (MQ=255)
                           gTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCt                                                                                                                            >  2:1289586/1‑150 (MQ=255)
                             tCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGa                                                                                                                          >  1:287543/1‑150 (MQ=255)
                              cGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAg                                                                                                                         <  1:647479/150‑1 (MQ=255)
                              cGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGACTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAg                                                                                                                         >  2:196508/1‑150 (MQ=255)
                                gTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCg                                                                                                                       >  2:1415601/1‑150 (MQ=255)
                                 tGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGAGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGa                                                                                                                      >  1:714708/1‑150 (MQ=255)
                                  gATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGAt                                                                                                                     <  1:363055/150‑1 (MQ=255)
                                   atatTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATg                                                                                                                    <  2:763191/150‑1 (MQ=255)
                                   atatTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATg                                                                                                                    <  2:1539288/150‑1 (MQ=255)
                                       tCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGaaa                                                                                                                <  2:57026/150‑1 (MQ=255)
                                        cTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGaaaa                                                                                                               >  2:378835/1‑150 (MQ=255)
                                        cTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGaaaa                                                                                                               >  2:378839/1‑150 (MQ=255)
                                          ggTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGt                                                                                                             >  1:1415096/1‑150 (MQ=255)
                                            tttACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTg                                                                                                           <  1:367688/150‑1 (MQ=255)
                                             ttACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTgg                                                                                                          <  1:608083/150‑1 (MQ=255)
                                             ttACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTgg                                                                                                          <  2:535687/150‑1 (MQ=255)
                                              tACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTggg                                                                                                         <  1:1380869/150‑1 (MQ=255)
                                              tACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTggg                                                                                                         <  1:1179962/150‑1 (MQ=255)
                                               aCCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCATCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTgggg                                                                                                        <  2:170379/150‑1 (MQ=255)
                                                  gAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGaaa                                                                                                     >  1:742552/1‑150 (MQ=255)
                                                  gAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGaaa                                                                                                     >  1:1162920/1‑150 (MQ=255)
                                                     tGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGaaaaaa                                                                                                  >  2:326104/1‑150 (MQ=255)
                                                     tGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGaaaaaa                                                                                                  <  1:572227/150‑1 (MQ=255)
                                                       aCATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAAcc                                                                                                >  1:671584/1‑150 (MQ=255)
                                                        cATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCg                                                                                               >  1:982163/1‑150 (MQ=255)
                                                               gcgcCAGCAGTTCGCCATGATTATCAGCGGTACGTTGCTCGTAAACTTAATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTTCGGCTGATCGCAATGTTCTGAGCGATGGAAAAGTTGGGTAAAAAACAGATATccc                                                                                        >  2:841909/1‑150 (MQ=255)
                                                                 gcCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCcgg                                                                                      >  1:1090170/1‑150 (MQ=255)
                                                                    agcagTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCcggcgg                                                                                   <  1:1628526/150‑1 (MQ=255)
                                                                    agcagTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCcggcgg                                                                                   <  1:765895/150‑1 (MQ=255)
                                                                      cagTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCt                                                                                 <  2:221220/150‑1 (MQ=255)
                                                                        gTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCttt                                                                               >  1:1424136/1‑150 (MQ=255)
                                                                         ttCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTc                                                                              <  2:742552/150‑1 (MQ=255)
                                                                            gcgTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGt                                                                           >  1:251853/1‑150 (MQ=255)
                                                                            gcgTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGt                                                                           >  1:307119/1‑150 (MQ=255)
                                                                             cgTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTg                                                                          <  1:1516603/150‑1 (MQ=255)
                                                                             cgTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTg                                                                          <  2:507038/150‑1 (MQ=255)
                                                                                tGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATg                                                                       <  2:1053733/150‑1 (MQ=255)
                                                                                 gATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGt                                                                      <  2:951400/150‑1 (MQ=255)
                                                                                 gATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGt                                                                      <  2:1415096/150‑1 (MQ=255)
                                                                                 gATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGt                                                                      >  2:709362/1‑150 (MQ=255)
                                                                                    tttCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAg                                                                   <  2:681794/150‑1 (MQ=255)
                                                                                     ttCAGCGATACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGc                                                                  >  2:879236/1‑150 (MQ=255)
                                                                                      tCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGCCTTTCAGAGATCTAAGcc                                                                 >  2:1389806/1‑150 (MQ=255)
                                                                                       cAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCa                                                                >  1:250037/1‑150 (MQ=255)
                                                                                         gCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATc                                                              <  1:922227/150‑1 (MQ=255)
                                                                                                gTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCttttt                                                       <  2:360842/150‑1 (MQ=255)
                                                                                                  tGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCtttttta                                                     >  2:710262/1‑149 (MQ=255)
                                                                                                    cTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCAc                                                   <  2:1090170/150‑1 (MQ=255)
                                                                                                    cTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCAc                                                   <  2:373804/150‑1 (MQ=255)
                                                                                                        tCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCa                                               >  1:877208/1‑150 (MQ=255)
                                                                                                         cAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCAt                                              >  2:1586431/1‑150 (MQ=255)
                                                                                                         cAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCAt                                              <  1:1359394/150‑1 (MQ=255)
                                                                                                          aaCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATc                                             >  1:1174510/1‑150 (MQ=255)
                                                                                                           aCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCa                                               >  1:116554/1‑147 (MQ=255)
                                                                                                           aCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCa                                               <  2:116554/147‑1 (MQ=255)
                                                                                                              tCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAgg                                         <  1:823398/150‑1 (MQ=255)
                                                                                                               cATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAggg                                        >  2:902095/1‑150 (MQ=255)
                                                                                                               cATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAggg                                        >  1:1656113/1‑150 (MQ=255)
                                                                                                                aTGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAgggg                                       <  1:879236/150‑1 (MQ=255)
                                                                                                                 tGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGt                                      <  2:78543/150‑1 (MQ=255)
                                                                                                                 tGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGACATCCCTTTTTTAACGTCATCCAgggtt                                      >  2:744740/1‑148 (MQ=255)
                                                                                                                 tGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGCCGCATTTCAGTGATGTAAGCAATCCATTTTTAAACGGCATCCAGGGGt                                      >  2:188501/1‑150 (MQ=255)
                                                                                                                  gAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTg                                     >  1:1250736/1‑150 (MQ=255)
                                                                                                                  gAAACAAGCCCATCGCGCCTTAATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTg                                     <  2:1070128/150‑1 (MQ=255)
                                                                                                                      cAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACg                                 <  1:966604/150‑1 (MQ=255)
                                                                                                                        aGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTg                               >  2:1463288/1‑150 (MQ=255)
                                                                                                                          cccATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGga                             <  1:1289586/150‑1 (MQ=255)
                                                                                                                          cccATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGga                             <  1:981225/150‑1 (MQ=255)
                                                                                                                             aTCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGgagag                          >  1:1500821/1‑150 (MQ=255)
                                                                                                                               cgcgCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGat                        <  2:287543/150‑1 (MQ=255)
                                                                                                                                gcgcCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGata                       <  2:307119/150‑1 (MQ=255)
                                                                                                                                      ttATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGATATCTTTg                 <  2:1656113/150‑1 (MQ=255)
                                                                                                                                       tataGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAATCCATCCCTTTTTTCACGCCATCCAGTGGTGAACGTGGAGAGATATCTTTGc                >  2:719227/1‑150 (MQ=255)
                                                                                                                                        ataGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGATATCTTTGcc               <  1:802243/150‑1 (MQ=255)
                                                                                                                                          aGCCGCAGTCTGAAAGACGGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCttt                                                         >  1:270161/1‑106 (MQ=255)
                                                                                                                                          aGCCGCAGTCTGAAAGACGGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCttt                                                         <  2:270161/106‑1 (MQ=255)
                                                                                                                                            ccgccgGCTTTCAGACTGCGGCTGATCGCGATGTTCTTAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGATATCTTTGCCATTa           <  2:750152/150‑1 (MQ=255)
                                                                                                                                             cgccgTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGATATCTTTGCCATTaa          >  1:812457/1‑150 (MQ=255)
                                                                                                                                                   cTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAgggg                                       >  1:1169358/1‑115 (MQ=255)
                                                                                                                                                   cTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAgggg                                       <  2:1169358/115‑1 (MQ=255)
                                                                                                                                                   cTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGATATCTTTGCCATTAAGACGGa    >  1:1006761/1‑150 (MQ=255)
                                                                                                                                                     ttCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGATATCTTTGCCATTAAGACGGAtt  <  1:1090939/150‑1 (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: