Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,982,830 A→T 5.1% I237I (ATT→ATA araH ← L‑arabinose ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,982,8300AT 5.1% 336.3 / 4.9 95I237I (ATT→ATAaraHL‑arabinose ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (44/46);  new base T (3/2);  total (47/48)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAAATTAAGATACCCGCCACCACATATGAGATTTTTCCGATGCCACCTTTCAGAGAAACGCCACCTAAAACGCAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAA  >  NC_000913/1982683‑1982977
                                                                                                                                                   |                                                                                                                                                   
tAAAATTAAGATACCCGCCACCACATATGAGATTTTTCCGATGCCACCTTTCAGAGAAACGCCACCTAAAACGCAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGaataa                                                                                                                                                    >  2:99107/1‑149 (MQ=255)
tAAAATTAAGATACCCGCCACCACATATGAGATTTTTCCGATGCCACCTTTCAGAGAAACGCCACCTAAAACGCAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGaataa                                                                                                                                                    >  3:356042/1‑149 (MQ=255)
tAAAATTAAGATACCCGCCACCACATATGAGATTTTTCCGATGCCACCTTTCAGAGAAACGCCACCTAAAACGCAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGaataa                                                                                                                                                    >  1:436616/1‑149 (MQ=255)
tAAAATTAAGATACCCGCCACCACATATGAGATTTTTCCGATGCCACCTTTCAGAGAAACGCCACCTAAAACGCAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGaataa                                                                                                                                                    >  1:458867/1‑149 (MQ=255)
  aaaTTAAGATACCCGCCACCACATATGAGATTTTTCCGATGCCACCTTTCAGAGAAACGCCACCTAAAACGCAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAAtt                                                                                                                                                  <  2:436616/149‑1 (MQ=255)
    aTTAAGATACCCGCCACCACATATGAGATTTTTCCGATGCCACCTTTCAGAGAAACGCCACCTAAAACGCAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTcc                                                                                                                                                <  3:6814/149‑1 (MQ=255)
        aGATACCCGCCACCACATATGAGATTTTTCCGATGCCACCTTTCAGAGAAACGCCACCTAAAACGCAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCt                                                                                                                                            >  1:913658/1‑149 (MQ=255)
         gATACCCGCCACCACATATGAGATTTTTCCGATGCCACCTTTCAGAGAAACGCCACCTAAAACGCAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTa                                                                                                                                           >  2:96591/1‑149 (MQ=255)
          aTACCCGCCACCACATATGAGATTTTTCCGATGCCACCTTTCAGAGAAACGCCACCTAAAACGCAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTAt                                                                                                                                          >  4:153850/1‑149 (MQ=255)
           tACCCGCCACCACATATGAGATTTTTCCGATGCCACCTTTCAGAGAAACGCCACCTAAAACGCAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATc                                                                                                                                         <  1:256807/149‑1 (MQ=255)
             cccGCCACCACATATGAGATTTTTCCGATGCCACCTTTCAGAGAAACGCCACCTAAAACGCAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGc                                                                                                                                       <  1:99107/149‑1 (MQ=255)
              ccGCCACCACATATGAGATTTTTCCGATGCCACCTTTCAGAGAAACGCCACCTAAAACGCAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCt                                                                                                                                      <  2:580111/149‑1 (MQ=255)
                gCCACCACATATGAGATTTTTCCGATGCCACCTTTCAGAGAAACGCCACCTAAAACGCAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGa                                                                                                                                    >  3:649233/1‑149 (MQ=255)
                  caccacATATGAGATTTTTCCGATGCCACCTTTCAGAGAAACGCCACCTAAAACGCAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTGATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATa                                                                                                                                  >  4:698331/1‑149 (MQ=255)
                  caccacATATGAGATTTTTCCGATGCCACCTTTCAGAGAAACGCCACCTAAAACGCAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATa                                                                                                                                  >  2:655805/1‑149 (MQ=255)
                   accacATATGAGATTTTTCCGATGCCACCTTTCAGAGAAACGCCACCTAAAACGCAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATAc                                                                                                                                 <  2:732191/149‑1 (MQ=255)
                      acaTATGAGATTTTTCCGATGCCACCTTTCAGAGAAACGCCACCTAAAACGCAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAg                                                                                                                              <  2:45933/149‑1 (MQ=255)
                            gagaTTTTTCCGATGCCACCTTTCAGAGAAACGCCACCTAAAACGCAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTga                                                                                                                        <  3:909295/149‑1 (MQ=255)
                            gagaTTTTTCCGATGCCACCTTTCAGAGAAACGCCACCTAAAACGCAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTga                                                                                                                        <  1:96591/149‑1 (MQ=255)
                                 ttttCCGATGCCACCTTTCAGAGAAACGCCACCTAAAACGCAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGaa                                                                                                                   <  3:53920/149‑1 (MQ=255)
                                     ccGATGCCACCTTTCAGAGAAACGCCACCTAAAACGCAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACaaa                                                                                                               <  3:580888/149‑1 (MQ=255)
                                       gATGCCACCTTTCAGAGAAACGCCACCTAAAACGCAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGa                                                                                                             <  2:681546/149‑1 (MQ=255)
                                          gCCACCTTTCAGAGAAACGCCACCTAAAACGCAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGataa                                                                                                          <  3:328452/149‑1 (MQ=255)
                                           ccACCTTTCAGAGAAACGCCACCTAAAACGCAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGataat                                                                                                         >  1:561787/1‑149 (MQ=255)
                                             aCCTTTCAGAGAAACGCCACCTAAAACGCAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGAtaataa                                                                                                       <  1:299518/149‑1 (MQ=255)
                                              ccTTTCAGAGAAACGCCACCTAAAACGCAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATaataat                                                                                                      >  1:606043/1‑149 (MQ=255)
                                                tttCAGAGAAACGCCACCTAAAACGCAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAAttt                                                                                                    >  1:429128/1‑149 (MQ=255)
                                                tttCAGAGAAACGCCACCTAAAACGCAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAAttt                                                                                                    >  1:175538/1‑149 (MQ=255)
                                                 ttCAGAGAAACGCCACCTAAAACGCAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAAtttt                                                                                                   <  4:356042/149‑1 (MQ=255)
                                                 ttCAGAGAAACGCCACCTAAAACGCAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAAtttt                                                                                                   <  2:606043/149‑1 (MQ=255)
                                                 ttCAGAGAAACGCCACCTAAAACGCAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAAtttt                                                                                                   >  4:748772/1‑149 (MQ=255)
                                                    agagAAACGCCACCTAAAACGCAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGt                                                                                                <  3:647490/149‑1 (MQ=255)
                                                    agagAAACGCCACCTAAAACGCAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGt                                                                                                >  1:147568/1‑149 (MQ=255)
                                                    agagAAACGCCACCTAAAACGCAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGt                                                                                                <  3:735548/149‑1 (MQ=255)
                                                     gagaAACGCCACCTAAAACGCAGGCGGAGATAACAATCAGCTCATAACCAATCGTCGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTg                                                                                               <  4:301824/149‑1 (MQ=255)
                                                     gagaAACGCCACCTAAAACGCAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTg                                                                                               <  3:637691/149‑1 (MQ=255)
                                                     gagaAACGCCACCTAAAACGCAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGATCATACGTGACGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTg                                                                                               >  1:434157/1‑149 (MQ=255)
                                                       gaAACGCCACCTAAAACGCAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTcc                                                                                                                                                >  3:775329/1‑98 (MQ=255)
                                                       gaAACGCCACCTAAAACGCAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTcc                                                                                                                                                <  4:775329/98‑1 (MQ=255)
                                                       gaAACGCCACCTAAAACGCAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCg                                                                                             >  4:599082/1‑149 (MQ=255)
                                                            gCCACCTAAAACGCAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGaacaa                                                                                        <  2:147568/149‑1 (MQ=255)
                                                            gCCACCTAAAACGCAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGaacaa                                                                                        >  2:824636/1‑149 (MQ=255)
                                                              cACCTAAAACGCAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAAcc                                                                                      >  2:159866/1‑149 (MQ=255)
                                                                ccTAAAACGCAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCgg                                                                                    <  4:286735/149‑1 (MQ=255)
                                                                  tAAAACGCAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTa                                                                                  <  4:2972/149‑1 (MQ=255)
                                                                         cAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGc                                                                           <  3:748772/149‑1 (MQ=255)
                                                                         cAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGc                                                                           <  3:153850/149‑1 (MQ=255)
                                                                           ggCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCa                                                                         <  4:649233/149‑1 (MQ=255)
                                                                           ggCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCa                                                                         <  1:893122/149‑1 (MQ=255)
                                                                            gCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAg                                                                        >  2:554389/1‑149 (MQ=255)
                                                                            gCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAg                                                                        >  3:702602/1‑149 (MQ=255)
                                                                                  aTAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGAcgcgc                                                                  >  4:274351/1‑149 (MQ=255)
                                                                                   tAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGAcgcgcg                                                                 >  2:926274/1‑149 (MQ=255)
                                                                                    aaCAATCAGCTCATACCCAGTCGACGTCATTGGCTCCATACACGTCATACGTGAAGCCAGAATTAATCCCACCACCTCTGATACCCCGCCAGTGAGAGCAAAAAAAACAAATTTAGTGCATACATCCGGAACACCCCCCAGCcgcgctc                                                                >  3:371478/1‑147 (MQ=255)
                                                                                       aaTCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCct                                                             >  2:123515/1‑149 (MQ=255)
                                                                                       aaTCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAAAAACCGGTACACACGCCAAACGCGCGGCct                                                             >  2:71637/1‑149 (MQ=255)
                                                                                         tCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCctct                                                           >  2:38177/1‑149 (MQ=255)
                                                                                           aGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTc                                                         >  2:784446/1‑149 (MQ=255)
                                                                                            gCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCg                                                        >  1:848928/1‑149 (MQ=255)
                                                                                              tCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGtt                                                      >  3:39009/1‑149 (MQ=255)
                                                                                                  aaCCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAAtttt                                                                                                   <  3:206578/100‑1 (MQ=255)
                                                                                                  aaCCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAAtttt                                                                                                   >  4:206578/1‑100 (MQ=255)
                                                                                                     cAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCa                                               >  1:525312/1‑149 (MQ=255)
                                                                                                      aaTCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCaa                                              >  3:146584/1‑149 (MQ=255)
                                                                                                             gTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAg                                       <  2:848928/149‑1 (MQ=255)
                                                                                                                 ttGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGgtg                                   >  4:384304/1‑149 (MQ=255)
                                                                                                                  tGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTTCACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGgtgt                                  <  1:605319/149‑1 (MQ=255)
                                                                                                                   ggCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGtt                                 <  1:40476/149‑1 (MQ=255)
                                                                                                                            cTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCaaa                        >  4:611667/1‑149 (MQ=255)
                                                                                                                            cTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCaaa                        <  2:913658/149‑1 (MQ=255)
                                                                                                                            cTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCaaa                        <  3:384304/149‑1 (MQ=255)
                                                                                                                               gTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAggt                     <  4:371478/149‑1 (MQ=255)
                                                                                                                               gTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGGTACGACCAAAggt                     <  4:702602/149‑1 (MQ=255)
                                                                                                                                  aTACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAggtggt                  <  1:655805/149‑1 (MQ=255)
                                                                                                                                   tACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGtt                 <  2:561787/149‑1 (MQ=255)
                                                                                                                                      gTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGATAGCAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAg                                       <  4:269602/124‑1 (MQ=255)
                                                                                                                                      gTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAg                                       >  3:269602/1‑124 (MQ=255)
                                                                                                                                      gTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTa                                >  1:737081/1‑131 (MQ=255)
                                                                                                                                      gTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTa                                <  2:737081/131‑1 (MQ=255)
                                                                                                                                      gTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTa              <  2:297034/149‑1 (MQ=255)
                                                                                                                                       tGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCTAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTAt             <  1:926274/149‑1 (MQ=255)
                                                                                                                                       tGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTAt             >  3:677236/1‑149 (MQ=255)
                                                                                                                                        gAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTAtt            <  2:934623/149‑1 (MQ=255)
                                                                                                                                          aGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTca          >  2:301154/1‑149 (MQ=255)
                                                                                                                                           gCCGGAATTATTCTGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTa              <  3:341021/144‑1 (MQ=255)
                                                                                                                                           gCCGGAATTATTCTGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTa              >  4:341021/1‑144 (MQ=255)
                                                                                                                                           gCCGGAATTATTCTGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTc           <  2:99725/147‑1 (MQ=255)
                                                                                                                                           gCCGGAATTATTCTGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTc           >  1:99725/1‑147 (MQ=255)
                                                                                                                                           gCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTcag         <  2:458867/149‑1 (MQ=255)
                                                                                                                                           gCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTcag         >  2:724806/1‑149 (MQ=255)
                                                                                                                                            ccAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTcagc        >  2:684402/1‑149 (MQ=255)
                                                                                                                                            ccAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTcagc        >  3:489665/1‑149 (MQ=255)
                                                                                                                                              aGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTcagcag      >  4:510090/1‑149 (MQ=255)
                                                                                                                                              aGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTcagcag      >  4:344389/1‑149 (MQ=255)
                                                                                                                                              aGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGACAGGGTGTTACGACCAAAGGTGGTTTTATAcagcag      >  3:616071/1‑149 (MQ=255)
                                                                                                                                               gAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGAtaataa                                                                                                       >  2:279671/1‑51 (MQ=255)
                                                                                                                                               gAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGAtaataa                                                                                                       <  1:279671/51‑1 (MQ=255)
                                                                                                                                               gAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAg                                       <  4:512772/115‑1 (MQ=255)
                                                                                                                                               gAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAg                                       >  3:512772/1‑115 (MQ=255)
                                                                                                                                                  taatTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCaaa  <  1:125694/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TAAAATTAAGATACCCGCCACCACATATGAGATTTTTCCGATGCCACCTTTCAGAGAAACGCCACCTAAAACGCAGGCGGAGATAACAATCAGCTCATAACCAATCGACGTCATTGGCTGCCCACTGGTCATACGTGAAGCCAGAATAATTCCGGCTATCGCTGATACCAGGCCTGAGAGAACAAAGATAATAATTTTGGTGCGAACAACCGGTACACCCGCCAGACGCGCGGCCTCTTCGTTCCCGCCAATCGCCAGGGTGTTACGACCAAAGGTGGTTTTATTCAGCAGCAAA  >  NC_000913/1982683‑1982977

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: