Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 959,667 A→C 7.5% T286P (ACC→CCC)  aroA → 5‑enolpyruvylshikimate‑3‑phosphate synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913959,6670AC 7.5% 154.9 / 2.9 62T286P (ACC→CCC) aroA5‑enolpyruvylshikimate‑3‑phosphate synthetase
Reads supporting (aligned to +/- strand):  ref base A (30/27);  new base C (2/3);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

TTATTTGGTCGAAGGCGATGCATCTTCGGCTTCTTACTTTCTGGCAGCAGCAGCAATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAATAT  >  NC_000913/959519‑959812
                                                                                                                                                    |                                                                                                                                                 
ttATTTGGTCGAAGGCGATGCATCTTCGGCTTCTTACTTTCTGGCAGCAGCAGCAATCAATGGCGGCACTTTAAAAGTGACCGGTATTGGACGTAACAGTAAGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGa                                                                                                                                                   <  1:789448/149‑1 (MQ=255)
  aTTTGGTCGAAGGCGATGCATCTTCGGCTTCTTACTTTCTGGCAGCAGCAGCAATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGTGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGAcc                                                                                                                                                 >  1:629044/1‑149 (MQ=255)
  aTTTGGTCGAAGGCGATGCATCTTCGGCTTCTTACTTTCTGGCAGCAGCAGCAATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGAcc                                                                                                                                                 >  2:1120802/1‑149 (MQ=255)
       gTCGAAGGCGATGCATCTTCGGCTTCTTACTTTCTGGCAGCAGCAGCAATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTg                                                                                                                                            >  2:140624/1‑149 (MQ=255)
       gTCGAAGGCGATGCATATTCGGCTTCTTACTTTCTGGCAGCAGCAGCCATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACTGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCAtttt                                                                                                                                            >  2:1099145/1‑148 (MQ=255)
        tCGAAGGCGATGCATCTTCGGCTTCTTACTTTCTGGCAGCAGCAGCAATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGc                                                                                                                                           <  2:552045/149‑1 (MQ=255)
         cGAAGGCGATGCATCTTCGGCTTCTTACTTTCTGGCAGCAGCAGCAATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGTCGCGCCCAtttttt                                                                                                                                          <  1:538414/149‑4 (MQ=255)
           aaGGCGATGCAGCTTCGGCTTCTTACTTTCTGGCAGCAGCAGCAATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGTCGCGCCCAttttttcc                                                                                                                                        <  1:222586/149‑6 (MQ=255)
             ggCGAAGCATCTACGGCTGCTTACTATCTGGAAGTAGCAGCTAGCAACGGCGGCAATGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGGGATATTCGCTTTGCTGATGTGCTGGAAAATATGGTCGCGCCCAttttttcgag                                                                                                                                      <  2:21035/149‑8 (MQ=255)
               cGATGCATCTTCGGCTTCTTACTTTCTGGCAGCAGAAGCAATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCTGGGTGATATTCGCTTTTCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCg                                                                                                                                    >  1:1101505/1‑149 (MQ=255)
                gATGCATCTTCGGCTTCTTACTTTCTGGCAGCAGCAGCAATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCga                                                                                                                                   >  1:363642/1‑149 (MQ=255)
                   gCATCTTCGGCTTCTTACTTTCTGGCAGCAGCAGCAATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCgatga                                                                                                                                >  1:468225/1‑149 (MQ=255)
                            gCTTCTTACTTTCTTGCATCATCAGCAATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTc                                                                                                                       <  2:792142/149‑1 (MQ=255)
                              ttcttACTTTCTGGCAGCAGCAGCAATCAAAGGCGGCACTGTAAATGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCt                                                                                                                     <  2:507591/149‑1 (MQ=255)
                              ttcttACTTTCTGGCAGCAGCAGCAATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCt                                                                                                                     <  1:523044/149‑1 (MQ=255)
                               tcttACTTTCTGGCAGCAGCAGCAATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTg                                                                                                                    <  1:965294/149‑1 (MQ=255)
                                cttACTTTCTGGCAGCAGCAGCAATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGc                                                                                                                   <  1:140624/149‑1 (MQ=255)
                                     tttCTGGCAGCAGCAGCAATCAAAGGCGGCACTGTAAAAGTGACCGGGATTGGACGTAACAGTATGCAGGGTGATATTCTCTTTGCTGATGTGCTGGATAAAATTGGCGCGACCATTTGCTGAGGCGATGATTATATTTCCTGCAcgag                                                                                                              >  2:784697/1‑147 (MQ=255)
                                         tGGCAGCAGCAGCAATCAAAGGCGGCACTGTAAATGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCgtggt                                                                                                          <  1:488780/149‑1 (MQ=255)
                                         tGGCAGCAGCAGCAATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCGTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCgtggt                                                                                                          >  2:988910/1‑149 (MQ=255)
                                         tGGCAGCAGCAGCAATCAAAGGCGGCACTGTAAAAGTGACCGGGAGTGGACGTAACAGTATGCAGGGTGATATTCGCTTGGCTGATGGGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCgtggt                                                                                                          <  2:236083/149‑1 (MQ=255)
                                          ggCAGCAGCAGCAATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCgtggtg                                                                                                         <  1:1120802/149‑1 (MQ=255)
                                          ggCAGCAGCAGCAATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCgtggtg                                                                                                         <  1:1071276/149‑1 (MQ=255)
                                           gcagcagcagcaATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGtga                                                                                                        >  1:24050/1‑149 (MQ=255)
                                            cagcagcagcaATCAAAGGCGGCACTGTAAAAGTGACCGGTATGGGACGTAACAGTATGCAGGGTGATATTCGCTTGGCTGAGGTGCTGGAAAAAATGCGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGtgaa                                                                                                       <  2:207885/149‑1 (MQ=255)
                                              gcagcagcaATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTTGCTTTGCTGATGTGCTGGAAAAAATGTGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGttaatt                                                                                                     <  1:1022390/149‑3 (MQ=255)
                                              gcagcagcaATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGtgaact                                                                                                     <  1:391051/149‑1 (MQ=255)
                                                 gcagcaATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGtgaactgaa                                                                                                  <  2:19279/149‑1 (MQ=255)
                                                 gcagcaATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGtgaactgaa                                                                                                  >  2:92646/1‑149 (MQ=255)
                                                       aTCAAATGCGGCACTGTAAAAGGGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTAt                                                                                            <  2:1065514/149‑1 (MQ=255)
                                                         cAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTg                                                                                          >  2:406618/1‑149 (MQ=255)
                                                             ggcgggACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTgatat                                                                                      <  1:326320/149‑1 (MQ=255)
                                                             ggcggcACTGTAACAGGGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTgatat                                                                                      <  1:321693/149‑1 (MQ=255)
                                                                     tGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGa                                                                              >  1:556784/1‑149 (MQ=255)
                                                                      gTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCCGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGaa                                                                             <  2:16140/149‑1 (MQ=255)
                                                                       tAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGCACTGAACGCTATTGATATGGATATGAAc                                                                            >  1:410543/1‑149 (MQ=255)
                                                                        aataGTGACCGGTATTGGACGTAACAGTATGCAGGGGGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAAcc                                                                           <  2:263294/146‑1 (MQ=255)
                                                                                    tATTGTACGTAACAGTATGTAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATg                                                               <  1:383814/149‑1 (MQ=255)
                                                                                         gACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATAGATATGAACCATATTCCTGATgcggcg                                                          >  2:998246/1‑149 (MQ=255)
                                                                                               aCAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGAcc                                                    >  2:612409/1‑149 (MQ=255)
                                                                                                cAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGACCa                                                   >  2:849334/1‑149 (MQ=255)
                                                                                                 cgTATGCAGGGGGATATGCGCTTTGCTGATGGGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATAATTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGACCAt                                                  <  1:857749/148‑1 (MQ=255)
                                                                                                   tATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGACCATTg                                                >  2:727794/1‑149 (MQ=255)
                                                                                                     tGCAGGGTGATATTCGCATTGATGAAGCGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATTAACCATATTCCTGATGCGGCGATGACCATTGcc                                              <  1:1070857/149‑1 (MQ=255)
                                                                                                        aGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGCCCATTGCCAcg                                           >  1:58859/1‑149 (MQ=255)
                                                                                                          ggTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGACCATTGCCAcggc                                         >  2:177004/1‑149 (MQ=255)
                                                                                                           gTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGCCCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGACCATTGCCTcggcg                                        >  1:330722/1‑149 (MQ=255)
                                                                                                             gATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGCACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGACCATTGCCACggcggc                                      >  2:411567/1‑149 (MQ=255)
                                                                                                                atTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGCGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATCGATATGAACCATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGtt                                   >  2:186636/1‑149 (MQ=255)
                                                                                                                 ttCGCTTTGCTGATGTTCTGGAAAAAATCGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTCATGCGGCGATGACCATTGCCACGGCGGCGTTa                                  >  2:103139/1‑149 (MQ=255)
                                                                                                                    gCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGCGGCGATGATTATATTTCCTGCACGCGTGGTGCACTGAACGCTATTGATATGGAGATGAACCATATTCCTGATGAGGCGATCACCATTGACACGGCGGAGTTAttt                               >  1:645400/1‑149 (MQ=255)
                                                                                                                          cTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCaaaa                         <  1:988910/149‑1 (MQ=255)
                                                                                                                            gATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGTTATGGATATGAACCATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAgg                       <  2:363642/149‑1 (MQ=255)
                                                                                                                            gATGTGCAGGAATAAAAGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGCATATGAACCATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAgg                       <  2:363420/149‑1 (MQ=255)
                                                                                                                             aTGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGc                      >  1:34213/1‑149 (MQ=255)
                                                                                                                                 gCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCTCGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGAGGACCATTGCCACGGCGGCGTTATTTGCTAAAGGcaccc                  >  1:1149619/1‑148 (MQ=255)
                                                                                                                                  cTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGcaccac                 <  1:406618/149‑1 (MQ=255)
                                                                                                                                    ggAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCAAGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGACCATTGACACGGCGGCGTTATTTGCAAAAGTaaccacca               >  1:18215/1‑149 (MQ=255)
                                                                                                                                      agcaaaTGGTCGCGCCCATTTGCGGGGGGGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGACCCTTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACg             >  2:1077077/4‑149 (MQ=255)
                                                                                                                                       aaaaaTGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGCATATGAACCATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACGc            >  2:1122775/1‑149 (MQ=255)
                                                                                                                                         aaaTGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGAACATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACGCTg          <  2:24050/149‑1 (MQ=255)
                                                                                                                                          aaTGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTGCTGATGCGGCGATGACCATTGCCACGGCGGCGTGATTTGCAAAAGGCACCACCACGCTgc         >  1:219934/1‑149 (MQ=255)
                                                                                                                                                 gcgACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAatat  >  2:447125/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TTATTTGGTCGAAGGCGATGCATCTTCGGCTTCTTACTTTCTGGCAGCAGCAGCAATCAAAGGCGGCACTGTAAAAGTGACCGGTATTGGACGTAACAGTATGCAGGGTGATATTCGCTTTGCTGATGTGCTGGAAAAAATGGGCGCGACCATTTGCTGGGGCGATGATTATATTTCCTGCACGCGTGGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAATAT  >  NC_000913/959519‑959812

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑