Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,262,293 T→C 5.6% P195P (CCA→CCG)  ispE ← 4‑diphosphocytidyl‑2‑C‑methylerythritol kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,262,2930TC 5.6% 235.9 / 3.3 71P195P (CCA→CCG) ispE4‑diphosphocytidyl‑2‑C‑methylerythritol kinase
Reads supporting (aligned to +/- strand):  ref base T (31/36);  new base C (2/2);  total (33/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

AGACACAGGCCCCTGTCCCAGTCAGGCGCGACGGGGCGTATTCTAACAGCCAGGAAAGCACCGCATCAACCTCGCGAAAACGTTTTCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGAATACTTACACCAGGGTGCGCCACCAGATACCACTTCTCTGGCGGATCCACCGGCGTTAGTATTTCACCAACGCCTTCGGCAAACGCGGCATGCCCCCGA  >  NC_000913/1262145‑1262440
                                                                                                                                                    |                                                                                                                                                   
aGACACAGGCCCCTGTCCCAGTCAGGCGCGACGGGGCGTATTCTAACAGCCAGGAAAGCACCGCATCAACCTCGCGAAAACGTTTTCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCttttt                                                                                                                                                     <  1:103891/149‑1 (MQ=255)
aGACACAGGCCCCTGTCCCAGTCAGGCGCGACGGGGCGTATTCTAACAGCCAGGAAAGCACCGCATCAACCTCGCGAAAACGTTTTCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCttttt                                                                                                                                                     <  2:547075/149‑1 (MQ=255)
  acacAGGCCCCTGTCCCAGTCAGGCGCGACGGGGCGTATTCTAACAGCCAGGAAAGCACCGCATCAACCTCGCGAAAACGTTTTCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTgg                                                                                                                                                   >  3:171248/1‑149 (MQ=255)
  acacAGGCCCCTGTCCCAGTCAGGCGCGACGGGGCGTATTCTAACAGCCAGGAAAGCACCGCATCAACCTCGCGAAAACGTTTTCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTgg                                                                                                                                                   >  2:903551/1‑149 (MQ=255)
  acacAGGCCCCTGTCCCAGTCAGGCGCGACGGGGCGTATTCTAACAGCCAGGAAAGCACCGCATCAACCTCGCGAAAACGTTTTCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTgg                                                                                                                                                   >  2:731637/1‑149 (MQ=255)
      aGGCCCCTGTCCCAGTCAGGCGCGACGGGGCGTATTCTAACAGCCAGGAAAGCACCGCATCAACCTCGCGAAAACGTTTTCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTa                                                                                                                                               >  2:699481/1‑149 (MQ=255)
           ccTGTCCCAGTCAGGCGCGACGGGGCGTATTCTAACAGCCAGGAAAGCACCGCATCAACCTCGCGAAAACGTTTTCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTCGGCGTATTgcg                                                                                                                                          >  1:150592/1‑149 (MQ=255)
            cTGTCCCAGTCAGGCGCGACGGGGCGTATTCTAACAGCCAGGAAAGCACCGCATCAACCTCGCGAAAACGTTTTCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTgcgc                                                                                                                                         >  4:138322/1‑149 (MQ=255)
               tCCCAGTCAGGCGCGACGGGGCGTATTCTAACAGCCAGGAAAGCACCGAATCAACCTCGCGAAAACGTTTTCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGATTCTATTGCCATTTTTGGCGTATTCCGCggg                                                                                                                                      >  4:193666/1‑149 (MQ=255)
                cccAGTCAGGCGCGACGGGGCGTATTCTAACAGCCAGGAAAGCACCGCATCAACCTCGCGAAAACGTTTTCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTGCGCGGGa                                                                                                                                     <  1:108179/149‑1 (MQ=255)
                cccAGTCAGGCGCGACGGGGCGTATTCTAACAGCCAGGAAAGCACCGCATCAACCTCGCGAAAACGTTTTCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTCGGCGTATTGCGCGGGa                                                                                                                                     <  2:150592/149‑1 (MQ=255)
                     tCAGGCGCGACGGGGCGTATTCTAACAGCCAGGAAAGCACCGCATCAACCTCGCGAAAACGTTTTCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCa                                                                                                                                <  1:238435/149‑1 (MQ=255)
                        ggCGCGACGGGGCGTATTCTAACAGCCAGGAAAGCACCGCATCAACCTCGCGAAAACGTTTTCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGa                                                                                                                             <  3:43925/149‑1 (MQ=255)
                        ggCGCGACGGGGCGTATTCTAACAGCCAGGAAAGCACCGCATCAACCTCGCGAAAACGTTTTCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGa                                                                                                                             <  4:476080/149‑1 (MQ=255)
                          cgcgACGGGGCGTATTCTAACAGCCAGGAAAGCACCGCATCAACCTCGCGAAAACGTTTTCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATc                                                                                                                           <  2:206934/149‑1 (MQ=255)
                                      tATTCTAACAGCCAGGAAAGCACCGCATCAACCTCGCGAAAACGTTTTCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCAc                                                                                                               <  2:379702/149‑1 (MQ=255)
                                          cTAACAGCCAGGAAAGCACCGCATCAACCTCGCGAAAACGTTTTCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGa                                                                                                           <  1:850629/149‑1 (MQ=255)
                                              cAGCCAGGAAAGCACCGCATCAACCTCGCGAAAACGTTTTCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCg                                                                                                       <  2:836611/149‑1 (MQ=255)
                                              cAGCCAGGAAAGCACCGCATCAACCTCGCGAAAACGTTTTCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCg                                                                                                       >  2:827342/1‑149 (MQ=255)
                                               aGCCAGGAAAGCACCGCATCAACCTCGCGAAAACGTTTTCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCgg                                                                                                      >  2:527541/1‑149 (MQ=255)
                                               aGCCAGGAAAGCACCGCATCAACCTCGCGAAAACGTTTTCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTGCACGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCgg                                                                                                      >  4:625857/1‑149 (MQ=255)
                                                   aGGAAAGCACCGCATCAACCTCGCGAAAACGTTTTCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGAATa                                                                                                  >  4:381353/1‑149 (MQ=255)
                                                    ggAAAGCACCGCATCAACCTCGCGAAAACGTTTTCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGAATAc                                                                                                 >  3:324787/1‑149 (MQ=255)
                                                    ggAAAGCACCGCATCAACCTCGCGAAAACGTTTTCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGAATAc                                                                                                 <  2:456524/149‑1 (MQ=255)
                                                     gAAAGCACCGCATCAACCTCGCGAAAACGTTTTCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGa                                                                                                     >  2:508676/1‑144 (MQ=255)
                                                     gAAAGCACCGCATCAACCTCGCGAAAACGTTTTCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGa                                                                                                     <  1:508676/144‑1 (MQ=255)
                                                     gAAAGCACCGCATCAACCTCGCGAAAACGTTTTCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTCGGCGTATTGCGCGGGAGTTCAGGa                                                                                                                             >  2:608789/1‑120 (MQ=255)
                                                     gAAAGCACCGCATCAACCTCGCGAAAACGTTTTCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTCGGCGTATTGCGCGGGAGTTCAGGa                                                                                                                             <  1:608789/120‑1 (MQ=255)
                                                       aaGCACCGCATCAACCTCGCGAAAACGTTTTCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGAATACTTa                                                                                              <  1:731637/149‑1 (MQ=255)
                                                        aGCACCGCATCAACCTCGCGAAAACGTTTTCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGAATACTTac                                                                                             <  2:747726/149‑1 (MQ=255)
                                                            ccGCATCAACCTCGCGAAAACGTTTTCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGAATACTTACACCa                                                                                         <  4:29259/149‑1 (MQ=255)
                                                             cGCATCAACCTCGCGAAAACGTTTTCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGAATACTTACACCAg                                                                                        <  1:727183/149‑1 (MQ=255)
                                                              gCATCAACCTCGCGAAAACGTTTTCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGCGTCGGAATACTTACACCAgg                                                                                       <  2:468010/149‑1 (MQ=255)
                                                                  cAACCTCGCGAAAACGTTTTCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGAATACTTACACCAGGGTgc                                                                                   <  3:809521/149‑1 (MQ=255)
                                                                  cAACCTCGCGAAAACGTTTTCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGAATACTTACACCAGGGTgc                                                                                   <  4:534804/149‑1 (MQ=255)
                                                                      cTCGCGAAAACGTTTTCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTAt                                                                                                                                              >  1:126534/1‑86 (MQ=255)
                                                                      cTCGCGAAAACGTTTTCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTAt                                                                                                                                              <  2:126534/86‑1 (MQ=255)
                                                                       tCGCGAAAACGTTTTCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGAATACTTACACCAGGGTGCGccac                                                                              >  4:208508/1‑149 (MQ=255)
                                                                            aaaaCGTTTTCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGAATACTTACACCAGGGTGCGCCACCAGAt                                                                         <  2:808846/149‑1 (MQ=255)
                                                                            aaaaCGTTTTCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGAATACTTACACCAGGGTGCGCCACCAGAt                                                                         <  1:675718/149‑1 (MQ=255)
                                                                             aaaCGTTTTCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGAATACTTACACCAGGGTGCTCCACCAGATa                                                                        >  4:493897/1‑149 (MQ=255)
                                                                             aaaCGTTTTCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGAATACTTACACCAGGGTGCGCCACCAGAtg                                                                        >  1:827370/1‑148 (MQ=255)
                                                                              aaCGTTTTCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGAATACTTACACCAGGGTGCGCCACCAGATAc                                                                       >  1:361764/1‑149 (MQ=255)
                                                                              aaCGTTTTCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGAATACTTACACCAGGGTGCGCCACCAGATAc                                                                       <  4:728367/149‑1 (MQ=255)
                                                                              aaCGTTTTCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGAATACTTACACCAGGGTGCGCCACCAGATAc                                                                       <  1:527541/149‑1 (MQ=255)
                                                                                     tCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGAATACTTACACCAGGGTGCGCCACCAGATACCACTtct                                                                <  3:625857/149‑1 (MQ=255)
                                                                                           gATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGAATACTTACACCAGGGTGCGCCACCAGATACCACTTCTCTGGCg                                                          <  2:549856/149‑1 (MQ=255)
                                                                                             tAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGAATACTTACACCAGGGTGCGACACCAGATACCACTTCTCTGGCGGa                                                        <  2:606425/149‑1 (MQ=255)
                                                                                               aCCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGAATACTTACACCAGGGTGCGCCACCAGATACCACTTCTCTGGCGGATc                                                      >  3:96579/1‑149 (MQ=255)
                                                                                                  tCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGAATACTTACACCAGGGTGCGCCACCAGATACCACTTCTCTGGCGGATCCAc                                                   >  2:686907/1‑149 (MQ=255)
                                                                                                  tCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAAGATCTTTAAAAATCACCGGAGTCGGAATACTTACACCAGGGTGCGCCACCAGATACCACTTCTCTGGCGGATCCAc                                                   >  2:868478/1‑149 (MQ=255)
                                                                                                       aTCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGAATACTTACACCAGGGTGCGCCACCAGATACCACTTCTCTGGCGGATCCACCGGCg                                              >  2:851054/1‑149 (MQ=255)
                                                                                                          aTTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGAATACTTACACCAGGGTGCGCCACCAGATACCACTTCTCTGGCGGATCCACCGGCGTTa                                           <  1:699481/149‑1 (MQ=255)
                                                                                                           ttGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGAATACTTACACCAGGGTGCGCCACCAGATACCACTTCTCTGGCGGATCCACCGGCGTTAg                                          <  3:138322/149‑1 (MQ=255)
                                                                                                              cTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGAATACTTACACCAGGGTGCGCCACCAGATACCACTTCTCTGGCGGATCCACCGGCGTTAGTAt                                       <  4:96579/149‑1 (MQ=255)
                                                                                                              cTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGAATACTTACACCAGGGTGCGCCACCAGATACCACTTCTCTGGCGGATCCACCGGCGTTAGTAt                                       <  3:412538/149‑1 (MQ=255)
                                                                                                                gAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATGTTTAAAAATCACCGGAGTCGGAATACTTACACCAGGGTGCGTCACCAGATACCACTTCTCTGGCGGATCCACCGGCGTTAGTAttt                                     <  4:568070/149‑1 (MQ=255)
                                                                                                                 aaTTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGAATACTTACACCAGGGTGCGCCACCAGATACCACTTCTCTGGCGGATCCACCGGCGTTAGTATTTc                                    >  2:936893/1‑149 (MQ=255)
                                                                                                                     cacaTTTTAGCAACGTTTCTATTGACCTTTTTGTCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGAATACTTACACCAGGGTGCGCCACCAGATACCACTTCTCTGGCGGATCCACCGGCGTGAGTATTTCACCa                                >  1:303925/1‑149 (MQ=255)
                                                                                                                      acaTTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGAATACTTACACCAGGGTGCGCCACCAGATACCACTTCTCTGGCGGATCCACCGGCGTTAGTATTTCACCaa                               <  1:936893/149‑1 (MQ=255)
                                                                                                                        aTTTTAGCAACGTTTCTATTGAACTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGAATACTTACACCAGGGTGCGCCACCAGATACCACTTCTCTGGCGGATCCACCGGCGTTAGTATTTCACCAACg                             >  1:568580/1‑149 (MQ=255)
                                                                                                                                 aCGTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGAATACTTACACCAGGGTGCGCCACCAGATACCACTTCTCTGGCGGATCCACCGGCGTTAGTATTTCACCAACGCCTTCGGCa                    >  3:429817/1‑149 (MQ=255)
                                                                                                                                   gTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGAATACTTACACCAGGGTGCGCCACCAGa                                                                          <  1:855910/93‑1 (MQ=255)
                                                                                                                                   gTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGAATACTTACACCAGGGTGCGCCACCAGa                                                                          >  2:855910/1‑93 (MQ=255)
                                                                                                                                   gTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGAATACTTACACCAGGGTGCGCCACCAGATACCACTTCTCTGGCGGATCCACCGGCGTTAGTATTTCACCAACGCCTTCGGCaaa                  >  4:69743/1‑149 (MQ=255)
                                                                                                                                    tttCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGAATACTTACACCAGGGTGCGCCACCAGATACCACTTCTCTGGCGGATCCACCGGCGTTAGTATTTCACCAACGCCTTCGGCAAAc                 >  4:907577/1‑149 (MQ=255)
                                                                                                                                    tttCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGAATACTTACACCAGGGTGCGCCACCAGATACCACTTCTCTGGCGGATCCACCGGCGTTAGTATTTCACCAACGCCTTCGGCAAAc                 <  3:315675/149‑1 (MQ=255)
                                                                                                                                    tttCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGAATACTTACACCAGGGTGCGCCACCAGATACCACTTCTCTGGCGGATCCACCGGCGTTAGTATTTCACCAACGCCTTCGGCAAAc                 <  1:40229/149‑1 (MQ=255)
                                                                                                                                      tCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGAATACTTACACCAGGGTGCGCCACCAGATACCACTTCTCTGGCGGATCCACCGGCGTTAGTATTTCACCAACGCCTTCGGCAAAcgc               <  2:450212/149‑1 (MQ=255)
                                                                                                                                          ttGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGAATACTTACACCAGGGTGCGCCACCAGATACCACTTCTCTGGCGGATCCACCGGCGTTAGTATTTCACCAACGCCTTCGGCAAACGCGGCa           >  1:642987/1‑149 (MQ=255)
                                                                                                                                             aCCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGAATACTTACACCAGGGTGCGCCACCAGATACCACTTCTCTGGCGGATCCACCgg                                                <  1:645019/109‑1 (MQ=255)
                                                                                                                                             aCCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGAATACTTACACCAGGGTGCGCCACCAGATACCACTTCTCTGGCGGATCCACCgg                                                >  2:645019/1‑109 (MQ=255)
                                                                                                                                               cTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGAATACTTACACCAGGGTGCGCCACCAGATACCACTTCTCTGGCGGATCCACCGGCGTTAGTATTTCACCAACGCCTTCGGCAAACGCGGCATGccc      >  2:883005/1‑149 (MQ=255)
                                                                                                                                                 ttttGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGAATACTTACACCAGGGTGCGCCACCAGATACCACTTCTCTGGCGGATCCACCGGCGTTAGTATTTCACCAACGCCTTCGGCAAACGCGGCATGccccc    <  3:193666/149‑1 (MQ=255)
                                                                                                                                                   ttGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGAATACTTACACCAGGGTGCGCCACCAGATACCACTTCTCTGGCGGATCCACCGGCGTTAGTATTTCACCAACGCCTTCGGCAAACGCGGCATGCCCCCGa  >  2:128318/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AGACACAGGCCCCTGTCCCAGTCAGGCGCGACGGGGCGTATTCTAACAGCCAGGAAAGCACCGCATCAACCTCGCGAAAACGTTTTCTTGCGATAACCTCGCAATCATTGCTGAATTCACATTTTAGCAACGTTTCTATTGACCTTTTTGGCGTATTGCGCGGGAGTTCAGGATCTTTAAAAATCACCGGAGTCGGAATACTTACACCAGGGTGCGCCACCAGATACCACTTCTCTGGCGGATCCACCGGCGTTAGTATTTCACCAACGCCTTCGGCAAACGCGGCATGCCCCCGA  >  NC_000913/1262145‑1262440

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: ‑