Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,260,904 A→G 6.4% intergenic (‑101/+24) dauA ← / ← prs C4‑dicarboxylic acid transporter/phosphoribosylpyrophosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,260,9040AG 6.4% 455.2 / 15.5 140intergenic (‑101/+24)dauA/prsC4‑dicarboxylic acid transporter/phosphoribosylpyrophosphate synthase
Reads supporting (aligned to +/- strand):  ref base A (81/50);  new base G (5/4);  total (86/54)
Fisher's exact test for biased strand distribution p-value = 7.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.43e-01

Only 100 of 141 total aligned reads displayed.

AGAGCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTGATTTGATTTCATCGCTCAGCGGAA  >  NC_000913/1260760‑1261052
                                                                                                                                                |                                                                                                                                                    
agagCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCgg                                                                                                                                                  <  2:851965/149‑1 (MQ=255)
agagCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCgg                                                                                                                                                  >  2:579862/1‑149 (MQ=255)
agagCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCgg                                                                                                                                                  <  2:54819/149‑1 (MQ=255)
   gcgcGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTc                                                                                                                                               >  2:254648/1‑149 (MQ=255)
     gcgGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCtt                                                                                                                                             >  1:394025/1‑149 (MQ=255)
      cgGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCTGCGGCGGGTtttt                                                                                                                                            <  4:219626/149‑1 (MQ=255)
        gAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGa                                                                                                                                          >  3:88010/1‑149 (MQ=255)
               aTCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGt                                                                                                                                   >  1:760056/1‑149 (MQ=255)
                tCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGtt                                                                                                                                  <  2:507353/149‑1 (MQ=255)
                tCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGtt                                                                                                                                  >  2:797401/1‑149 (MQ=255)
                    aTGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGAt                                                                                                                              >  2:435888/1‑149 (MQ=255)
                    aTGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGAt                                                                                                                              >  1:69597/1‑149 (MQ=255)
                    aTGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATAATCTGTCAACTAAATGGCAGGTAAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGAt                                                                                                                              <  4:269250/149‑1 (MQ=255)
                        gAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAgt                                                                                                                          <  2:931348/149‑1 (MQ=255)
                           gAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAGACCCGCTGCGGCGGGTtttttt                                                                                                                                          >  1:97332/1‑128 (MQ=255)
                           gAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAGACCCGCTGCGGCGGGTtttttt                                                                                                                                          <  2:97332/130‑3 (MQ=255)
                           gAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGtt                                                                                                                       >  3:637604/1‑149 (MQ=255)
                           gAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGtt                                                                                                                       >  3:148383/1‑149 (MQ=255)
                           gAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGtt                                                                                                                       >  3:265768/1‑149 (MQ=255)
                               atatTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGaa                                                                                                                   <  1:333876/149‑1 (MQ=255)
                                   tttGTTCACAATAATGTTTAATCCGTGAGCGCATCATCTGTCAACTAAAGGCCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATg                                                                                                               <  3:773799/149‑1 (MQ=255)
                                          acaATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAgaga                                                                                                        >  1:449328/1‑149 (MQ=255)
                                            aataatGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGt                                                                                                                                   >  3:49287/1‑120 (MQ=255)
                                            aataatGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGt                                                                                                                                   <  4:49287/120‑1 (MQ=255)
                                            aataatGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATc                                                                                                      >  4:394054/1‑149 (MQ=255)
                                               aatGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGAt                                                                                                   <  1:170079/149‑1 (MQ=255)
                                                   tttCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGAttctt                                                                                               >  4:671478/1‑149 (MQ=255)
                                                   tttCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGAttctt                                                                                               >  4:63195/1‑149 (MQ=255)
                                                    ttCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGAttcttc                                                                                              <  2:464983/149‑1 (MQ=255)
                                                      cATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGt                                                                                            <  2:172985/149‑1 (MQ=255)
                                                      cATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGt                                                                                            >  1:546088/1‑149 (MQ=255)
                                                       aTCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAgt                                                                                                                          >  1:721244/1‑118 (MQ=255)
                                                       aTCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAgt                                                                                                                          <  2:721244/118‑1 (MQ=255)
                                                       aTCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAgt                                                                                                                          >  2:709587/1‑118 (MQ=255)
                                                       aTCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGa                                                                                                    >  2:344178/1‑140 (MQ=255)
                                                       aTCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGtt                                                                                           >  4:75249/1‑149 (MQ=255)
                                                       aTCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGtt                                                                                           >  2:494559/1‑149 (MQ=255)
                                                         ccGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGc                                                                                         <  3:424459/149‑1 (MQ=255)
                                                          cGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCt                                                                                        <  1:59665/149‑1 (MQ=255)
                                                             gAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGAt                                                                                     <  4:265768/149‑1 (MQ=255)
                                                               gcgcATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATac                                                                                   <  4:492910/149‑1 (MQ=255)
                                                                  catcatCGGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATacgac                                                                                >  2:505797/1‑149 (MQ=255)
                                                                   atcatcTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATacgacg                                                                               >  1:915744/1‑149 (MQ=255)
                                                                      atcTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAAt                                                                            <  3:814506/149‑1 (MQ=255)
                                                                      atcTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAAt                                                                            >  1:219009/1‑149 (MQ=255)
                                                                           tCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCtt                                                                       >  1:6450/1‑149 (MQ=255)
                                                                            cAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTc                                                                      <  1:466401/149‑1 (MQ=255)
                                                                             aaCTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGa                                                                                                                               <  3:215195/91‑1 (MQ=255)
                                                                              aCTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCgg                                                                    >  2:66157/1‑149 (MQ=255)
                                                                              aCTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCgg                                                                    >  3:235011/1‑149 (MQ=255)
                                                                              aCTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAAAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCgg                                                                    >  4:281555/1‑149 (MQ=255)
                                                                               cTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGc                                                                   >  1:746964/1‑149 (MQ=255)
                                                                                 aaaTGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCa                                                                 <  4:694498/149‑1 (MQ=255)
                                                                                  aaTGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAg                                                                >  1:714288/1‑149 (MQ=255)
                                                                                  aaTGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAg                                                                >  2:275967/1‑149 (MQ=255)
                                                                                      gCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATa                                                            >  4:143310/1‑149 (MQ=255)
                                                                                        aGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATAcc                                                          <  1:302712/149‑1 (MQ=255)
                                                                                         ggTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCt                                                         <  4:134544/149‑1 (MQ=255)
                                                                                           tGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGa                                                       >  3:177543/1‑149 (MQ=255)
                                                                                            gAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGAc                                                      >  2:775469/1‑149 (MQ=255)
                                                                                               ggaggCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAgg                                                   >  3:488941/1‑149 (MQ=255)
                                                                                                 aggCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGt                                                 <  3:521632/149‑1 (MQ=255)
                                                                                                  ggCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTc                                                <  2:449328/149‑1 (MQ=255)
                                                                                                     aTAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCaga                                             >  4:554348/1‑149 (MQ=255)
                                                                                                     aTAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCaga                                             <  3:433144/149‑1 (MQ=255)
                                                                                                       aGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGt                                           >  2:813077/1‑149 (MQ=255)
                                                                                                        ggTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTa                                          <  3:599300/149‑1 (MQ=255)
                                                                                                          tCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACg                                        >  3:405929/1‑149 (MQ=255)
                                                                                                             tACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCAc                                     >  1:150508/1‑149 (MQ=255)
                                                                                                                aaaTGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGtt                                  >  2:479541/1‑149 (MQ=255)
                                                                                                                 aaTGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTc                                 <  4:637604/149‑1 (MQ=255)
                                                                                                                 aaTGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTc                                 >  2:560162/1‑149 (MQ=255)
                                                                                                                  aTGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGAttcttc                                                                                              <  4:603886/87‑1 (MQ=255)
                                                                                                                  aTGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCg                                >  4:778034/1‑149 (MQ=255)
                                                                                                                   tGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCgg                               <  3:647824/149‑1 (MQ=255)
                                                                                                                    ggATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGc                              >  3:414297/1‑149 (MQ=255)
                                                                                                                    ggATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGCCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGc                              >  2:470609/1‑149 (MQ=255)
                                                                                                                       tatTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGt                           >  3:87201/1‑149 (MQ=255)
                                                                                                                          tACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGtgat                        <  2:760056/149‑1 (MQ=255)
                                                                                                                             aGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGtgatttg                     >  1:525061/1‑149 (MQ=255)
                                                                                                                               aCAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGtgatttgat                   >  1:877256/1‑149 (MQ=255)
                                                                                                                               aCAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGtgatttgat                   >  1:45804/1‑149 (MQ=255)
                                                                                                                                cAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGtgatttgatt                  <  2:915744/149‑1 (MQ=255)
                                                                                                                                 tcaaagaCCCGCTGCGGCGGGTTTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGAttcttc                                                                                              <  1:876902/70‑1 (MQ=25)
                                                                                                                                 tcaaagaCCCGCTGCGGCGGGTTTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGAttcttc                                                                                              >  2:876902/3‑72 (MQ=37)
                                                                                                                                 tcaaagaCCCGCTGCGGCGGGTTTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTgatttgattt                 >  3:787143/7‑149 (MQ=255)
                                                                                                                                 tcaaagaCCCGCTGCGGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTgatttgattt                 >  4:514719/3‑149 (MQ=255)
                                                                                                                                 aaaaaaaCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTgatttgattt                 <  1:579862/149‑1 (MQ=255)
                                                                                                                                  caaaaaCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGtgatttga                    >  1:385304/2‑145 (MQ=255)
                                                                                                                                   aaagaCCCGCTGCGGCGGGTTTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTGATTTGATTTCa               >  1:686230/1‑149 (MQ=255)
                                                                                                                                   aaaaaCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGtgatttgat                   <  2:356572/145‑1 (MQ=255)
                                                                                                                                   aaaaaCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGtgatttgat                   >  1:356572/1‑145 (MQ=255)
                                                                                                                                      aaCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTc                                                                                                                      >  4:356361/1‑43 (MQ=255)
                                                                                                                                      aaCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTc                                                                                                                      <  3:356361/43‑1 (MQ=255)
                                                                                                                                      aaCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCCTAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTGATTTGATTTCAtag            >  2:271831/1‑147 (MQ=255)
                                                                                                                                      aaCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGAttcttc                                                                                              >  3:641772/1‑67 (MQ=255)
                                                                                                                                      aaCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGAttcttc                                                                                              <  4:641772/67‑1 (MQ=255)
                                                                                                                                        cccGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTGATTTGATTTCATCGCt          >  3:559423/1‑149 (MQ=255)
                                                                                                                                         ccgccgCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTGATTTGATTTCATCGCTc         >  2:626559/1‑149 (MQ=255)
                                                                                                                                                aGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTGATTTGATTTCATCGCTCAGCGGaa  <  4:230328/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
AGAGCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTGATTTGATTTCATCGCTCAGCGGAA  >  NC_000913/1260760‑1261052

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: