Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,415,725 C→A 7.3% L327M (CTG→ATG) ‡ pta → phosphate acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,415,7250CA 7.3% 380.0 / 15.3 110L327M (CTG→ATG) ‡ptaphosphate acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (41/61);  new base A (3/5);  total (44/66)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.18e-01

Only 100 of 111 total aligned reads displayed.

GACTTCCGCAGACCGTCCTGACGTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAATC  >  NC_000913/2415580‑2415873
                                                                                                                                                 |                                                                                                                                                    
gACTTCCGCAGACCGTCCTGACGTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGc                                                                                                                                                   >  4:660520/1‑149 (MQ=255)
gACTTCCGCAGACCGTCCTGACGTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGc                                                                                                                                                   >  3:365352/1‑149 (MQ=255)
 aCTTCCGCAGACCGTCCTGACGTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGcc                                                                                                                                                  <  4:126906/149‑1 (MQ=255)
  cTTCCGCAGACCGTCCTGACGTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCg                                                                                                                                                 <  2:298385/149‑1 (MQ=255)
  cTTCCGCAGACCGTCCTGACGTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCg                                                                                                                                                 <  4:34726/149‑1 (MQ=255)
   ttCCGCAGACCGTCCTGACGTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCgg                                                                                                                                                >  1:148746/1‑149 (MQ=255)
    tCCGCAGACCGTCCTGACGTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGt                                                                                                                                               <  1:107367/149‑1 (MQ=255)
     ccGCAGACCGTCCTGACGTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCAGGCCGGTa                                                                                                                                              <  1:515079/149‑1 (MQ=255)
      cGCAGACCGTCCTGACGTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTAt                                                                                                                                             <  1:311528/149‑1 (MQ=255)
      cGCAGACCGTCCTGACGTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCAGGCCGGTag                                                                                                                                             <  1:337293/149‑2 (MQ=255)
       gCAGACCGTCCTGACGTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTAtt                                                                                                                                            >  4:450925/1‑149 (MQ=255)
       gCAGACCGTCCTGACGTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGGTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTAtt                                                                                                                                            >  4:575284/1‑149 (MQ=255)
           aCCGTCCTGACGTGCTGGTGGCCGCTTGCGTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGTGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCAGGCCGGTagcgaaa                                                                                                                                        <  2:322537/149‑7 (MQ=255)
           aCCGTCCTGACGTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCAGGCCGGTagcgaaa                                                                                                                                        <  3:334957/149‑7 (MQ=255)
            ccGTCCTGACGTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATgg                                                                                                                                       <  2:248145/149‑1 (MQ=255)
            ccGTCCTGACGTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATgg                                                                                                                                       <  1:597598/149‑1 (MQ=255)
            ccGTCCTGACGTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATgg                                                                                                                                       <  4:226524/149‑1 (MQ=255)
             cGTCCTGACGTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGt                                                                                                                                      <  4:196353/149‑1 (MQ=255)
             cGTCCTGACGTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGt                                                                                                                                      <  4:122756/149‑1 (MQ=255)
                   gACGTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGaacac                                                                                                                                <  4:563795/149‑1 (MQ=255)
                      gTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGcc                                                                                                                                                  >  3:62747/1‑128 (MQ=255)
                      gTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGcc                                                                                                                                                  <  4:62747/128‑1 (MQ=255)
                      gTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGaacaccaa                                                                                                                             <  4:76720/149‑1 (MQ=255)
                          tggtggCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGaacaccaacacc                                                                                                                         >  3:665977/1‑149 (MQ=255)
                           ggtggCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCt                                                                                                                        <  1:465862/149‑1 (MQ=255)
                            gtggCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTg                                                                                                                       <  3:288961/149‑1 (MQ=255)
                             tggCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTgg                                                                                                                      <  4:384108/149‑1 (MQ=255)
                               gCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCa                                                                                                                    >  2:531848/1‑149 (MQ=255)
                                 cGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGa                                                                                                                  >  2:668002/1‑149 (MQ=255)
                                 cGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGa                                                                                                                  >  2:526155/1‑149 (MQ=255)
                                  gCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGAc                                                                                                                 >  1:118975/1‑149 (MQ=255)
                                    ttGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACct                                                                                                               <  2:499490/149‑1 (MQ=255)
                                        cTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACctctct                                                                                                           >  4:316688/1‑149 (MQ=255)
                                        cTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACctctct                                                                                                           <  4:359130/149‑1 (MQ=255)
                                         tGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCAGGCCGGTag                                                                                                                                             >  3:360424/1‑113 (MQ=255)
                                         tGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCAGGCCGGTag                                                                                                                                             <  4:360424/114‑2 (MQ=255)
                                          ggCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGa                                                                                                         <  3:202057/149‑1 (MQ=255)
                                               ccATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGa                                                                                                         <  1:394965/144‑1 (MQ=255)
                                               ccATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGa                                                                                                         >  2:394965/1‑144 (MQ=255)
                                                cATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGc                                                                                                   >  1:544776/1‑149 (MQ=255)
                                                   gAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCaga                                                                                                <  1:189271/149‑1 (MQ=255)
                                                      cGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCt                                                                                             <  1:438075/149‑1 (MQ=255)
                                                      cGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCt                                                                                             <  2:544776/149‑1 (MQ=255)
                                                      cGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCt                                                                                             <  3:660520/149‑1 (MQ=255)
                                                          gTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCaa                                                                                         >  2:157933/1‑149 (MQ=255)
                                                           tAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAAc                                                                                        <  3:524806/149‑1 (MQ=255)
                                                              aaaTCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTg                                                                                     <  1:624116/149‑1 (MQ=255)
                                                               aaTCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTgg                                                                                    >  1:148357/1‑149 (MQ=255)
                                                                 tCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGaa                                                                                  <  2:334831/149‑1 (MQ=255)
                                                                   ggTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGc                                                                                              >  3:293853/1‑135 (MQ=255)
                                                                   ggTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGc                                                                                              <  4:293853/135‑1 (MQ=255)
                                                                             tgctgACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGAc                                                                      <  4:606961/149‑1 (MQ=255)
                                                                               ctgACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGa                                                                    >  2:417454/1‑149 (MQ=255)
                                                                                tgACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGAt                                                                   <  2:118975/149‑1 (MQ=255)
                                                                                  aCTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCa                                                                 >  2:265511/1‑149 (MQ=255)
                                                                                    tGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACg                                                               >  1:409114/1‑149 (MQ=255)
                                                                                    tGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACg                                                               >  1:76681/1‑149 (MQ=255)
                                                                                     ggCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGa                                                              <  1:265511/149‑1 (MQ=255)
                                                                                          ttACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTa                                                         >  3:142371/1‑149 (MQ=255)
                                                                                           tACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTAt                                                        <  1:492419/149‑1 (MQ=255)
                                                                                            aCGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATc                                                       <  4:365352/149‑1 (MQ=255)
                                                                                              gAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCga                                                     >  3:255976/1‑149 (MQ=255)
                                                                                              gAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCga                                                     >  2:234737/1‑149 (MQ=255)
                                                                                               aaaTGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCgag                                                    <  4:255976/149‑1 (MQ=255)
                                                                                                  tGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGaaa                                                 <  3:524738/149‑1 (MQ=255)
                                                                                                       gcgcgcATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATgg                                                                                                                                       >  3:475978/1‑58 (MQ=255)
                                                                                                       gcgcgcATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCa                                            <  1:531848/149‑1 (MQ=255)
                                                                                                       gcgcgcATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCa                                            >  2:535484/1‑149 (MQ=255)
                                                                                                       gcgcgcATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCa                                            <  3:677466/149‑1 (MQ=255)
                                                                                                       gcgcgcATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCa                                            <  4:111467/149‑1 (MQ=255)
                                                                                                       gcgcgcATTTCTAAACTGTGCGAACGCGCTTTCGCTACCGGCCTGCCGGTATTTATgg                                                                                                                                       <  4:475978/58‑1 (MQ=255)
                                                                                                         gcgcATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAg                                                <  3:155460/143‑1 (MQ=255)
                                                                                                         gcgcATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAg                                                >  4:155460/1‑143 (MQ=255)
                                                                                                           gcATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGaa                                        >  2:295134/1‑149 (MQ=255)
                                                                                                            cATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAAt                                       >  4:274543/1‑149 (MQ=255)
                                                                                                                 cTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCaga                                                                                                >  1:60686/1‑87 (MQ=255)
                                                                                                                 cTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCaga                                                                                                <  2:60686/87‑1 (MQ=255)
                                                                                                                 cTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGtt                                  <  4:665977/149‑1 (MQ=255)
                                                                                                                  tAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGa                                                                    <  2:580964/114‑1 (MQ=255)
                                                                                                                  tAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGa                                                                    <  2:368379/114‑1 (MQ=255)
                                                                                                                   aaaCTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGc                                <  3:410696/149‑1 (MQ=255)
                                                                                                                      cTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTaa                             <  4:142371/149‑1 (MQ=255)
                                                                                                                         tgCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTa                          <  3:275211/149‑1 (MQ=255)
                                                                                                                           cGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCCACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACa                        <  3:129686/149‑1 (MQ=255)
                                                                                                                           cGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACCAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACa                        <  2:409114/149‑1 (MQ=255)
                                                                                                                             aaCGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATc                      <  2:76681/149‑1 (MQ=255)
                                                                                                                                gTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCt                                                                                             <  1:410023/75‑1 (MQ=255)
                                                                                                                                gTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCt                                                                                             >  2:410023/1‑75 (MQ=255)
                                                                                                                                gTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAAc                   <  3:450925/149‑1 (MQ=255)
                                                                                                                                   cTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCt                <  1:157933/149‑1 (MQ=255)
                                                                                                                                   cTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCt                >  1:672554/1‑149 (MQ=255)
                                                                                                                                    tttCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGtt                                  >  4:567790/1‑130 (MQ=255)
                                                                                                                                    tttCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGtt                                  <  3:567790/130‑1 (MQ=255)
                                                                                                                                    tttCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTg               >  1:138666/1‑149 (MQ=255)
                                                                                                                                     ttCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGa              <  2:635646/149‑1 (MQ=255)
                                                                                                                                      tCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGAc             <  1:87637/149‑1 (MQ=255)
                                                                                                                                         atACCGGCAGGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTga          >  1:660971/2‑148 (MQ=255)
                                                                                                                                           aCCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGGCCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGAt        <  4:626335/149‑1 (MQ=255)
                                                                                                                                            ccGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATc       >  1:78771/1‑149 (MQ=255)
                                                                                                                                                 cTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAAtc  >  3:184257/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GACTTCCGCAGACCGTCCTGACGTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAATC  >  NC_000913/2415580‑2415873

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