Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,535,964 A→T 8.9% E44V (GAA→GTA)  crr → glucose‑specific enzyme IIA component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,535,9640AT 8.9% 672.2 / 41.1 203E44V (GAA→GTA) crrglucose‑specific enzyme IIA component of PTS
Reads supporting (aligned to +/- strand):  ref base A (89/95);  new base T (6/12);  total (96/107)
Fisher's exact test for biased strand distribution p-value = 3.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

Only 100 of 207 total aligned reads displayed.

TGCTTAGGAGAAGATCATGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGG  >  NC_000913/2535818‑2536111
                                                                                                                                                  |                                                                                                                                                   
tGCTTAGGAGAAGATCATGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGtaa                                                                                                                                                   >  3:388593/1‑149 (MQ=255)
   ttAGGAGAAGATCATGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGaaaaaa                                                                                                                                                <  4:75793/149‑1 (MQ=255)
       gagaAGATCATGAGTTTGTTCGATAAACTGAAATCTCTCGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGACGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGTAAAAATCGt                                                                                                                                            <  2:168851/149‑1 (MQ=255)
        agaagaTCATGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATAGCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGtt                                                                                                                                           <  4:857301/149‑1 (MQ=255)
         gaagaTCATGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTg                                                                                                                                          <  3:503904/149‑1 (MQ=255)
          aagaTCATGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTgg                                                                                                                                         <  1:687938/149‑1 (MQ=255)
             aTCATGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGa                                                                                                                                      >  1:145816/1‑149 (MQ=255)
              tCATGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGAt                                                                                                                                     >  1:60077/1‑149 (MQ=255)
              tCATGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGAt                                                                                                                                     <  2:1373646/149‑1 (MQ=255)
               cATGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATg                                                                                                                                    <  4:1425179/149‑1 (MQ=255)
                aTGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATgg                                                                                                                                   >  4:479794/1‑149 (MQ=255)
                aTGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATg                                                                                                                                    <  3:330647/148‑1 (MQ=255)
                  gggTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGTAAAAATCGtt                                                                                                                                           >  4:1315260/1‑139 (MQ=255)
                  gggTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTa                                                                                                                                 >  3:336335/1‑149 (MQ=255)
                   ggTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTAt                                                                                                                                <  1:1111283/149‑1 (MQ=255)
                   ggTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGCGATGGTAt                                                                                                                                <  2:229246/149‑1 (MQ=255)
                    gTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTAtt                                                                                                                               >  2:936665/1‑149 (MQ=255)
                     tttGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTg                                                                                                                              >  3:1175686/1‑149 (MQ=255)
                     tttGTTAGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTg                                                                                                                              >  4:660327/1‑149 (MQ=255)
                        gTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTa                                                                                                                           >  3:1082969/1‑149 (MQ=255)
                         ttCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTAt                                                                                                                          <  4:422676/149‑1 (MQ=255)
                          tCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATc                                                                                                                         <  1:1100998/149‑1 (MQ=255)
                                aaCTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCa                                                                                                                   >  3:1451696/1‑149 (MQ=255)
                                  cTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAAc                                                                                                                 >  3:923876/1‑149 (MQ=255)
                                   tGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGtt                                                                                                                                           >  2:1402065/1‑122 (MQ=255)
                                       aTCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTa                                                                                                            >  4:288440/1‑149 (MQ=255)
                                       aTCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTa                                                                                                            <  2:526052/149‑1 (MQ=255)
                                        tctcTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTaa                                                                                                           >  2:1011650/1‑149 (MQ=255)
                                            tGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACaaa                                                                                                       >  3:939273/1‑149 (MQ=255)
                                             ggTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGTAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACaaaa                                                                                                      <  1:250947/149‑1 (MQ=255)
                                              gTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGTAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAAt                                                                                                     <  2:1196123/149‑1 (MQ=255)
                                              gTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAAt                                                                                                     <  2:937113/149‑1 (MQ=255)
                                              gTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAAt                                                                                                     <  4:698043/149‑1 (MQ=255)
                                               tttCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATg                                                                                                    >  4:308734/1‑149 (MQ=255)
                                                  ccGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTc                                                                                                 <  2:676814/149‑1 (MQ=255)
                                                    gacgacAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTcgc                                                                                               >  4:614785/1‑149 (MQ=255)
                                                     acgacAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTcgcg                                                                                              <  2:846111/149‑1 (MQ=255)
                                                          aagaagGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACGAACGGGTAACAAAATGGTCGCTCGAGt                                                                                         >  4:1281610/1‑149 (MQ=255)
                                                          aagaagGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAAc                                                                                                                 <  4:1233442/125‑1 (MQ=255)
                                                              agGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGAc                                                                                     <  1:473446/149‑1 (MQ=255)
                                                               ggATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACg                                                                                    >  1:52495/1‑149 (MQ=255)
                                                               ggATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACg                                                                                    >  2:74857/1‑149 (MQ=255)
                                                                gATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGTAAAAATCGTTGGTGATGGTATTGc                                                                                                                             <  1:1228578/107‑1 (MQ=255)
                                                                gATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGTAAAAATCGTTGGTGATGGTATTGc                                                                                                                             >  2:1228578/1‑107 (MQ=255)
                                                                gATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACgg                                                                                   <  1:810460/149‑1 (MQ=255)
                                                                 aTACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGc                                                                                  <  3:795612/149‑1 (MQ=255)
                                                                    ctggAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACTTGCCGGATGTCGTTTTTGCGGCAAAAATCGTTGCTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCAcc                                                                               >  4:551845/3‑149 (MQ=255)
                                                                       gAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCTCGCCAGTAGACGGCACCAtt                                                                            >  2:1169861/1‑149 (MQ=255)
                                                                        aaCTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTg                                                                           <  1:804170/149‑1 (MQ=255)
                                                                          cTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGt                                                                         <  3:200714/149‑1 (MQ=255)
                                                                          cTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGt                                                                         <  1:1292890/149‑1 (MQ=255)
                                                                                  aTCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACgg                                                                                   >  3:521615/1‑131 (MQ=255)
                                                                                  aTCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACgg                                                                                   <  4:521615/131‑1 (MQ=255)
                                                                                  aTCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCtt                                                                 <  3:697563/149‑1 (MQ=255)
                                                                                  aTCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCtt                                                                 >  4:1126805/1‑149 (MQ=255)
                                                                                  aTCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCtt                                                                 <  4:793192/149‑1 (MQ=255)
                                                                                   tCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCttt                                                                <  3:288440/149‑1 (MQ=255)
                                                                                   tCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCttt                                                                <  1:1011783/149‑1 (MQ=255)
                                                                                   tCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCttt                                                                <  1:943892/149‑1 (MQ=255)
                                                                                   tCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCttt                                                                <  2:1092858/149‑1 (MQ=255)
                                                                                   tCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCttt                                                                <  1:207072/149‑1 (MQ=255)
                                                                                     aTTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGa                                                              >  1:1010452/1‑149 (MQ=255)
                                                                                       tGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGaaa                                                            <  1:1483296/149‑1 (MQ=255)
                                                                                          tCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAaacca                                                         >  2:1512579/1‑149 (MQ=255)
                                                                                           ccGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGGTGGTGATTGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGTCAGTAGAGGGCACCATTGGTAAAAGCTTTGAaaccaa                                                        >  4:368824/1‑149 (MQ=255)
                                                                                             gCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAaaccaacc                                                      >  1:666367/1‑149 (MQ=255)
                                                                                                          aTCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGTAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTAt                                         >  2:1269487/1‑149 (MQ=255)
                                                                                                          aTCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTAt                                         <  4:1006969/149‑1 (MQ=255)
                                                                                                           tCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGTAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCttt                                                                >  1:902618/1‑125 (MQ=255)
                                                                                                           tCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATc                                        <  2:1141846/149‑1 (MQ=255)
                                                                                                             gTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTg                                                               <  4:172961/124‑1 (MQ=255)
                                                                                                             gTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTg                                                               >  3:172961/1‑124 (MQ=255)
                                                                                                             gTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGa                                      >  1:1116659/1‑149 (MQ=255)
                                                                                                             gTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGa                                      >  2:1323699/1‑149 (MQ=255)
                                                                                                              tCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGaa                                     >  1:137633/1‑149 (MQ=255)
                                                                                                                 atatCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCt                                  >  2:1213084/1‑149 (MQ=255)
                                                                                                                     cGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATa                              <  1:1069929/149‑1 (MQ=255)
                                                                                                                     cGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATa                              <  3:1016033/149‑1 (MQ=255)
                                                                                                                     cGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATa                              <  2:1429356/149‑1 (MQ=255)
                                                                                                                     cGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATa                              <  2:479638/149‑1 (MQ=255)
                                                                                                                      gAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATAGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAg                             >  1:741369/1‑149 (MQ=255)
                                                                                                                             tGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGtt                      <  2:999239/149‑1 (MQ=255)
                                                                                                                             tGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGtt                      <  2:4967/149‑1 (MQ=255)
                                                                                                                                cGGATGTCGTTTTTGCGGTAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGaa                   <  1:1269487/149‑1 (MQ=255)
                                                                                                                                cGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGaa                   <  3:122248/149‑1 (MQ=255)
                                                                                                                                 ggATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAAc                  >  4:690709/1‑149 (MQ=255)
                                                                                                                                    tGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGt               <  4:923876/149‑1 (MQ=255)
                                                                                                                                     gTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGtt              <  3:556110/149‑1 (MQ=255)
                                                                                                                                     gTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGtt              >  1:1136352/1‑149 (MQ=255)
                                                                                                                                      tCGTTTTTGCGGTAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTc             <  2:1495700/149‑1 (MQ=255)
                                                                                                                                      tCGTTTTTGCGGAGAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTc             <  1:1422328/149‑1 (MQ=255)
                                                                                                                                      tCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCtt                                                                 <  1:552865/97‑1 (MQ=255)
                                                                                                                                      tCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCtt                                                                 >  2:552865/1‑97 (MQ=255)
                                                                                                                                      tCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCAtt                                               >  4:392166/1‑115 (MQ=255)
                                                                                                                                      tCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTc             <  3:408265/149‑1 (MQ=255)
                                                                                                                                       cGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCg            >  2:597507/1‑149 (MQ=255)
                                                                                                                                         tttttGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTc          >  3:1348751/1‑149 (MQ=255)
                                                                                                                                                ggAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATGGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCg   >  3:1289318/1‑149 (MQ=255)
                                                                                                                                                ggAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCg   <  3:649391/149‑1 (MQ=255)
                                                                                                                                                 gAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACAGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCgg  <  4:1175686/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TGCTTAGGAGAAGATCATGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGG  >  NC_000913/2535818‑2536111

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