Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 786,241 C→T 6.4% A203A (GCC→GCT aroG → 3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, phenylalanine repressible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913786,2410CT 6.4% 213.9 / 3.9 62A203A (GCC→GCTaroG3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, phenylalanine repressible
Reads supporting (aligned to +/- strand):  ref base C (27/31);  new base T (2/2);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAACTGGCATCAGGGCTTTCTTGTCCGGTCGGCTTCAAAAATGGCACCGACGGTACGATTAAAGTGGCTATCGATGCCATTAATGCCGCCGGTGCGCCGCACTGCTTCCTGTCCGTAACGAAATGGGGGCATTCGGCGATTGTGAATACCAGCGGTAACGGCGATTGCCATATCATTCTGCGCGGCGGTAAAGAGCCTAACTACAGCGCGAAGCACGTTGCTGAAGTGA  >  NC_000913/786106‑786380
                                                                                                                                       |                                                                                                                                           
cTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAACTGGCATCAGGGCTTTCTTGTCCGGTCGGCTTCAAAAATGGCACCGACGGTACGATTAAAGTGGCTATCGATGCCATTAATGCCGCCGGTg                                                                                                                                         >  1:1484737/1‑140 (MQ=255)
cTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAACTGGCATCAGGGCTTTCTTGTCCGGTCGGCTTCAAAAATGGCACCGACGGTACGATTAAAGTGGCTATCGATGCCATTAATGCCGCCGGTg                                                                                                                                         >  1:1484761/1‑140 (MQ=255)
cTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAACTGGCATCAGGGCTTTCTTGTCCGGTCGGCTTCAAAAATGCCACCGACGGTACGATTAAAGTGGCTATCGATGCCATTAATGCCGCCGGTg                                                                                                                                         >  1:24033/1‑140 (MQ=255)
  ggggCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAACTGGCATCAGGGCTTTCTTGTCCGGTCGGCTTCAAAAATGGCACCGACGGTACGATTAAAGTGGCTATCGATGCCATTAATGCCGCCGGTgcg                                                                                                                                       >  1:551204/1‑140 (MQ=255)
  ggggCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAACTGGCATCAGGGCTTTCTTGTCCGGTCGGCTTCAAAAATGGCACCGACGGTACGATTAAAGTGGCTATCGATGCCATTAATGCCGCCGGTgcg                                                                                                                                       >  2:380682/1‑140 (MQ=255)
  ggggCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAACTGGCATCAGGGCTTTCTTGTCCGGTCGGCTTCAAAAATGGCACCGACGGTACGATTAAAGTGGCTATCGATGCCATTAATGCCGCCGGTgcg                                                                                                                                       >  1:526915/1‑140 (MQ=255)
  ggggCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAACTGGCATCAGGGCTTTCTTGTCCGGTCGGCTTCAAAAATGGCACCGACGGTACGATTAAAGTGGCTATCGATGCCAGTAATGCCGCCGGTgcg                                                                                                                                       >  2:501544/1‑140 (MQ=255)
  ggggCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAACTGGCATCAGGGCTTTCTTGTCCGGTCGGCGTCAAAAATGGCACCGACGGTACGATTAAAGTGGCTAGCGATGCCATTAATGCCGCCGGTgcg                                                                                                                                       >  2:501522/1‑140 (MQ=255)
     gcgcAATTGGCGCACGTACCACCGAATCGCTGGTGCACCGCGAACTGGCATCAGGGCTTTCTTGTCCGGTCGGCTTCTAAAATGGCACCGACGGTACGATTAAAGTGGCTATCGATGCCATTAATGCCGCCGGTGcgccg                                                                                                                                    >  2:312128/1‑140 (MQ=255)
     gcgcAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAACTGGCATCAGGGCTTTCTTGTCCGGTCGGCTTCAAAAATGGCACCGACGGTACGATTAAAGTGGCTATCGATGCCATTAATGCCGCCGGTGcgccg                                                                                                                                    >  1:481083/1‑140 (MQ=255)
        cAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAACTGGCATCAGGGCTTTCTTGTCCGGTCGGCTTCAAAAATGGCACCGACGGTACGATTAAAGTGGCTATCGATGCCATTAATGCCGCCGGTGCGCCGCAc                                                                                                                                 <  1:326920/140‑1 (MQ=255)
        cAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAACTGGCATCAGGGCTTTCTTGTCCGGTCGGCTTCAAAAATGGCACCGACGGTACGATTAAAGTGGCTATCGATGCCATTAATGCCGCCGGTGCGCCGCAc                                                                                                                                 >  1:1679937/1‑140 (MQ=255)
              gcgcACGTACCACCGAATCGCAGGTGCACCGCGAACTGGCATCAGGGCTTTCTTGTCCGGTCGGCTTCAAAAATGGCACCGACGGTACGATTAAAGTGGCTATCGATGCCATTAATGCCGCCGGTGCGCCGCACTGCTTc                                                                                                                           <  2:650579/140‑1 (MQ=255)
                        caccGAATCGCAGGTGCACCGCGAACTGGCATCAGGGCTTTCTTGTCCGGTCGGCTTCAAAAATGGCACCGACGGTACGATTAAAGTGGCTATCGATGCCATTAATGCCGCCGGTGCGCCGCACTGCTTCCTGTCCGTaa                                                                                                                 <  1:1008053/140‑1 (MQ=255)
                        caccGAATCGCAGGTGCACCGCGAACTGGCATCAGGGCTTTCTTGTCCGGTCGGCTTCAAAAATGGCACCGACGGTACGATTAAAGTGGCTATCGATGCCATTAATGCCGCCGGTGCGCCGCACTGCTTCCTGTCCGTaa                                                                                                                 <  2:1305940/140‑1 (MQ=255)
                        caccGAATCGCAGGTGCACCGCGAACTGGCATCAGGGCTTTCTTGTCCGGTCGGCTTCAAAAATGGCACCGACGGTACGATTAAAGTGGCTATCGATGCCATTAATGCCGCCGGTGCGCCGCACTGCTTCCTGTCCGTaa                                                                                                                 <  2:399825/140‑1 (MQ=255)
                             aaTCGCAGGTGCACCGCGAACTGGCATCAGGGCTTTCTTGTCCGGTCGGCTTCAAAAATGGCACCGACGGTACGATTAAAGTGGCTATCGATGCCATTAATGCCGCTGGTGCGCCGCACTGCTTCCTGTCCGTAACGaaa                                                                                                            >  1:1306171/1‑140 (MQ=255)
                             aaTCGCAGGTGCACCGCGAACTGGCATCAGGGCTTTCTTGTCCGGTCGGCTTCAAAAATGGCACCGACGGTACGATTAAAGTGGCTATCGATGCCATTAATGCCGCTGGTGCGCCGCACTGCTTCCTGTCCGTAACGaaa                                                                                                            >  1:1305805/1‑140 (MQ=255)
                               tCGCAGGTGCACCGCGAACTGGCATCAGGGCTTTCTTGTCCGGTCGGCTTCAAAAATGGCACCGACGGTACGATTAAAGTGGCTATCGATGCCATTAATGCCGCCGGTGCGCCGCACTGCTTCCTGTCCGTAACGAAATg                                                                                                          >  1:356319/1‑140 (MQ=255)
                               tCGCAGGTGCACCGCGAACTGGCATCAGGGCTTTCTTGTCCGGTCGGCTTCAAAAATGGCACCGACGGTACGATTAAAGTGGCTATCGATGCCATTAATGCCGCCGGTGCGCCGCACTGCTTCCTGTCCGTAACGAAATg                                                                                                          >  1:520166/1‑140 (MQ=255)
                                     gTGCACCGCGAACTGGCATCAGGGCTTTCTTGTCCGGTCGGCTTCAAAAATGGCACCGACGGTACGATTAAAGTGGCTATCGATGCCATTAATGCCGCCGGTGCGCCGCACTGCTTCCTGTCCGTAACGAAATGGGGGCa                                                                                                    <  2:757373/140‑1 (MQ=255)
                                        cACCGCGAACTGGCATCAGGGCTTTCTTGTCCGGTCGGCTTCAAAAATGGCACCGACGGTACGATTAAAGTGGCTATCGATGCCATTAATGCCGCCGGTGCGCCGCACTGCTTCCTGTCCGTAACGAAATGGGGGCATTc                                                                                                 >  1:89317/1‑140 (MQ=255)
                                             cgAACTGGCATCAGGGCTTTCTTGTCCGGTCGGCTTCAAAAATGGCACCGACGGTACGATTAAAGTGGCTATCGATGCCATTAATGCCGCCGGTGCGCCGCACTGCTTCCTGTCCGTAACGAAATGGGGGCATTCGGCGa                                                                                            <  1:128271/140‑1 (MQ=255)
                                               aaCTGGCATCAGGGCTTTCTTGTCCGGTCGGCTTCAAAAATGGCACCGACGGTACGATTAAAGTGGCTATCGATGCCATTAATGCCGCCGGTgcg                                                                                                                                       >  2:48523/1‑95 (MQ=255)
                                               aaCTGGCATCAGGGCTTTCTTGTCCGGTCGGCTTCAAAAATGGCACCGACGGTACGATTAAAGTGGCTATCGATGCCATTAATGCCGCCGGTgcg                                                                                                                                       <  1:48523/95‑1 (MQ=255)
                                                     cATCAGGGCTTTCTTGTCCGGTCGGCTTCAAAAATGGCACCGACGGTACGATTAAAGTGGCTATCGATGCCATTAATGCCGCCGGTGCGCCGCACTGCTTCCTGTCCGTAACGAAATGGGGGCATTCGGCGATTGTGAAt                                                                                    >  2:913506/1‑140 (MQ=255)
                                                     cATCAGGGCTTTCTTGTCCGGTCGGCTTCAAAAATGGCACCGACGGTACGATTAAAGTGGCTATCGATGCCATTAATGCCGCCGGTGCGCCGCACTGCTTCCTGTCCGTAACGAAATGGGGGCATTCGGCGATTGTGAAt                                                                                    >  2:913964/1‑140 (MQ=255)
                                                     cATCAGGGCTTTCTTGTCCGGTCGGCTTCAAAAATGGCACCGACGGTACGATTAAAGTGGCTATCGATGCCATTAATGCCGCCGGTGCGCCGCACTGCTTCCTGTCCGTAACGAAATGGGGGCATTCGGCGATTGTGAAt                                                                                    >  1:1239676/1‑140 (MQ=255)
                                                        cAGGGCTTTCTTGTCCGGTCGGCTTCAAAAATGGCACCGACGGTACGATTAAAGTGGCTATCGATGCCATTAATGCCGCCGGTGCGCCGCACTGCTTCCTGTCCGTAACGAAATGGGGGCATTCGGCGATTGTGAATAcc                                                                                 <  2:1663385/140‑1 (MQ=255)
                                                          gggCTTTCTTGTCCGGTCGGCTTCAAAAATGGCACCGACGGTACGATTAAAGTGGCTATCGATGCCATTAATGCCGCCGGTGCGCCGCACTGCTTCCTGTCCGTAACGAAATGGGGGCATTCGGCGATTGTGAATACCAg                                                                               >  1:1594908/1‑140 (MQ=255)
                                                          cggCTTTCTTGTCCGGTCGGCTTCACAAATGGCACCGACGGTACGAATACAGTGGCTAACGAAGCCAGTAAAGCCGCCGGTGCGCCGCACGGCTTCCTGTCCGGAACGACATGGGGGCATTCCGCGAGTGTGAATACCAg                                                                               <  2:362146/139‑1 (MQ=255)
                                                                       cGGTCGGCTTCAAAAATGGCACCGACGGTACTATTAAAGTGGCTATCGATGCCATTAAAGCCGCCGCAGAGCCGCAATGCTACCTGTCCGTAACGAAATGTGGGCAAACTGCGATTGTCAATACCACCGGTAACGGCGAt                                                                  >  1:1590002/1‑140 (MQ=255)
                                                                       cGGTCGGCTTCAAAAATGGCACCGACGGTACGATTAAAGTGGCTATCGATGCCATTAATGCCGCCGGTGCGCCGCACTGCTTCCTGTCCGTAACGAAATGGGGGCATTCGGCGATTGTGAATACCAGCGGTAACGGCGAt                                                                  >  1:55586/1‑140 (MQ=255)
                                                                       cGGTCGGCTTCAAAAATGGCACCGACGGTACGATTAAAGTGGCTATCGATGCCATTAATGCCGCCGGTGCGCCGCACTGCTTCCTGTCCGTAACGAAATGGGGGCATTCGGCGATTGTGAATACCAGCGGTAACGGCGAt                                                                  >  1:58295/1‑140 (MQ=255)
                                                                       cGGTCGGCTTCAAAAATGGCACCGACGGTACGATTAAAGTGGCTATCGATGCCATTAATGCCGCCGGTGCGCCGCACTGCTTCCTGTCCGTAACGAAATGGGGGCATTCGGCGATTGTGAATACCAGCGGTAACGGCGAt                                                                  >  1:751764/1‑140 (MQ=255)
                                                                       cGGTCGGCTTCAAAAATGGCACCGACGGTACGATTAAAGTGGCTATCGATGCCATTAATGCCGCCGGTGCGCCGCACTGCTTCCTGTCCGTAACGAAATGGGGGCATTCGGCGATTGTGAATACCAGCGGTAACGGCGAt                                                                  >  1:751875/1‑140 (MQ=255)
                                                                         gTCGGCTTCAAAAATGGCACCGACGGTACGATTAAAGTGGCTATCGATGCCATTTATGCCGCCGGTGCGCCGCACTGCTTACTGTCCTTAACGAAATGGGGGCATTCGGCGATTGTTAATACCAGCGGTAACGGCGATTg                                                                <  1:1094298/140‑1 (MQ=255)
                                                                         gTCGGCTTCAAAAATGGCACCGACGGTACGATTAAAGTGGCTATCGATGCCATTAATGCCGCCGGTGCGCCGCACTGCTTCCTGTCCGTAACGAAATGGGGGCATTCGGCGATTGTGAATACCAGCGGTAACGGCGATTg                                                                <  1:1085597/140‑1 (MQ=255)
                                                                         gTCGGCTTCAAAAATGGCACCGACGGTACGATTAAAGTGGCTATCGATGCCATTAATGCCGCCGGTGCGCCGCACTGCTTCCTGTCCGTAACGAAATGGGGGCATTCGGCGATTGTGAATACCAGCGGTAACGGCGATTg                                                                <  1:1452444/140‑1 (MQ=255)
                                                                         gTCGGCTTCAAAAATGGCACCGACGGTACGATTAAAGTGGCTATCGATGCCATTAATGCCGCCGGTGCGCCGCACTGCTTCCTGTCCGTAACGAAATGGGGGCATTCGGCGATTGTGAATACCAGCGGTAACGGCGATTg                                                                <  1:1529628/140‑1 (MQ=255)
                                                                               ttCAAAAATGGCACCGACGGTACGATTAAAGTGGCTATCGATGCCATTAATGCCGCCGGTGCGCCGCACTGCTTCCTGTCCGTAACGAAATGGGGGCATTCGGCGATTGTGAATACCAGCGGTAACGGCGATTGCCatat                                                          >  2:1640613/1‑140 (MQ=255)
                                                                                 cAAAAATGGCACCGACGGTACGATTAAAGTGGCTATCGATGCCATTAATGCCGCCGGTGCGCCGCTCTGCTTCCTGTCCGTAACGAAATGGGGGCATTCGGCGATTGTGAATACCAGCGGTAACGGCGATTGCCATATCa                                                        <  2:142884/140‑1 (MQ=255)
                                                                                        ggCACCGACGGTACGATTAAAGTGGCTATCGATGCCATTAATGCCGCCGGTGCGCCGCACTGCTTCCTGTCCGTAACGAAATGGGGGCATTCGGCGATTGTGAATACCAGCGGTAACGGCGATTGCCATATCATTCTgcg                                                 <  1:1413815/140‑1 (MQ=255)
                                                                                        ggCACCGACGGTACGATTAAAGTGGCTATCGATGCCATTAATGCCGCCGGTGCGCCGCACTGCTTCCTGTCCGTAACGAAATGGGGGCATTCGGCGATTGTGAATACCAGCGGTAACGGCGATTGCCATATCATTCTgcg                                                 <  1:33441/140‑1 (MQ=255)
                                                                                                   tACGATTAACGTGGCTATCGATGCCATCTATGCCGCCGGAGCGCCGCACTGCTTCCTGTCCGAAACGACATGGGGGCATTCGGCGATTGTGAATACCAGCGGTAACGGCGATTGCCATATCATTCTGCGCGGCGGTAAag                                      <  2:1435322/140‑1 (MQ=255)
                                                                                                   tACGATTAAAGTGGCTATCGATGCCATTAATGCCGCCGGTGCGCCGCACTGCTTCCTGTCCGTAACGAAATGGGGGCATTCGGCGATTGTGAATACCAGCGGTAACGGCGATTGCCATATCATTCTGCGCGGCGGTAAag                                      <  2:186933/140‑1 (MQ=255)
                                                                                                   tACGATTAAAGTGGCTATCGATGCCATTAATGCCGCCGGTGCGCCGCACTGCTTCCTGTCCGTAACGAAATGGGGGCATTCGGCGATTGTGAATACCAGCGGTAACGGCGATTGCCATATCATTCTGCGCGGCGGTAAag                                      <  2:187042/140‑1 (MQ=255)
                                                                                                        ttAAAGTGGCTATCGTTGCCATTAATGCCGCCGGTGCGCCGCACTGCTTCCTGTCCGTAACGAAATGGGGGCATTCGGCGATTGTGAATACCAGCGGTAACGGCGATTGCCATATCATTCTGCGCGGCGGTAAAGAGCCt                                 <  2:1309631/140‑1 (MQ=255)
                                                                                                        ttAAAGTGGCTATCGATGCCATTAATGCCGCCGGTGCGCCGCACTGCTTCCTGTCCGTAACGAAATGGGGGCATTCGGCGATTGTGAATACCAGCGGTAACGGCGATTGCCATATCATTCTGCGCGGCGGTAAAGAGCCt                                 <  2:1309868/140‑1 (MQ=255)
                                                                                                        ttAAAGTGGCTATCGATGACATTAATGCCGCCGGTGCGCCGCACTGCTTCCTGTCCGTAACGAAATGGGGGCATTCGGCGATTGTGAATACCAGCGGTAACGGCGATTGCCATATCATTCTGCGCGGCGGTAAAGAGCCt                                 <  2:979080/140‑1 (MQ=255)
                                                                                                           aaGTGGCTATCGATGCCATTAATGCCGCCGGTGCGCCGCACTGCTTCCTGTCCGTAACGAAATGGGGGCATTCGGCGATTGTGAATACCAGCGGTAACGGCGATTGCCATATCATTCTGCGCGGCGGTAAAGAGCCTAAc                              >  1:595866/1‑140 (MQ=255)
                                                                                                             gTGGCTATCGATGCCATTAATGCCGCTGGTGCGCCGCACTGCTTCCTGTCCGTAACGAAATGGGGGCATTCGGCGATTGTGAATACCAGCGGTAACGGCGATTGCCATATCATTCTGCGCGGCGGTAAAGAGCCTAACTa                            <  1:534994/140‑1 (MQ=255)
                                                                                                             gTGGCTATCGATGCCATTAATGCCGCTGGTGCGCCGCACTGCTTCCTGTCCGTAACGAAATGGGGGCATTCGGCGATTGTGAATACCAGCGGTAACGGCGATTGCCATATCATTCTGCGCGGCGGTAAAGAGCCTAACTa                            <  1:534982/140‑1 (MQ=255)
                                                                                                                    tCGATGCCATTAATGCCGCCGGTGCGCCGCACTGCTTCCTGTCCGTAACGAAATGGGGGCATTCGGCGATTGTGAATACCAGCGGTAACGGCGATTGCCATATCATTCTGCGCGGCGGTAAAGAGCCTAACTACAGcgcg                     <  2:1181889/140‑1 (MQ=255)
                                                                                                                    tCGATGCCATTAATGCCGCCGGTGCGCCGCACTGCTTCCTGTCCGTAACGAAATGGGGGCATTCGGCGATTGTGAATACCAGCGGTAACGGCGATTGCCATATCATTCTGCGCGGCGGTAAAGAGCCTAACTACAGcgcg                     <  2:416967/140‑1 (MQ=255)
                                                                                                                    tCGATGCCATTAATGCCGCCGGTGCGCCGCACTGCTTCCTGTCCGTAACGAAATGGGGGCATTCGGCGATTGTGAATACCAGCGGTAACGGCGATTGCCATATCATTCTGCGCGGCGGTAAAGAGCCTAACTACAGcgcg                     <  2:1420592/140‑1 (MQ=255)
                                                                                                                    tCGATGCCATTAATGCCGCCGGTGCGCCGCACTGCTTCCTGTCCGTAACGAAATGGGGGCATTCGGCGATTGTGAATACCAGCGGTAACGGCGATTGCCATATCATTCTGCGCGGCGGTAAAGAGCCTAACTACAGcgcg                     <  2:1420492/140‑1 (MQ=255)
                                                                                                                    tCGATGCCATTAATGCCGCCGGTGCGCCGCACTGCTTCCTGTCCGTAACGAAATGGGGGCATTCGGCGATTGTGAATACCAGCGGTAACGGCGATTGCCATATCATTCTGCGCGGCGGTAAAGAGCCTAACTACAGcgcg                     <  2:1420141/140‑1 (MQ=255)
                                                                                                                    tCGATGCCATTAATGCCGCCGGTGCGCCGCACTGCTTCCTGTCCGTAACGAAATGGGGGCATTCGGCGATTGTGAATACCAGCGGTAACGGCGATTGCCATATCATTATGCGCGGCGGTAAAGAGCCTAACTACAGcgcg                     <  2:20800/140‑1 (MQ=255)
                                                                                                                      gATGCCATTAATGCCGCCGGTGCGCCGCACTGCTTCCTGTCCGTAACGAAATGGGGGCATTCGGCGATTGTGAATACCAGCGGTAACGGCGATTGCCATATCATTCTGCGCGGCGGTAAAGAGCCTAACTACAGCGCGaa                   <  2:71071/140‑1 (MQ=255)
                                                                                                                       aTGCCATTAATGCCGCCGGTGCGCCGCACTGCTTCCTGTCCGTAACGAAATGGGGGCATTCGGCGATTGTGAATACCAGCGGTAACGGCGATTGCCATATCATTCTGCGCGGCGGTAAAGAGCCTAACTACAGCGCGAAg                  >  1:1600642/1‑140 (MQ=255)
                                                                                                                                 tGCCGCCGGTGCGCCGCACTGCTTCCTGTCCGTAACGAAATGGGGGCATTCGGCGATTGTGAATACCAGCGGTAACGGCGATTGCCATATCATTCTGCGCGGCGGTAAAGAGCCTAACTACAGCGCGAAGCACGTTGCTg        <  1:1314787/140‑1 (MQ=255)
                                                                                                                                     gccgGTGCGCCGCACTGCTTCCTGTCCGTAACGAAATGGGGGCATTCGGCGATTGTGAATACCAGCGGTAACGGCGATTGCCATATCATTCTGCGCGGCGGTAAAGAGCCTAACTACAGCGCGAAGCACGTTGCTGAAGt    <  2:675651/140‑1 (MQ=255)
                                                                                                                                     gccgGTGCGCCGCACTGCTTCCTGTCCGTAACGAAATGGGGGCATTCGGCGATTGTGAATACCAGCGGTAACGGCGATTGCCATATCATTCTGCGCGGCGGTAAAGAGCCTAACTACAGCGCGAAGCACGTTGCTGAAGt    <  2:675676/140‑1 (MQ=255)
                                                                                                                                     gccgGTGCGCCGCACTGCTTCCTGTCCGTAACGAAATGGGGGCATTCGGCGATTGTGAATACCAGCGGTAACGGCGATTGCCATATCATTCTGCGCGGCGGTAAAGAGCCTAACTACAGCGCGAAGCACGTTGCTGAAGt    <  1:1420781/140‑1 (MQ=255)
                                                                                                                                       cgGTGCGCCGCACTGCTTCCTGTCCGTAACGAAATGGGGGCATTCGGCGATTGTGAATACCAGCGGTAACGGCGATTGCCATATCATTCTGCGCGGCGGTAAAGAGCCTAACTACAGCGCGAAGCACGTTGCTGAAGTGa  >  2:1428892/1‑140 (MQ=255)
                                                                                                                                       |                                                                                                                                           
CTGGGGCGCAATTGGCGCACGTACCACCGAATCGCAGGTGCACCGCGAACTGGCATCAGGGCTTTCTTGTCCGGTCGGCTTCAAAAATGGCACCGACGGTACGATTAAAGTGGCTATCGATGCCATTAATGCCGCCGGTGCGCCGCACTGCTTCCTGTCCGTAACGAAATGGGGGCATTCGGCGATTGTGAATACCAGCGGTAACGGCGATTGCCATATCATTCTGCGCGGCGGTAAAGAGCCTAACTACAGCGCGAAGCACGTTGCTGAAGTGA  >  NC_000913/786106‑786380

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