Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 641,355 G→A 6.0% intergenic (+37/+84) ahpF → / ← uspG alkyl hydroperoxide reductase, F52a subunit, FAD/NAD(P)‑binding/universal stress protein UP12

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913641,3550GA 6.0% 414.8 / 11.2 133intergenic (+37/+84)ahpF/uspGalkyl hydroperoxide reductase, F52a subunit, FAD/NAD(P)‑binding/universal stress protein UP12
Reads supporting (aligned to +/- strand):  ref base G (68/56);  new base A (5/3);  total (73/60)
Fisher's exact test for biased strand distribution p-value = 7.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.50e-01

Only 100 of 140 total aligned reads displayed.

GTGACTGTACGACGGTTCCGTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGAACCTAA  >  NC_000913/641212‑641501
                                                                                                                                               |                                                                                                                                                  
gTGACTGTACGACGGTTCCGTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGt                                                                                                                                                >  3:324010/1‑148 (MQ=255)
     tGTACGACGGTTCCGTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAAcc                                                                                                                                          >  5:315090/1‑149 (MQ=255)
       tACGACGGTTCCGTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTg                                                                                                                                        >  7:349060/1‑149 (MQ=255)
       tACGACGGTTCCGTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTg                                                                                                                                        >  4:63468/1‑149 (MQ=255)
                   gTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTa                                                                                                                            <  1:193663/149‑1 (MQ=255)
                   gTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTa                                                                                                                            <  8:158689/149‑1 (MQ=255)
                   gTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTa                                                                                                                            <  1:278579/149‑1 (MQ=255)
                   gTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTa                                                                                                                            >  5:216675/1‑149 (MQ=255)
                    tACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTAt                                                                                                                           >  3:282103/1‑149 (MQ=255)
                    tACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTAt                                                                                                                           <  3:63468/149‑1 (MQ=255)
                      cAAGCATATCATCATCGCCACTGGCGAAGGTGCCAAAGGCTATCTGAGTTCTTTTGACTACCTGAATCGCACAAAAAATGGATAAGAAGAAGTAAGATTAACCTGAAATTGAATAGCCGCCTGGGTCAAACCTGGCGGCTTTTTTATgg                                                                                                                         <  3:344341/149‑1 (MQ=255)
                      cAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATg                                                                                                                          <  4:324010/148‑1 (MQ=255)
                         gCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCAt                                                                                                                      <  6:216675/149‑1 (MQ=255)
                          cAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTggc                                                                                                                                      >  2:286906/1‑132 (MQ=255)
                          cAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTggc                                                                                                                                      <  1:286906/132‑1 (MQ=255)
                            tatcatcatcGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTa                                                                                                                    >  7:332710/2‑148 (MQ=255)
                             atcatcatcGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTaa                                                                                                                   >  2:310424/1‑148 (MQ=255)
                              tcatcatcGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTaaaa                                                                                                                 <  6:315090/149‑1 (MQ=255)
                                   atcGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGccc                                                                                                            >  1:374362/1‑149 (MQ=255)
                                    tcGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGcccc                                                                                                           >  2:251762/1‑149 (MQ=255)
                                    tcGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGcccc                                                                                                           >  8:355714/1‑149 (MQ=255)
                                       ccACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCccgcc                                                                                                        <  6:404161/149‑1 (MQ=255)
                                        cACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGACGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCccgccg                                                                                                       <  1:169965/149‑1 (MQ=255)
                                           tGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGa                                                                                                    >  2:279396/1‑149 (MQ=255)
                                             gCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGAt                                                                                                   <  3:70510/148‑1 (MQ=255)
                                                    tGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCgg                                                                                           <  4:296311/149‑1 (MQ=255)
                                                        aaaGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTggc                                                                                                                                      >  6:265856/1‑102 (MQ=255)
                                                        aaaGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTggc                                                                                                                                      <  5:265856/102‑1 (MQ=255)
                                                        aaaGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCa                                                                                                                       <  7:341104/117‑1 (MQ=255)
                                                        aaaGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAg                                                                                       >  1:106578/1‑149 (MQ=255)
                                                        aaaGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAg                                                                                       <  5:314055/149‑1 (MQ=255)
                                                        aaaGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAg                                                                                       >  4:146160/1‑149 (MQ=255)
                                                        aaaGCCGCTCTGAGTGCCTTTGACTACCTGATTCCCACCAAAACTGCATAAGAAGAACTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGGCAAACCTGTCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCGCCGGTCAg                                                                                       <  6:367143/149‑1 (MQ=255)
                                                         aaGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAgg                                                                                      <  8:332710/149‑1 (MQ=255)
                                                          aGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAggg                                                                                     <  7:420640/149‑1 (MQ=255)
                                                          aGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAggg                                                                                     >  1:302883/1‑149 (MQ=255)
                                                             ctctctGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCgg                                                                                  <  2:374362/149‑1 (MQ=255)
                                                             ctctctGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCgg                                                                                  <  8:179188/149‑1 (MQ=255)
                                                               ctctGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGa                                                                                >  4:147259/1‑149 (MQ=255)
                                                                    aGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGGCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATa                                                                              <  1:241363/146‑1 (MQ=255)
                                                                      tGCTTTTGACTACCTGATTCGCACCAAAAATGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTAt                                                                         >  4:316069/1‑149 (MQ=255)
                                                                         ttttGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGa                                                                                <  6:448216/139‑1 (MQ=255)
                                                                         ttttGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGa                                                                                >  5:448216/1‑139 (MQ=255)
                                                                         ttttGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCt                                                                      >  7:135437/1‑149 (MQ=255)
                                                                         ttttGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCt                                                                      >  4:98714/1‑149 (MQ=255)
                                                                                tACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTgg                                                                    >  5:201390/1‑144 (MQ=255)
                                                                                  ccTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTaa                                                             <  4:336491/149‑1 (MQ=255)
                                                                                   cTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAAc                                                            <  1:279396/149‑1 (MQ=255)
                                                                                      aTTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACg                                                           <  4:22671/147‑1 (MQ=255)
                                                                                          gCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAAc                                                     >  1:86935/1‑149 (MQ=255)
                                                                                          gCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAAc                                                     <  6:452518/149‑1 (MQ=255)
                                                                                           cACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAATGCACAAc                                                     >  4:244067/1‑148 (MQ=255)
                                                                                           cACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGAAATTAACGCACAAc                                                     >  4:90683/1‑148 (MQ=255)
                                                                                           cACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAACCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAAc                                                     >  8:414754/1‑148 (MQ=255)
                                                                                            aCCAAAACTGCGTAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAGAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAAc                                                     >  3:332042/1‑147 (MQ=255)
                                                                                            aCCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAAc                                                     >  6:258352/1‑147 (MQ=255)
                                                                                             ccAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAg                                                  <  2:447631/149‑1 (MQ=255)
                                                                                                aaaCTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCAc                                               >  1:401255/1‑149 (MQ=255)
                                                                                                aaaCTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCAc                                               <  8:51923/149‑1 (MQ=255)
                                                                                                    tGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCgg                                            >  2:448860/1‑148 (MQ=255)
                                                                                                    tGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCgg                                            >  2:7904/1‑148 (MQ=255)
                                                                                                       aTAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCgg                                            >  1:442016/1‑145 (MQ=255)
                                                                                                        tAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGAt                                       >  6:341418/1‑149 (MQ=255)
                                                                                                          agaagaagTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGAtt                                      <  1:254648/148‑1 (MQ=255)
                                                                                                            aagaagTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATtggc                                   <  3:98714/149‑1 (MQ=255)
                                                                                                                agTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATtggcgtg                                <  2:442016/148‑1 (MQ=255)
                                                                                                                 gTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATtggcgtggc                              >  3:198258/1‑149 (MQ=255)
                                                                                                                  tAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACaa                                                      >  2:283934/1‑124 (MQ=255)
                                                                                                                  tAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATtggcgtggcg                             >  2:71406/1‑149 (MQ=255)
                                                                                                                     gATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATAc                                                                             >  3:356176/1‑98 (MQ=255)
                                                                                                                     gATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATAc                                                                             <  4:356176/98‑1 (MQ=255)
                                                                                                                      aTTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATc                         <  1:464369/149‑1 (MQ=255)
                                                                                                                      aTTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCACCACCGGCAGATTGGCGTGGCGGATc                         <  5:341418/149‑1 (MQ=255)
                                                                                                                       ttCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCa                        >  6:200055/1‑149 (MQ=255)
                                                                                                                          aCCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGc                     >  7:170757/1‑149 (MQ=255)
                                                                                                                          aCCTGCAATAGATTAGCCGCTGGGGTCAAACATGCCGGCTTTTTTATCGCAATAAAAAGCCCCGCCGGGATGGCTCCCGGCAGGGCGGAATAATTAATCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGTCGGATAACGc                     >  8:313007/1‑149 (MQ=255)
                                                                                                                           ccTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCg                                             <  3:226365/124‑1 (MQ=255)
                                                                                                                           ccTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCg                                             >  4:226365/1‑124 (MQ=255)
                                                                                                                            cTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTc                   <  8:349060/149‑1 (MQ=255)
                                                                                                                               cAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAg                >  8:411835/1‑149 (MQ=255)
                                                                                                                                 ataaaaaaGCCGCCAGGTTTGACCCCGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGc              >  7:90414/8‑149 (MQ=255)
                                                                                                                                 aTTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATc                         >  8:50095/1‑138 (MQ=255)
                                                                                                                                 aTTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATc                         <  7:50095/138‑1 (MQ=255)
                                                                                                                                 aTTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGc              <  1:71406/149‑1 (MQ=255)
                                                                                                                                 aTTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGc              <  7:411835/149‑1 (MQ=255)
                                                                                                                                  ttGCTTAGCCGCCGGTGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCg             >  4:78896/1‑149 (MQ=255)
                                                                                                                                  ttGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCg             >  4:51264/1‑149 (MQ=255)
                                                                                                                                     cTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGt            >  7:393734/1‑147 (MQ=255)
                                                                                                                                     cTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTa          <  5:68069/149‑1 (MQ=255)
                                                                                                                                     cTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTa          <  3:316729/149‑1 (MQ=255)
                                                                                                                                     aaaaGCCGCCAGGTTTGACCCCGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGt            >  2:127454/4‑147 (MQ=255)
                                                                                                                                     aaaaGCCGCCAGGTTTGACCCCGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTa          <  2:24804/146‑1 (MQ=255)
                                                                                                                                      ttAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCTCCGGGATGGCTCCGGGCAGGGAGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATAGGCGTGGCGGATCACGCTCGAGGCGt            >  6:232991/1‑146 (MQ=255)
                                                                                                                                      ttAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGt            >  1:223627/1‑146 (MQ=255)
                                                                                                                                      ttAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGAAGGCTCCGGGCAGGGCGGGATACTTAATCTGGCAATTAACGCAAAACCATCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGt            >  8:421354/1‑146 (MQ=255)
                                                                                                                                            gccgGGGTCAGACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGc                     >  5:334297/1‑131 (MQ=255)
                                                                                                                                            gccgGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGAACCTa   <  2:106578/149‑1 (MQ=255)
                                                                                                                                             ccgGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCTGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGAACCTaa  <  3:146160/149‑1 (MQ=255)
                                                                                                                                             ccgGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGAACCTaa  <  7:439337/149‑1 (MQ=255)
                                                                                                                                             ccgGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGAACCTaa  <  3:147259/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
GTGACTGTACGACGGTTCCGTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGAACCTAA  >  NC_000913/641212‑641501

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: