Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 362,830 C→A C117F (TGT→TTT)  lacY ← lactose permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913362,8300CA97.4% 300.5 / 0.5 77coding (350/1254 nt)lacYlactose permease
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (34/41);  minor base . (0/2);  total (34/43)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

ATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTT  >  NC_000913/362685‑362975
                                                                                                                                                  |                                                                                                                                                 
attGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAaaaaaaaa                                                                                                                                               <  1:1160431/150‑1 (MQ=255)
   gATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGcc                                                                                                                                            <  1:71257/150‑1 (MQ=255)
      ggTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTa                                                                                                                                          <  1:317050/149‑1 (MQ=255)
       gTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGAt                                                                                                                                       <  2:926453/151‑1 (MQ=255)
        tGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGAt                                                                                                                                       >  1:570035/1‑150 (MQ=255)
          aaCATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATaa                                                                                                                                     >  2:382031/1‑150 (MQ=255)
           aCATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATaaa                                                                                                                                    <  2:75703/150‑1 (MQ=255)
            cATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGAAGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATa                                                                                                                                  <  1:1066830/151‑1 (MQ=255)
             atgatgCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATac                                                                                                                                 <  1:157738/151‑1 (MQ=255)
               gatgCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATacc                                                                                                                                <  2:533927/150‑1 (MQ=255)
                atgCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATaccac                                                                                                                              <  1:844417/151‑1 (MQ=255)
                     gACAATAGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGAACAAATTCGAAATTACTGCGACCGCTGGCTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAAGCCTAGATAAATACCACcaac                                                                                                                          <  2:650906/151‑1 (MQ=255)
                        aaTCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATc                                                                                                                      <  2:389327/151‑1 (MQ=255)
                        aaTCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGATGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATc                                                                                                                      <  2:1078055/151‑1 (MQ=255)
                           cGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTTACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCAAGATAAATACCACCAACAATcta                                                                                                                    <  2:776711/150‑3 (MQ=255)
                           cGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGa                                                                                                                    <  1:1004238/150‑1 (MQ=255)
                           cGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGAt                                                                                                                   <  2:989132/151‑1 (MQ=255)
                               gCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCct                                                                                                                >  1:597996/1‑150 (MQ=255)
                                 acacaGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCctact                                                                                                             <  1:847863/151‑1 (MQ=255)
                                  cacaGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCctacta                                                                                                            <  2:62960/151‑1 (MQ=255)
                                  cacaGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCctact                                                                                                             <  2:69511/150‑1 (MQ=255)
                                    caGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTaa                                                                                                           >  2:526741/1‑150 (MQ=255)
                                       cgcCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAAtg                                                                                                       <  2:1185818/151‑1 (MQ=255)
                                             gCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGttgtat                                                                                                 >  1:584185/1‑151 (MQ=255)
                                               cAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGttgtattg                                                                                               >  1:1145207/1‑151 (MQ=255)
                                                aaCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGttgtattgt                                                                                              <  1:92812/151‑1 (MQ=255)
                                                  cacaGCCAAACATCCGCGCGCGACCAAATTCGAAATTAATGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGttgtattgt                                                                                              >  2:162992/1‑149 (MQ=255)
                                                    caGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAAc                                                                                           >  1:209552/1‑150 (MQ=255)
                                                       ccAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTg                                                                                       <  2:717873/151‑1 (MQ=255)
                                                         aaaCATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGc                                                                                     >  2:44320/1‑151 (MQ=255)
                                                                cgcgcgcgACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGa                                                                                 >  2:81833/1‑149 (MQ=255)
                                                                cgcgcgcgACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGa                                                                              <  2:1077299/151‑1 (MQ=255)
                                                                  cgcgcgACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATa                                                                            <  2:774768/151‑1 (MQ=255)
                                                                   gcgcgACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATa                                                                            >  2:568028/1‑150 (MQ=255)
                                                                    cgcgACCAAATTCGAAATTACTGCGACGGCTGACTTTCTAAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATaaa                                                                          <  1:503467/151‑1 (MQ=255)
                                                                     gcgACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATaaaa                                                                         <  2:464590/151‑1 (MQ=255)
                                                                     gcgACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATaaa                                                                          >  1:750568/1‑150 (MQ=255)
                                                                      cgACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATaa                                                                           <  2:60535/148‑1 (MQ=255)
                                                                          cAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATaaa                                                                    <  1:815454/151‑1 (MQ=255)
                                                                          cAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATaa                                                                     >  1:936489/1‑150 (MQ=255)
                                                                          cAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATa                                                                      <  1:557838/149‑1 (MQ=255)
                                                                           aaaTTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATaa                                                                     >  2:395337/1‑149 (MQ=255)
                                                                           aaaTTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAg                                                                   <  1:537529/151‑1 (MQ=255)
                                                                               tCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAAc                                                                <  2:689910/150‑1 (MQ=255)
                                                                                cGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACg                                                               <  1:512123/150‑1 (MQ=255)
                                                                                   aaTTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGgcgc                                                           >  1:583440/1‑151 (MQ=255)
                                                                                   aaTTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGgcgc                                                           <  2:1162180/151‑1 (MQ=255)
                                                                                     ttACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCa                                                          >  1:142023/1‑150 (MQ=255)
                                                                                      tACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCaa                                                         >  1:832462/1‑150 (MQ=255)
                                                                                       aCTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCaaa                                                        >  1:842676/1‑150 (MQ=255)
                                                                                       aCTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAAc                                                       >  1:923674/1‑151 (MQ=255)
                                                                                          gCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACa                                                      <  2:596312/149‑1 (MQ=255)
                                                                                           cGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATc                                                    >  2:758304/1‑150 (MQ=255)
                                                                                           cGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATc                                                    >  2:921917/1‑150 (MQ=255)
                                                                                           cGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCa                                                   <  1:962115/150‑1 (MQ=255)
                                                                                                   gACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATg                                           >  1:915065/1‑151 (MQ=255)
                                                                                                      tttCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCg                                        >  1:673597/1‑151 (MQ=255)
                                                                                                        tctcAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCgg                                       >  2:870511/1‑150 (MQ=255)
                                                                                                        tctcAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCgg                                       <  2:640673/150‑1 (MQ=255)
                                                                                                           cAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGta                                     <  2:840297/149‑1 (MQ=255)
                                                                                                                 aTGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCca                             <  1:125238/151‑1 (MQ=255)
                                                                                                                      tctACTGCCGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCAcagc                         >  1:653419/1‑150 (MQ=255)
                                                                                                                       ctACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCAcagc                         <  1:834526/149‑1 (MQ=255)
                                                                                                                          ctgctgGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAgg                      >  2:97365/1‑149 (MQ=255)
                                                                                                                          ctgctgGCGCACCGGCGTTAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAgg                      <  1:452602/148‑1 (MQ=255)
                                                                                                                           tgctgGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTAt                   >  2:828753/1‑151 (MQ=255)
                                                                                                                            gctgGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTAtt                  >  2:827678/1‑151 (MQ=255)
                                                                                                                               ggCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTgc               <  1:526756/151‑1 (MQ=255)
                                                                                                                               ggCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTg                <  1:946146/150‑1 (MQ=255)
                                                                                                                                gcgcACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTgcg              <  2:1020939/151‑1 (MQ=255)
                                                                                                                                  gcACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCa            <  1:915214/151‑1 (MQ=255)
                                                                                                                                   cACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTgcgc             >  1:643092/1‑149 (MQ=255)
                                                                                                                                   cACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTgcgc             >  1:1087736/1‑149 (MQ=255)
                                                                                                                                     ccGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGcc         >  1:909182/1‑151 (MQ=255)
                                                                                                                                     ccGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGc          <  2:135766/150‑1 (MQ=255)
                                                                                                                                      cGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGccc        >  2:831848/1‑151 (MQ=255)
                                                                                                                                           ttAAAAAAAAAGCCTAGATAAATACCACAAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGtt   >  2:449160/1‑151 (MQ=255)
                                                                                                                                            tAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGttt  >  2:18195/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
ATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTT  >  NC_000913/362685‑362975

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