Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 558,385 T→G 5.9% R58R (CGT→CGG sfmA → FimA homolog, function unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913558,3850TG 5.9% 375.3 / 9.7 118R58R (CGT→CGGsfmAFimA homolog, function unknown
Reads supporting (aligned to +/- strand):  ref base T (54/57);  new base G (5/2);  total (59/59)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

Only 100 of 120 total aligned reads displayed.

CGCTGGCTGCCGCACTATTCGCCGCTACGGGTAGTTATGCTGCCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCTGATGCAATCAATGATAATTTATTGGCCATTG  >  NC_000913/558237‑558533
                                                                                                                                                    |                                                                                                                                                    
cgcTGGCTGCCGCACTATTCGCCGCTACGGGTAGTTATGCTGCCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGt                                                                                                                                                      >  2:1386431/1‑149 (MQ=255)
 gctggctgCCGCACTATTCGCCGCTACGGGTAGTTATGCTGCCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTa                                                                                                                                                     >  1:1169824/1‑149 (MQ=255)
  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)
                aTTCGCCGCTACGGGTAGTTATGCTGCCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATg                                                                                                                                      >  1:443442/1‑149 (MQ=255)
                 ttCGCCGCTACGGGTAGTTATGCTGCCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGc                                                                                                                                     >  1:1424035/1‑149 (MQ=255)
                 ttCGCCGCTACGGGTAGTTATGCTGCCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGc                                                                                                                                     >  1:1537500/1‑149 (MQ=255)
                 ttCGCCGCTACGGGGAGTTATGCTGCCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGc                                                                                                                                     <  1:766007/149‑1 (MQ=255)
                    gccgcTACGGGTAGTTATGCTGCCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGt                                                                                                                                  >  1:844502/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)
                     ccgcTACGGGTAGTTATGCTGCCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGtt                                                                                                                                 >  1:37348/1‑149 (MQ=255)
                       gcTACGGGTAGTTATGCTGCCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTgg                                                                                                                               <  1:347106/149‑1 (MQ=255)
                       gcTACGGGTAGTTATGCTGCCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTgg                                                                                                                               <  1:872771/149‑1 (MQ=255)
                       gcTACGGGTAGTTATGCTGCCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTgg                                                                                                                               <  2:28220/149‑1 (MQ=255)
                       gcTACGGGTAGTTATGCTGCCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTgg                                                                                                                               <  2:727123/149‑1 (MQ=255)
                       gcTACGGGTAGTTATGCTGCCGTTGTAGATGGCGGTACAATTCACTATGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTgg                                                                                                                               <  2:512884/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)
                                        tgcCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGat                                                                                                                                                >  1:1326870/1‑115 (MQ=255)
                                          ccGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAAtt                                                                                                            <  1:936215/149‑1 (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)
                                                       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)
                                                                    tttGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTg                                                                                  >  2:630509/1‑149 (MQ=255)
                                                                             taaCTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCtgt                                                                         <  1:904940/148‑1 (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                                                                <  1:356643/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)
                                                                                                      tGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCa                                                >  2:1145643/1‑149 (MQ=255)
                                                                                                       ggaTACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCAt                                               >  1:1290002/3‑149 (MQ=255)
                                                                                                        aaTACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATCCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCAtt                                              >  1:568309/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                                     >  2:487605/1‑149 (MQ=255)
                                                                                                                 tCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCa                                     >  1:1552446/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)
                                                                                                                        aCCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCCATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCTGAt                              >  1:1555580/1‑149 (MQ=255)
                                                                                                                          cAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCTGATGc                            >  1:821150/1‑149 (MQ=255)
                                                                                                                          cAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCTGATGc                            >  2:1600756/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)
                                                                                                                                gTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTATTTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCTGATGcaatcaa                      <  1:236815/149‑1 (MQ=255)
                                                                                                                                gTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCTGATGcaatcaa                      <  2:765135/149‑1 (MQ=255)
                                                                                                                                  cacaCTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCTGATGCAATCAATg                    >  1:656143/1‑149 (MQ=255)
                                                                                                                                     acTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCTGATGCAATCAATGATa                 >  1:56064/1‑149 (MQ=255)
                                                                                                                                        ttgaaaaTATCGGTACGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCTGATGCAATCAATGATAAt               >  2:884440/6‑148 (MQ=255)
                                                                                                                                        ttgaaaaTATCGGTACGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCTGATGCAATCAATGATAAt               <  1:884440/143‑1 (MQ=255)
                                                                                                                                          gaaaaTATCGGTACGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCTGATGCAATCAATGATAATTTa            >  2:157088/4‑149 (MQ=255)
                                                                                                                                          gTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCTGATGCAATCAATGATAATTTa            >  2:620461/1‑149 (MQ=255)
                                                                                                                                          gTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCTGATGCAATCAATGATAATTTa            >  1:391763/1‑149 (MQ=255)
                                                                                                                                           tCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCTGAGGCAATCAATGATAATTTAt           >  2:312127/1‑149 (MQ=255)
                                                                                                                                           aaaaTATCGGTACGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCTGATGCAATCAATGATAATTTAt           >  2:1501719/3‑149 (MQ=255)
                                                                                                                                           aaaaTATCGGTACGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCTGATGCAATCAATGATAATTTAt           >  2:1343347/3‑149 (MQ=255)
                                                                                                                                           aaaaTATCGGTACGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCTGATGCAATCAATGATAATTTAt           >  2:1500608/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)
                                                                                                                                                atCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCAGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCTGATGCAATCAATGATAATTTATTGGcc      >  1:1450107/1‑149 (MQ=255)
                                                                                                                                                 tCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCTGATGCAATCAATGATAATTTATTGGCCa     <  2:405752/149‑1 (MQ=255)
                                                                                                                                                   gTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCTGATGCAATCAATGATAATTTATTGGCCAtt   <  1:1315008/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGCTGGCTGCCGCACTATTCGCCGCTACGGGTAGTTATGCTGCCGTTGTAGATGGCGGTACAATTCACTTTGAAGGCGAACTGGTGAATGCTGCCTGTTCAGTGAATACTGACTCGGCAGACCAGGTTGTCACACTCGGTCAATATCGTACCGATATTTTCAATGCTGTTGGTAATACCTCTGCATTAATTCCATTCACCATTCAGTTGAACGACTGCGATCCTGTTGTTGCCGCTAATGCTGCCGTTGCATTTTCTGGTCAGGCTGATGCAATCAATGATAATTTATTGGCCATTG  >  NC_000913/558237‑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: