Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,720,798 G→T 6.0% R40M (AGG→ATG) ‡ insJ → IS150 transposase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,720,7980GT 6.0% 589.3 / 20.4 166R40M (AGG→ATG) ‡insJIS150 transposase A
Reads supporting (aligned to +/- strand):  ref base G (74/82);  new base T (2/8);  total (76/90)
Fisher's exact test for biased strand distribution p-value = 1.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.31e-01

Only 100 of 167 total aligned reads displayed.

GTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCTGCTGCTCACTT  >  NC_000913/3720650‑3720945
                                                                                                                                                    |                                                                                                                                                   
gttgGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAg                                                                                                                                                     <  4:968758/149‑1 (MQ=255)
   ggaggaTGAAATGGAATAGCCCCTAATATGTCAATGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACACGATCGTCTCGGCACGTTTTGGTTTCCCCTGAACCCAGGTCAGGAc                                                                                                                                                  >  4:1028728/1‑149 (MQ=255)
       gATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCTTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATgg                                                                                                                                              >  4:623884/1‑149 (MQ=255)
       gATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCATGTCCTGACCTgg                                                                                                                                              <  3:346244/149‑1 (MQ=255)
       gATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATgg                                                                                                                                              <  2:898714/149‑1 (MQ=255)
             aTGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCATGTCCTGACCTGGGttcgg                                                                                                                                        <  3:76259/149‑4 (MQ=255)
             aTGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGcc                                                                                                                                        >  2:116473/1‑149 (MQ=255)
              tGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCATGTCCTGACCTGGGttcggg                                                                                                                                       <  1:895594/149‑5 (MQ=255)
              tGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGccc                                                                                                                                       >  2:1008780/1‑149 (MQ=255)
               ggAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCATGTCCTGACCTGGGttcgggg                                                                                                                                      <  3:149472/149‑6 (MQ=255)
               ggAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCATGTCCTGACCTGGGttcgggg                                                                                                                                      <  2:576644/149‑6 (MQ=255)
               ggAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCATGTCCTGACCTGGGttcgggg                                                                                                                                      <  4:665296/149‑6 (MQ=255)
                 aaTAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCctct                                                                                                                                    <  2:577764/149‑1 (MQ=255)
                  aTAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTa                                                                                                                                   <  2:656660/149‑1 (MQ=255)
                       cccTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGaaa                                                                                                                              <  3:5122/149‑1 (MQ=255)
                        ccTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGaaaa                                                                                                                             >  1:362079/1‑149 (MQ=255)
                            atatGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACAt                                                                                                                         <  4:547566/149‑1 (MQ=255)
                            atatGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACAt                                                                                                                         <  4:864447/149‑1 (MQ=255)
                              atGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATgg                                                                                                                       >  1:1053966/1‑149 (MQ=255)
                               tGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGga                                                                                                                      <  3:753787/149‑1 (MQ=255)
                                gTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGgag                                                                                                                     <  2:906953/149‑1 (MQ=255)
                                  cacggCCAAAATACCCTTTTGAAAAGCGCCATGAAGTCGTAAATCACTACTTCACAACTGATGATGCTTACAGGTTCATCGCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGTACATCGGTTGCCCTCTGTGATTAACATGGAGaa                                                                                                                   >  2:1080361/5‑149 (MQ=255)
                                     aGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGaaaaa                                                                                                                <  3:348047/149‑1 (MQ=255)
                                     aGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGaaaaa                                                                                                                <  2:362079/149‑1 (MQ=255)
                                      gCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAg                                                                                                               >  1:782969/1‑149 (MQ=255)
                                        cAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGt                                                                                                             >  3:995459/1‑149 (MQ=255)
                                         aaaaTACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGtt                                                                                                            >  3:883213/1‑149 (MQ=255)
                                              aCCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACCGGATCATCTCGCCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAAtt                                                                                                       <  2:786006/149‑1 (MQ=255)
                                                  ttttGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCATGTCCTGACCTGGGttcggg                                                                                                                                       <  3:36683/113‑5 (MQ=255)
                                                  ttttGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCATGTCCTGACCTGGGttcggg                                                                                                                                       >  4:36683/1‑109 (MQ=255)
                                                  ttttGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCa                                                                                                   >  4:653362/1‑149 (MQ=255)
                                                      gAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAAcc                                                                                               >  3:268799/1‑149 (MQ=255)
                                                       aaaaGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCt                                                                                              >  3:1069211/1‑149 (MQ=255)
                                                          aGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTaaa                                                                                           <  4:607085/149‑1 (MQ=255)
                                                          aGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTaaa                                                                                           <  1:116473/149‑1 (MQ=255)
                                                           gcgcCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAg                                                                                          <  4:180630/149‑1 (MQ=255)
                                                                 tGAAGTAGTGATTCACGACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGt                                                                                                             <  1:97474/124‑1 (MQ=255)
                                                                   aaGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCTGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGt                                                                                  <  3:677057/149‑1 (MQ=255)
                                                                        gTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGa                                                                             <  4:714271/149‑1 (MQ=255)
                                                                              cACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCaga                                                                       <  4:789006/149‑1 (MQ=255)
                                                                              cACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCaga                                                                       <  2:548280/149‑1 (MQ=255)
                                                                                    ttCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGt                                                                                                             <  1:205569/105‑1 (MQ=255)
                                                                                    ttCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCg                                                                 >  3:251153/1‑149 (MQ=255)
                                                                                    ttCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCg                                                                 >  4:757947/1‑149 (MQ=255)
                                                                                    ttCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCg                                                                 >  1:553386/1‑149 (MQ=255)
                                                                                    ttCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCg                                                                 >  2:240242/1‑149 (MQ=255)
                                                                                     tCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGc                                                                  >  2:630481/1‑147 (MQ=255)
                                                                                     tCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGt                                                                <  2:1053966/149‑1 (MQ=255)
                                                                                      cacaACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTa                                                               <  1:151497/149‑1 (MQ=255)
                                                                                      cacaACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATAGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTa                                                               <  4:640584/149‑1 (MQ=255)
                                                                                       acaACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTAt                                                              >  4:849737/1‑149 (MQ=255)
                                                                                        caACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTAtt                                                             <  3:504503/149‑1 (MQ=255)
                                                                                         aaCTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCctc                                                                                                                                     >  2:389267/1‑76 (MQ=255)
                                                                                         aaCTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGtt                                                                                                            >  1:63375/1‑101 (MQ=255)
                                                                                          aCTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTaa                                                           >  3:252152/1‑149 (MQ=255)
                                                                                              atgatgGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTc                                                       <  1:697758/149‑1 (MQ=255)
                                                                                                 atgGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCt                                                                               <  3:172631/122‑1 (MQ=255)
                                                                                                   ggTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGaa                                                  >  3:229119/1‑149 (MQ=255)
                                                                                                     ttACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAg                                                >  2:999119/1‑149 (MQ=255)
                                                                                                       aCAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCt                                              <  2:1024122/149‑1 (MQ=255)
                                                                                                         aGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCtgt                                            <  3:887051/149‑1 (MQ=255)
                                                                                                           gATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGaaaaa                                                                                                                <  3:837820/79‑1 (MQ=255)
                                                                                                            atcatcTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGAt                                         <  4:151992/149‑1 (MQ=255)
                                                                                                             tcatcTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATc                                        <  4:268799/149‑1 (MQ=255)
                                                                                                             tcatcTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATc                                        <  1:1011494/149‑1 (MQ=255)
                                                                                                             tcatcTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATc                                        >  2:886414/1‑149 (MQ=255)
                                                                                                               atcTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGttt                                                                                                           >  3:766807/1‑80 (MQ=255)
                                                                                                               atcTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGttt                                                                                                           <  4:766807/80‑1 (MQ=255)
                                                                                                               atcTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGa                                      <  3:156977/149‑1 (MQ=255)
                                                                                                                 ctcGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGagc                                    <  2:169319/149‑1 (MQ=255)
                                                                                                                 ctcGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGagc                                    <  2:276877/149‑1 (MQ=255)
                                                                                                                  tcGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGagca                                   >  1:933767/1‑149 (MQ=255)
                                                                                                                  tcGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGagca                                   >  3:65725/1‑149 (MQ=255)
                                                                                                                  tcGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGagca                                   <  4:783867/149‑1 (MQ=255)
                                                                                                                  tcGGCACGATTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGagca                                   <  2:903055/149‑1 (MQ=255)
                                                                                                                    ggCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAAtt                                                                                                       >  4:808694/1‑79 (MQ=255)
                                                                                                                    ggCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAAtt                                                                                                       <  3:808694/79‑1 (MQ=255)
                                                                                                                     gCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAgcagca                                >  2:831760/1‑149 (MQ=255)
                                                                                                                      cACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGTCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTTAAAGCTGTGATCGAGCAGcac                               <  2:236631/149‑1 (MQ=255)
                                                                                                                      cACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGcac                               <  1:831760/149‑1 (MQ=255)
                                                                                                                       aCGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGcaca                              <  1:604895/149‑1 (MQ=255)
                                                                                                                       aCGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGcaca                              <  4:245469/149‑1 (MQ=255)
                                                                                                                          ttttGGTGTCCCCCGGACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGt                           <  3:543933/149‑1 (MQ=255)
                                                                                                                            ttGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTcc                         <  1:749171/149‑1 (MQ=255)
                                                                                                                                    ccccGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGc                                                                                        >  2:694591/1‑78 (MQ=255)
                                                                                                                                    ccccGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGc                                                                                        <  1:694591/78‑1 (MQ=255)
                                                                                                                                    ccccGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGt                                                     <  1:816582/113‑1 (MQ=255)
                                                                                                                                     cccGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGc                                                                  <  1:906603/99‑1 (MQ=255)
                                                                                                                                     cccGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAg                <  3:523710/149‑1 (MQ=255)
                                                                                                                                     cccGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAAGTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAg                >  3:625046/1‑149 (MQ=255)
                                                                                                                                      ccGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAgg               >  3:543550/1‑149 (MQ=255)
                                                                                                                                      ccGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGAAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAg                                                          >  2:84560/1‑106 (MQ=255)
                                                                                                                                      ccGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGAAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAg                                                          <  1:84560/106‑1 (MQ=255)
                                                                                                                                      ccGAACCCAGGTCAGGACATGGGTTGCCCTCAATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAgg               <  4:111618/149‑1 (MQ=255)
                                                                                                                                             cAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGgctgctgc        >  1:405515/1‑149 (MQ=255)
                                                                                                                                                gTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCaga                                                                       <  3:780911/83‑1 (MQ=255)
                                                                                                                                                gTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCTGCTGCTCa     >  4:19092/1‑149 (MQ=255)
                                                                                                                                                   aGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCCGCACATGTCCCTTAATCAGGCTGCTGCTCACtt  <  3:920321/149‑1 (MQ=255)
                                                                                                                                                   aGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCTGCTGCTCACtt  <  4:480696/149‑1 (MQ=255)
                                                                                                                                                   aGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCTGCTGCTCACtt  >  4:909145/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCTGCTGCTCACTT  >  NC_000913/3720650‑3720945

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