Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 562,582 A→G T83A (ACT→GCT)  sfmH → FimA homolog, function unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913562,5820AG100.0% 305.2 / NA 79T83A (ACT→GCT) sfmHFimA homolog, function unknown
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (30/49);  total (30/49)

ATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACGTTTCGCAAAACCTTCCATTCG  >  NC_000913/562435‑562729
                                                                                                                                                   |                                                                                                                                                   
atCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCt                                                                                                                                                   >  2:539443/1‑150 (MQ=255)
 tCTTTAATAGCACCAATAACCAGCCAGGGCGGGTTGTTGTTTTACCTGAAAAATCCGGCTGGGTAGGTGTCTCGGAAATTTGTCCTCCCGGTACGCTGGTGAATTATACATACAGTAGTTATGTGACCAACTTTATTGTTCAGGTAGCTAt                                                                                                                                                 >  2:846395/1‑151 (MQ=255)
 tCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTAt                                                                                                                                                 >  2:1005238/1‑151 (MQ=255)
    ttAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATc                                                                                                                                                <  2:511717/149‑1 (MQ=255)
    ttAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATc                                                                                                                                                <  1:809594/149‑1 (MQ=255)
     tAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGa                                                                                                                                              >  2:660258/1‑150 (MQ=255)
      aaTAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGAt                                                                                                                                             <  1:796090/150‑1 (MQ=255)
        tAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAAt                                                                                                                                          >  2:126373/1‑151 (MQ=255)
         aGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAAtt                                                                                                                                         >  1:250904/1‑151 (MQ=255)
         aGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAAt                                                                                                                                          <  2:987066/150‑1 (MQ=255)
            aCCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATtat                                                                                                                                       >  1:798020/1‑150 (MQ=255)
              cAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATaaa                                                                                                                                    <  2:872705/151‑1 (MQ=255)
               aaTAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAat                                                                                                                                   <  2:238912/151‑1 (MQ=255)
                aTAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAat                                                                                                                                   >  2:97202/1‑150 (MQ=255)
                  aaCCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAatat                                                                                                                                 <  1:1025526/150‑1 (MQ=255)
                  aaCCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAatat                                                                                                                                 >  2:34516/1‑150 (MQ=255)
                   aCCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAatatat                                                                                                                               <  2:510928/151‑1 (MQ=255)
                   aCCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAatatat                                                                                                                               >  2:1046941/1‑151 (MQ=255)
                    ccagccagGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATg                                                                                                                              >  2:480568/1‑151 (MQ=255)
                       gccagGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAAGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCaa                                                                                                                           <  1:413074/151‑1 (MQ=255)
                         cagGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAAtt                                                                                                                         <  2:623320/151‑1 (MQ=255)
                          agGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTa                                                                                                                        <  1:1061951/151‑1 (MQ=255)
                            ggCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACa                                                                                                                      <  1:533902/151‑1 (MQ=255)
                              cAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATg                                                                                                                    >  1:563981/1‑151 (MQ=255)
                               aGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATg                                                                                                                    >  1:128708/1‑150 (MQ=255)
                                   tgttgtTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTa                                                                                                                <  2:706387/150‑1 (MQ=255)
                                   tgttgtTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTa                                                                                                                <  2:500941/150‑1 (MQ=255)
                                       gttTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTAt                                                                                                           <  2:340332/151‑1 (MQ=255)
                                          ttACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTAtt                                                                                                          >  1:690314/1‑149 (MQ=255)
                                           tACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAg                                                                                                        <  1:45391/150‑1 (MQ=255)
                                                 aaaaaTCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCTCCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGAt                                                                                                                  <  1:609028/134‑1 (MQ=255)
                                                 aaaaaTCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGAt                                                                                                                  >  2:609011/1‑134 (MQ=255)
                                                  aaaaTCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATg                                                                                                <  2:342819/151‑1 (MQ=255)
                                                   aaaTCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGa                                                                                               <  2:261775/151‑1 (MQ=255)
                                                        cGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTg                                                                                          <  1:666776/151‑1 (MQ=255)
                                                             gggTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGa                                                                                     <  2:804691/151‑1 (MQ=255)
                                                              ggTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGAt                                                                                    <  2:595756/151‑1 (MQ=255)
                                                                tAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGAt                                                                                    >  2:880235/1‑149 (MQ=255)
                                                                 aGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAgt                                                                                 <  1:80532/151‑1 (MQ=255)
                                                                 aGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAgt                                                                                 >  2:935291/1‑151 (MQ=255)
                                                                 aGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAg                                                                                  <  2:237639/150‑1 (MQ=255)
                                                                         cAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATgg                                                                          <  1:802524/150‑1 (MQ=255)
                                                                          aGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGat                                                                        >  2:1032424/1‑151 (MQ=255)
                                                                          aGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGat                                                                        <  1:937095/151‑1 (MQ=255)
                                                                            cAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGatat                                                                      <  2:120/151‑1 (MQ=255)
                                                                            cAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGatat                                                                      <  1:177661/151‑1 (MQ=255)
                                                                             aaTTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATAtt                                                                     <  1:815436/151‑1 (MQ=255)
                                                                                 tGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAg                                                                  <  2:69872/150‑1 (MQ=255)
                                                                                    ccACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTc                                                               <  1:97204/150‑1 (MQ=255)
                                                                                         cGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCg                                                         <  1:788746/151‑1 (MQ=255)
                                                                                          ggTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGc                                                        <  1:349038/151‑1 (MQ=255)
                                                                                           gTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCa                                                       <  2:525543/151‑1 (MQ=255)
                                                                                               gCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAAt                                                   <  2:190451/151‑1 (MQ=255)
                                                                                                  ggTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAAtt                                                  <  1:14930/149‑1 (MQ=255)
                                                                                                      aaTTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTc                                            >  2:1027879/1‑151 (MQ=255)
                                                                                                          atacatacCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGa                                         <  1:260379/150‑1 (MQ=255)
                                                                                                             catacCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATggg                                     >  2:489654/1‑151 (MQ=255)
                                                                                                              atacCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATggg                                     >  2:6093/1‑150 (MQ=255)
                                                                                                                acCGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATAAGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAAc                                  >  1:405042/1‑151 (MQ=255)
                                                                                                                  cGTAGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAg                                <  2:5137/151‑1 (MQ=255)
                                                                                                                     aGTTATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATc                             <  2:390580/151‑1 (MQ=255)
                                                                                                                       ttATGTCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCt                           >  1:677533/1‑151 (MQ=255)
                                                                                                                            tCACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACGt                      >  1:464005/1‑151 (MQ=255)
                                                                                                                             cACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACGtt                     <  2:509300/151‑1 (MQ=255)
                                                                                                                             cACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACGtt                     <  2:1001335/151‑1 (MQ=255)
                                                                                                                             cACCAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACGtt                     <  1:603692/151‑1 (MQ=255)
                                                                                                                              aCCAACTTTATTGTTCAGGCAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACGttt                    >  2:341122/1‑151 (MQ=255)
                                                                                                                               ccAACTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACGttt                    <  1:132142/150‑1 (MQ=255)
                                                                                                                                 aaCTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACGTTTCg                  >  1:307781/1‑150 (MQ=255)
                                                                                                                                 aaCTTTATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACGTTTCGc                 >  2:408058/1‑151 (MQ=255)
                                                                                                                                     ttATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACGTTTCGCaa               <  1:250850/149‑1 (MQ=255)
                                                                                                                                      tATTGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACGTTTCGCaaaa             <  1:24998/150‑1 (MQ=255)
                                                                                                                                         tGTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACGTTTCGCAAAACCt          <  2:419721/150‑1 (MQ=255)
                                                                                                                                          gTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACGTTTCGCAAAACCtt         <  1:1011167/150‑1 (MQ=255)
                                                                                                                                          gTTCAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACGTTTCGCAAAACCTTc        >  2:284288/1‑151 (MQ=255)
                                                                                                                                             cAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACGTTTCGCAAAACCTTCCAt     <  2:507523/151‑1 (MQ=255)
                                                                                                                                             cAGGAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACGTTTCGCAAAACCTTCCAt     >  1:957371/1‑151 (MQ=255)
                                                                                                                                               ggAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACGTTTCGCAAAACCTTCCATTc   <  2:641523/151‑1 (MQ=255)
                                                                                                                                                gAAGCTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACGTTTCGCAAAACCTTCCATTCg  >  1:910010/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGAAACTATCGATAATTATAAATATATGCAATTACATGATTATCTATTAGGTGCGATGAGTCTGGTTGATAGTGTGATGGATATTCAGTTCCCCCCGCAAAATTATATTCGGATGGGAACAGATCCTAACGTTTCGCAAAACCTTCCATTCG  >  NC_000913/562435‑562729

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