Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,062,371 A→T 7.8% N42I (AAT→ATT)  yihN → MFS transporter family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,062,3710AT 7.8% 351.5 / 16.6 102N42I (AAT→ATT) yihNMFS transporter family protein
Reads supporting (aligned to +/- strand):  ref base A (47/47);  new base T (6/2);  total (53/49)
Fisher's exact test for biased strand distribution p-value = 2.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.13e-01

Only 100 of 105 total aligned reads displayed.

ATCGCTTAGGGGATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTGGGTGTTTATTTGACG  >  NC_000913/4062230‑4062519
                                                                                                                                              |                                                                                                                                                    
aTCGCTTAGGGGATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCaa                                                                                                                                               <  4:33568/149‑1 (MQ=255)
 tCGCTTAGGGGATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCaaa                                                                                                                                              >  3:488406/1‑149 (MQ=255)
    cTTAGGGGATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTg                                                                                                                                           <  2:310475/149‑1 (MQ=255)
    cTTAGGGGATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTg                                                                                                                                           >  1:310475/1‑149 (MQ=255)
            aTGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCt                                                                                                                                   <  4:488406/149‑1 (MQ=255)
                           aaaaGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTa                                                                                                                                        <  1:775676/129‑1 (MQ=255)
                           aaaaGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTa                                                                                                                                        >  2:775676/1‑129 (MQ=255)
                                aaaTGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAAGTTATAAATTACCGTCGCTGAAAGATGCGTTTTACATCCCGATGCAGGAATATTCCCATTTGACCAATGGTCAAATTGGTCATGCTATGTCGGTAAACTCATTCGt                                                                                                               <  4:911035/149‑1 (MQ=255)
                                  aTGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGACGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTca                                                                                                             >  4:336383/1‑149 (MQ=255)
                                                gTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTgg                                                                                                                                          <  2:1280407/106‑1 (MQ=255)
                                                gTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTgg                                                                                                                                          >  1:1280407/1‑106 (MQ=255)
                                                gTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCttt                                                                                               >  3:304521/1‑149 (MQ=255)
                                                gTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCttt                                                                                               >  2:42094/1‑149 (MQ=255)
                                                   tATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGt                                                                                                                                         <  4:1162081/104‑1 (MQ=255)
                                                   tATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGt                                                                                                                                         >  3:1162081/1‑104 (MQ=255)
                                                    aTTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAAt                                                                                                                                                     >  2:1431012/1‑91 (MQ=255)
                                                    aTTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAAt                                                                                                                                                     <  1:1431012/91‑1 (MQ=255)
                                                     ttAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAACATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCt                                                                                          >  4:377437/1‑149 (MQ=255)
                                                       aaCACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGACGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGt                                                                                        <  3:336383/149‑1 (MQ=255)
                                                        acacTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTc                                                                                       >  2:1310995/1‑149 (MQ=255)
                                                          acTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCta                                                                                     >  2:541126/1‑149 (MQ=255)
                                                          acTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCta                                                                                     >  4:287060/1‑149 (MQ=255)
                                                           cTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCtat                                                                                    <  1:1310995/149‑1 (MQ=255)
                                                             gtgtgGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCttt                                                                                               <  3:392514/136‑1 (MQ=255)
                                                             gtgtgGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCttt                                                                                               >  4:392514/1‑136 (MQ=255)
                                                             gtgtgGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGt                                                                                        >  1:259372/1‑143 (MQ=255)
                                                             gtgtgGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGt                                                                                        <  2:259372/143‑1 (MQ=255)
                                                                 ggCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTatttattt                                                                              <  1:541126/149‑1 (MQ=255)
                                                                  gCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGt                                                                                                                                         <  2:318363/89‑1 (MQ=255)
                                                                  gCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGt                                                                                                                                         >  1:318363/1‑89 (MQ=255)
                                                                     gTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCt                                                                                          <  1:1155078/133‑1 (MQ=255)
                                                                     gTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCt                                                                                          >  2:1155078/1‑133 (MQ=255)
                                                                     gTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGcc                                                                          <  1:42094/149‑1 (MQ=255)
                                                                      tACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCg                                                                         >  3:118129/1‑149 (MQ=255)
                                                                       aCAATTTATAAATTACCGTCGCTGAAACATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGa                                                                        <  3:377437/149‑1 (MQ=255)
                                                                         aaTTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTACCAATTTGACCATTGGTCAAAt                                                                                                                                             >  3:1395904/1‑79 (MQ=37)
                                                                          aTTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCtttttt                                                                                            <  3:632843/126‑1 (MQ=255)
                                                                          aTTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCtttttt                                                                                            >  4:632843/1‑126 (MQ=255)
                                                                           tttATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAAt                                                                                                                                      <  3:96169/83‑1 (MQ=255)
                                                                                  aTTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTg                                                                                                                                           <  1:705414/71‑1 (MQ=255)
                                                                                  aTTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTg                                                                                                                                           >  2:705414/1‑71 (MQ=255)
                                                                                    tACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCgcgc                                                           >  3:1466918/1‑149 (MQ=255)
                                                                                       cGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGa                                                        <  4:1466918/149‑1 (MQ=255)
                                                                                             tGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCtattta                                                                                 <  2:992797/118‑1 (MQ=255)
                                                                                             tGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCtattta                                                                                 >  1:992797/1‑118 (MQ=255)
                                                                                               aaaGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCAt                                                <  4:304521/149‑1 (MQ=255)
                                                                                                   aTGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAAt                                                                                                                                                     >  3:1268056/1‑44 (MQ=255)
                                                                                                   aTGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAAt                                                                                                                                                     <  4:1268056/44‑1 (MQ=255)
                                                                                                    tGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCt                                                                                          <  4:1108473/102‑1 (MQ=255)
                                                                                                    tGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCt                                                                                          <  2:1435154/102‑1 (MQ=255)
                                                                                                    tGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCt                                                                                          >  1:1435154/1‑102 (MQ=255)
                                                                                                    tGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCt                                                                                          >  3:1108473/1‑102 (MQ=255)
                                                                                                          ttATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTg                                                                                                     >  2:239235/1‑85 (MQ=255)
                                                                                                          ttATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTg                                                                                                     <  1:239235/85‑1 (MQ=255)
                                                                                                              atCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTAtttt                                                                             >  3:1385879/1‑105 (MQ=255)
                                                                                                              atCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTAtttt                                                                             <  4:1385879/105‑1 (MQ=255)
                                                                                                              atCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAAc                                                                   <  4:439252/115‑1 (MQ=255)
                                                                                                              atCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAAc                                                                   >  3:439252/1‑115 (MQ=255)
                                                                                                               tCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCAtt                                <  2:640693/149‑1 (MQ=255)
                                                                                                               tCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCAtt                                >  1:175691/1‑149 (MQ=255)
                                                                                                               tCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCAtt                                >  3:233186/1‑149 (MQ=255)
                                                                                                               tCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCAtt                                >  1:640693/1‑149 (MQ=255)
                                                                                                                 ccGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGc                              <  4:118129/149‑1 (MQ=255)
                                                                                                                  cGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCg                             >  3:452190/1‑149 (MQ=255)
                                                                                                                   gATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCAtt                                >  2:333132/1‑145 (MQ=255)
                                                                                                                   gATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCAtt                                <  1:333132/145‑1 (MQ=255)
                                                                                                                      gCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTggg                                                                                                   <  1:500005/75‑1 (MQ=255)
                                                                                                                        aGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTctca                                      >  1:589968/1‑134 (MQ=255)
                                                                                                                        aGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTctca                                      >  3:465975/1‑134 (MQ=255)
                                                                                                                        aGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTctca                                      <  4:465975/134‑1 (MQ=255)
                                                                                                                        aGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTctca                                      <  2:589968/134‑1 (MQ=255)
                                                                                                                        aGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTctca                                      <  2:621002/134‑1 (MQ=255)
                                                                                                                          gAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGAtt                     >  1:125890/1‑149 (MQ=255)
                                                                                                                            atatTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCt                                                                                          >  3:817612/1‑78 (MQ=255)
                                                                                                                            atatTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCt                                                                                          <  4:817612/78‑1 (MQ=255)
                                                                                                                               tttCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATtct                                        >  1:521009/1‑125 (MQ=255)
                                                                                                                               tttCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATtct                                        <  2:521009/125‑1 (MQ=255)
                                                                                                                               tttCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTgg                <  4:233186/149‑1 (MQ=255)
                                                                                                                               tttCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTgg                <  2:125890/149‑1 (MQ=255)
                                                                                                                                ttCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAACCTACCGCGCAGATACACCATGTCATTCTCACTCAt                                 >  1:1256696/1‑131 (MQ=255)
                                                                                                                                ttCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCAt                                 <  2:1256696/131‑1 (MQ=255)
                                                                                                                                ttCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTggg               <  4:452190/149‑1 (MQ=255)
                                                                                                                                   cATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATacac                                                   <  1:105829/110‑1 (MQ=255)
                                                                                                                                   cATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATacac                                                   >  2:105829/1‑110 (MQ=255)
                                                                                                                                    aaTTTGACCATTGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTGGgtgt            >  3:600909/2‑149 (MQ=255)
                                                                                                                                     aTTTGACCATTGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCAtt                                          >  4:721689/1‑118 (MQ=255)
                                                                                                                                     aTTTGACCATTGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCAtt                                          <  3:721689/118‑1 (MQ=255)
                                                                                                                                     aTTTGACCATTGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTGGGTGtt           >  4:202083/1‑149 (MQ=255)
                                                                                                                                     aTTTGACCATTGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTGGGTGtt           >  1:229750/1‑149 (MQ=255)
                                                                                                                                     aTTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTctca                                      <  4:141153/122‑1 (MQ=255)
                                                                                                                                     aTTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTctca                                      >  3:141153/1‑122 (MQ=255)
                                                                                                                                      tttGACCATTGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTGGGTGttt          >  2:1221510/1‑149 (MQ=255)
                                                                                                                                      tttGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTGGGTGttt          >  3:923189/1‑149 (MQ=255)
                                                                                                                                         gACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATAGACCATGTCATTCTCACTCATTGCGACAg                         >  3:1024018/1‑131 (MQ=255)
                                                                                                                                         gACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAg                         <  4:1024018/131‑1 (MQ=255)
                                                                                                                                            cAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATaa                                                                     <  2:727722/84‑1 (MQ=255)
                                                                                                                                            cAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATaa                                                                     >  1:727722/1‑84 (MQ=255)
                                                                                                                                              aTGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTGGgt              >  3:1427591/1‑137 (MQ=255)
                                                                                                                                              aTGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTGGgt              <  4:1427591/137‑1 (MQ=255)
                                                                                                                                              aTGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTGGGTGTTTATTTgacg  >  4:229052/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
ATCGCTTAGGGGATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCTATTAACACTGTGTGGCGGTACAATTTATAAATTACCGTCGCTGAAAGATGCGTTTTATATCCCGATGCAGGAATATTTCCATTTGACCAATGGTCAAATTGGTAATGCTATGTCGGTAAACTCATTTGTCACCACAGTGGGCTTTTTTCTGTCTATTTATTTTGCCGATAAACTACCGCGCAGATACACCATGTCATTCTCACTCATTGCGACAGGATTACTGGGTGTTTATTTGACG  >  NC_000913/4062230‑4062519

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