Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 558,387 C→A 6.7% T59N (ACC→AAC) ‡ sfmA → FimA homolog, function unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913558,3870CA 6.7% 347.6 / 10.4 113T59N (ACC→AAC) ‡sfmAFimA homolog, function unknown
Reads supporting (aligned to +/- strand):  ref base C (51/54);  new base A (6/2);  total (57/56)
Fisher's exact test for biased strand distribution p-value = 2.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.94e-01

Only 100 of 117 total aligned reads displayed.

CTGGCTGCCGCACTATTCGCCGCTACGGGTAGTTATGCTGCCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCTGATGCAATCAATGATAATTTATTGGCCATTG  >  NC_000913/558239‑558533
                                                                                                                                                    |                                                                                                                                                  
ctggctgCCGCACTATTCGCCGCTACGGGTAGTTATGCTGCCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTAc                                                                                                                                                    <  2:423138/149‑1 (MQ=255)
ctggctgCCGCACTATTCGCCGCTACGGGTAGTTATGCTGCCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTAc                                                                                                                                                    <  1:996518/149‑1 (MQ=255)
         gCACTATTCGCCGCTACGGGTAGTTATGCTGCCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATAtttt                                                                                                                                           >  1:840503/1‑149 (MQ=255)
         gCACTATTCGCCGCTACGGGTAGTTATGCTGCCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATAtttt                                                                                                                                           <  2:77756/149‑1 (MQ=255)
         gCACTATTCGCCGCTACGGGTAGTTATGCTGCCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGGTACGATAttga                                                                                                                                           <  1:1560763/149‑3 (MQ=255)
           aCTATTCGCCGCTACGGGTAGTTATGCTGCCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCa                                                                                                                                         <  1:1378824/149‑1 (MQ=255)
            cTATTCGCCGCTACGGGTAGTTATGCTGCCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCaa                                                                                                                                        <  1:800370/149‑1 (MQ=255)
            cTATTCGCCGCTACGGGTAGTTATGCTGCCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCaa                                                                                                                                        <  2:1237841/149‑1 (MQ=255)
             tATTCGCCGCTACGGGTAGTTATGCTGCCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAAt                                                                                                                                       >  2:1199855/1‑149 (MQ=255)
             tATTCGCCGCTACGGGTAGTTATGCTGCCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAAt                                                                                                                                       >  1:13714/1‑149 (MQ=255)
              aTTCGCCGCTACGGGTAGTTATGCTGCCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATg                                                                                                                                      <  2:456628/149‑1 (MQ=255)
              aTTCGCCGCTACGGGTAGTTATGCTGCCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATg                                                                                                                                      >  1:443442/1‑149 (MQ=255)
               ttCGCCGCTACGGGTAGTTATGCTGCCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGc                                                                                                                                     >  1:1424035/1‑149 (MQ=255)
               ttCGCCGCTACGGGTAGTTATGCTGCCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGc                                                                                                                                     >  1:1537500/1‑149 (MQ=255)
                  gccgcTACGGGTAGTTATGCTGCCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGt                                                                                                                                  <  2:807965/149‑1 (MQ=255)
                  gccgcTACGGGTAGTTATGCTGCCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGt                                                                                                                                  <  1:314182/149‑1 (MQ=255)
                  gccgcTACGGGTAGTTATGCTGCCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGt                                                                                                                                  <  1:1120162/149‑1 (MQ=255)
                  gccgcTACGGGTAGTTATGCTGCCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGt                                                                                                                                  >  1:844502/1‑149 (MQ=255)
                   ccgcTACGGGTAGTTATGCTGCCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGtt                                                                                                                                 >  1:37348/1‑149 (MQ=255)
                   ccgcTACGGGTAGTTATGCTGCCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGtt                                                                                                                                 >  2:356643/1‑149 (MQ=255)
                    cgcTACGGGTAGTTATGCTGCCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGGTGTCACACTCGGGCAATATCGTACCGATATTTTCAATGCTGTTg                                                                                                                                <  1:668521/149‑1 (MQ=255)
                     gcTACGGGTAGTTATGCTGCCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTgg                                                                                                                               <  2:727123/149‑1 (MQ=255)
                     gcTACGGGTAGTTATGCTGCCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTgg                                                                                                                               <  1:347106/149‑1 (MQ=255)
                     gcTACGGGTAGTTATGCTGCCGTTGTAGATGGCGGTACAATTCACTATGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTgg                                                                                                                               <  2:512884/149‑1 (MQ=255)
                     gcTACGGGTAGTAATGCTGCCGTTGTAGATGGCTGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCATGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACGCGGTCAATATCGTACCGATATTTTCAATGCTGTTgg                                                                                                                               <  1:993887/149‑1 (MQ=255)
                      cTACGGGTAGTTATGCTGCCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGt                                                                                                                              >  1:765135/1‑149 (MQ=255)
                      cTACGGGTAGTTATGCTGCCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGt                                                                                                                              >  1:587857/1‑149 (MQ=255)
                       tACGGGTAGTTATGCTGCCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTa                                                                                                                             <  1:1024780/149‑1 (MQ=255)
                            gTAGTTATGCTGCCGTTGTATATGGAGGTACAATTCAATATGAAGGAGAAATGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACACGGTCAATATCGTACCTATATTTTAAATGCTGTTGGTAATAcc                                                                                                                        >  2:236815/1‑149 (MQ=255)
                             tAGTTATGCTGCCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACct                                                                                                                       <  2:693986/149‑1 (MQ=255)
                               gTTATGCTGCCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGTTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACctct                                                                                                                     <  1:90263/149‑1 (MQ=255)
                                  aTGCTGCCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCa                                                                                                                  <  1:213974/149‑1 (MQ=255)
                                   tgctgcCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCAt                                                                                                                 <  2:1481291/149‑1 (MQ=255)
                                      tgcCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGat                                                                                                                                                >  1:1326870/1‑115 (MQ=255)
                                        ccGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAAtt                                                                                                            <  1:936215/149‑1 (MQ=255)
                                           ttGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTcca                                                                                                         >  1:1488109/1‑149 (MQ=255)
                                            tGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTccat                                                                                                        <  1:405419/149‑1 (MQ=255)
                                            tGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTccat                                                                                                        <  1:1107557/149‑1 (MQ=255)
                                                gATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTccattcac                                                                                                    <  1:1454442/149‑1 (MQ=255)
                                                   ggCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTCGGTAATACCTCTGCATTAATTccattcaccat                                                                                                 <  2:840503/149‑1 (MQ=255)
                                                    gCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTccattcaccatt                                                                                                >  1:536680/1‑149 (MQ=255)
                                                     cGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTccattcaccattc                                                                                               >  2:1044568/1‑149 (MQ=255)
                                                      ggTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTccattcaccattca                                                                                              >  2:1084689/1‑149 (MQ=255)
                                                       gTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAg                                                                                             <  1:1336333/149‑1 (MQ=255)
                                                         aCAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGtt                                                                                           <  2:13714/149‑1 (MQ=255)
                                                           aaTTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGa                                                                                         >  2:732680/1‑149 (MQ=255)
                                                            aTTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTGACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGaa                                                                                        <  1:100180/149‑1 (MQ=255)
                                                            aTTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCTTTAATTCCATTCACCATTCAGTTGaa                                                                                        <  2:925252/149‑1 (MQ=255)
                                                               cACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGa                                                                                     <  1:1458680/149‑1 (MQ=255)
                                                                aCTGTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCAGACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTATTACCTCTGCATTAATTCCATTCCCCATTCAGTTGAACGAc                                                                                    >  2:141950/1‑149 (MQ=255)
                                                                  tttGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTAc                                                                                                                                                    >  1:1028287/1‑83 (MQ=255)
                                                                  tttGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTAc                                                                                                                                                    <  2:1028287/83‑1 (MQ=255)
                                                                  tttGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAactactt                                                                                  >  1:17593/1‑148 (MQ=255)
                                                                  tttGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTg                                                                                  <  1:1502335/149‑1 (MQ=255)
                                                                  tttGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTg                                                                                  >  1:28234/1‑149 (MQ=255)
                                                                           gAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCtgt                                                                         <  1:898815/149‑1 (MQ=255)
                                                                             aCTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCtgttg                                                                       <  1:893340/149‑1 (MQ=255)
                                                                                 gTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGc                                                                   >  2:1394676/1‑149 (MQ=255)
                                                                                    aaTGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGc                                                                >  2:357231/1‑149 (MQ=255)
                                                                                    aaTGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGc                                                                >  2:1169528/1‑149 (MQ=255)
                                                                                    aaTGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGc                                                                <  1:356643/149‑1 (MQ=255)
                                                                                         tgcCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAAtg                                                           >  1:2479/1‑149 (MQ=255)
                                                                                            cTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAAtgctg                                                        <  1:88851/149‑1 (MQ=255)
                                                                                            cTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTAACGATATTTTCACTGCTGTTGGTTATACCTCTGAATTAATTCCATTCACCATTCTGTTGAACGACTTCGATCCTGTTGTTGCCGCTAAtgctc                                                        >  1:66536/1‑148 (MQ=255)
                                                                                             tGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAAtgctgc                                                       <  2:596004/149‑1 (MQ=255)
                                                                                               ttCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCg                                                     >  1:167489/1‑149 (MQ=255)
                                                                                                  aGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTg                                                  >  2:200053/1‑149 (MQ=255)
                                                                                                       aTACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCAttt                                             >  2:345775/1‑149 (MQ=255)
                                                                                                           tGACTCGGCAGACCAGGTTGTCACACTCGGGCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTg                                         <  2:368902/149‑1 (MQ=255)
                                                                                                             aCTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGt                                       <  2:1169824/149‑1 (MQ=255)
                                                                                                             aCTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACATCTGCATTAATTACATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGt                                       >  1:489228/1‑149 (MQ=255)
                                                                                                               tCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCa                                     >  1:1552446/1‑149 (MQ=255)
                                                                                                               tCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCa                                     >  2:487605/1‑149 (MQ=255)
                                                                                                                 ggCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCAgg                                   >  1:405752/1‑149 (MQ=255)
                                                                                                                 ggCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACATCTGCATTAATTCCATTCACCATTCAGTTGAACGAATGCGATCCTGTTGTTGCCGCTAATGCTGCCTTTGCATTTTATGGTCAgg                                   >  2:1303639/1‑149 (MQ=255)
                                                                                                                   cAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCt                                 >  1:784188/1‑149 (MQ=255)
                                                                                                                     gACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATAAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCTGa                               <  2:693177/149‑1 (MQ=255)
                                                                                                                        cAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCTGATGc                            >  2:1600756/1‑149 (MQ=255)
                                                                                                                        cAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCTGATGc                            >  1:821150/1‑149 (MQ=255)
                                                                                                                         aGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCTGATGca                           >  1:691390/1‑149 (MQ=255)
                                                                                                                           gTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACTACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCTGATGcaat                         <  2:37348/149‑1 (MQ=255)
                                                                                                                           gTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCTGATGcaat                         >  2:1208604/1‑149 (MQ=255)
                                                                                                                              gTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCTGATGcaatcaa                      <  2:765135/149‑1 (MQ=255)
                                                                                                                                cacaCTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCTGATGCAATCAATg                    >  1:656143/1‑149 (MQ=255)
                                                                                                                                   acTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCTGATGCAATCAATGATa                 >  1:56064/1‑149 (MQ=255)
                                                                                                                                     tCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCTGATGCAATCAATGATAAt               >  1:805946/1‑149 (MQ=255)
                                                                                                                                      ttgaaaaTATCGGTACGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCTGATGCAATCAATGATAAt               <  1:884440/143‑1 (MQ=255)
                                                                                                                                      ttgaaaaTATCGGTACGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCTGATGCAATCAATGATAAt               >  2:884440/6‑148 (MQ=255)
                                                                                                                                       ggTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTTAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCTGATGCAATCAATGATAAttt             >  1:1577135/1‑149 (MQ=255)
                                                                                                                                        gaaaaTATCGGTACGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCTGATGCAATCAATGATAATTTa            >  2:157088/4‑149 (MQ=255)
                                                                                                                                        gTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCTGATGCAATCAATGATAATTTa            >  2:620461/1‑149 (MQ=255)
                                                                                                                                         tCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCTGAGGCAATCAATGATAATTTAt           >  2:312127/1‑149 (MQ=255)
                                                                                                                                         aaaaTATCGGTACGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCTGATGCAATCAATGATAATTTAt           >  2:1500608/3‑149 (MQ=255)
                                                                                                                                         aaaaTATCGGTACGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCTGATGCAATCAATGATAATTTAt           >  2:1343347/3‑149 (MQ=255)
                                                                                                                                          cAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGACGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCTGATGCAATCAATGATAATTTAtt          <  1:20857/149‑1 (MQ=255)
                                                                                                                                           aaTATCGTACCGATATTTTCAATGCGGGTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGAGCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCTGATGCAATCAATGATAATTTATTg         <  2:239065/149‑1 (MQ=255)
                                                                                                                                            atatCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCTGATGCAATCAATGATAATTTATTgg        <  2:1537500/149‑1 (MQ=255)
                                                                                                                                                 gTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCTGATGCAATCAATGATAATTTATTGGCCAtt   <  1:1315008/149‑1 (MQ=255)
                                                                                                                                                 gTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCTGATGCAATCAATGATAATTTATTGGCCAtt   <  1:1169528/149‑1 (MQ=255)
                                                                                                                                                  tACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCTGATGCAATCAATGATAATTTATTGGCCACTg  <  2:656143/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CTGGCTGCCGCACTATTCGCCGCTACGGGTAGTTATGCTGCCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCTGATGCAATCAATGATAATTTATTGGCCATTG  >  NC_000913/558239‑558533

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: