Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,431,978 T→C 5.3% K322R (AAA→AGA)  folC ← bifunctional folylpolyglutamate synthase/ dihydrofolate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,431,9780TC 5.3% 231.0 / 2.9 75K322R (AAA→AGA) folCbifunctional folylpolyglutamate synthase/ dihydrofolate synthase
Reads supporting (aligned to +/- strand):  ref base T (37/34);  new base C (2/2);  total (39/36)
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.08e-01

TGCTCAAGCAGTTGTTCTGCCGTGGCACCGCGCGGCCCTTCCAGTGGCGCACAATACCAGTCATCAACCACGCTTTTCAACCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTGAGATATTCCGCCGCATGTGGATTATGCGCGACATCAAAAATAACGCGTGGCGACTCGCTCACAATCTGGAAACGTCCCGGCAAAATTGCGCTGGCAATCCCGTCGCGAATGGCATTTTCA  >  NC_000913/2431831‑2432124
                                                                                                                                                   |                                                                                                                                                  
tGCTCAAGCAGTTGTTCTGCCGTGGCACCGCGCGGCCCTTCCAGTGGCGCACAATACCAGTCATCAACCACGCTTTTCAACCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGttttt                                                                                                                                                   <  5:35692/149‑1 (MQ=255)
tGCTCAAGCAGTTGTTCTGCCGTGGCACCGCGCGGCCCTTCCAGTGGCGCACAATACCAGTCATCAACCACGCTTTTCAACCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGttttt                                                                                                                                                   <  7:266561/149‑1 (MQ=255)
  cTCAAGCAGTTGTTCTGCCGTGGCACCGCGCGGCCCTTCCAGTGGCGCACAATACCAGTCATCAACCACGCTTTTCAACCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTc                                                                                                                                                  <  7:3271/148‑1 (MQ=255)
  cTCAAGCAGTTGTTCTGCCGTGGCACCGCGCGGCCCTTCCAGTGGCGCACAATACCAGTCATCAACCACGCTTTTCAACCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTc                                                                                                                                                  <  1:169121/148‑1 (MQ=255)
  cTCAAGCAGTTGTTCTGCCGTGGCACCGCGCGGCCCTTCCAGTGGCGCACAATACCAGTCATCAACCACGCTTTTCAACCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCg                                                                                                                                                 >  5:350363/1‑149 (MQ=255)
  cTCAAGCAGTTGTTCTGCCGTGGCACCGCGCGGCCCTTCCAGTGGCGCACAATACCAGTCATCAACCACGCTTTTCAACCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCACCCGTTTTTc                                                                                                                                                  <  7:195941/148‑1 (MQ=255)
   tCAAGCAGTTGTTCTGCCGTGGCACCGCGCGGCCCTTCCAGTGGCGCACAATACCAGTCATCAACCACGCTTTTCAACCAGGCCAGAGATCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCgg                                                                                                                                                >  3:263113/1‑149 (MQ=255)
     aaGCAGTTGTTCTGCCGTGGCACCGCGCGGCCCTTCCAGTGGCGCACAATACCAGTCATCAACCACGCTTTTCAACCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCgg                                                                                                                                                <  3:149122/147‑1 (MQ=255)
     aaGCAGTTGTTCTGCCGTGGCACCGCGCGGCCCTTCCAGTGGCGCACAATACCAGTCATCAACCACGCTTTTCAACCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCgg                                                                                                                                                >  4:308987/1‑147 (MQ=255)
     aaGCAGTTGTTCTGCCGTGGCACCGCGCGGCCATTCCAGTGGCGCACCATACCAGACATCAACCACGCTTTTCAACCAGCCCAGAGTTACGGCAATATCTTTATCATGTAGCATAACGATAACCGCAAGCACTCGCCCGTATTTCgg                                                                                                                                                >  8:29721/1‑147 (MQ=255)
       gCAGTTGTTCTGCCGTGGCACCGCGCGGCCCTTCCAGTGGCGCACAATACCAGTCATCAACCACGCTTTTCAACCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAgc                                                                                                                                            >  1:151932/1‑149 (MQ=255)
        cAGTTGTTCTGCCGGGGCACCGCGCGGCCCTTCCAGTGGCGCACAATACCAGTCATCAACCACGCTTTTCACCCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAgcg                                                                                                                                           <  6:123678/149‑1 (MQ=255)
              ttCTGCCGTGGCACCGCGCGGCCCTTCCAGTGGCGCACAATACCAGTCATCAACCACGCTTTTCAACCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCa                                                                                                                                     <  5:338772/149‑1 (MQ=255)
              ttCTGCCGTGGCACCGCGCGGCCCTTCCAGTGGCGCACAATACCAGTCATCAACCACGCTTTTCAACCAGGCCAGAGATCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTc                                                                                                                                      <  4:263113/148‑1 (MQ=255)
              ttCTGCCGTGGCACCGCGCGGCCCTTCCAGTGGCGCACAATACAAGTCATCAACCACGCTTTTCAACCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCa                                                                                                                                     >  4:29658/1‑149 (MQ=255)
                 tGCCGTGGCACCGCGCGGCCCTTCCAGTGGCGCACAATACCAGTCATCAACCACGCTTTTCAACCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATAc                                                                                                                                  >  2:342639/1‑149 (MQ=255)
                   ccGTGGCACCGCGCGGCCCTTCCAGTGGCGCACAATACCAGTCATCAACCACGCTTTTCAACCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGc                                                                                                                                >  8:329443/1‑149 (MQ=255)
                      tGGCACCGCGCGGCCCTTCCAGTGGCGCACAATACCAGTCATCAACCACGCTTTTCAACCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAgcg                                                                                                                                           >  1:10239/1‑135 (MQ=255)
                      tGGCACCGCGCGGCCCTTCCAGTGGCGCACAATACCAGTCATCAACCACGCTTTTCAACCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAgcg                                                                                                                                           <  2:10239/135‑1 (MQ=255)
                          aCCGCGCGGCCCTTCCAGTGGCGCACAATACCAGTCATCAACCACGCTTTTCAACCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTg                                                                                                                          >  2:256181/1‑148 (MQ=255)
                           ccGCGCGGCCCTTCCAGTGGCGCACAATACCAGTCATCAACCACGCTTTTCAACCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGt                                                                                                                           >  6:109680/1‑146 (MQ=255)
                           ccGCGCGGCCCTTCCAGTGGCGCACAATACCAGTCATCAACCACGCTTTTCAACCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTgag                                                                                                                        <  4:318597/149‑1 (MQ=255)
                           ccGCGCGGCCCTTCCAGTGGCGCACAATACCAGTCATCAACCACGCTTTTCAACCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTgag                                                                                                                        <  3:308987/149‑1 (MQ=255)
                            cgcgcgGCCCTTCCAGTGGCGCACAATACCAGTCATCAACCACGCTTTTCAACCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTgaga                                                                                                                       <  2:151932/149‑1 (MQ=255)
                               gcgGCCCTTCCAGTGGCGCACAATACCAGTCATCAACCACGCTTTTCAACCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTGAGat                                                                                                                      <  2:154855/147‑1 (MQ=255)
                               gcgGCCCTTCCAGTGGCGCACAATACCAGTCATCAACCACGCTTTTCAACCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTGAGat                                                                                                                      >  1:154855/1‑147 (MQ=255)
                                  gCCCTTCCAGTGGCGCACAATACCAGTCATCAACCACGCTTTTCAACCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTGAGATATTcc                                                                                                                 >  8:254194/1‑149 (MQ=255)
                                        ccAGTGGCGCACAATACCAGTCATCAACCACGCTTTTCAACCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTGAGATATTCCGCCGCa                                                                                                           >  8:52483/1‑149 (MQ=255)
                                         cAGTGGCGCACAATACCAGTCATCAACCACGCTTTTCAACCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTGAGATATTCCGCCGCAt                                                                                                          <  5:46455/149‑1 (MQ=255)
                                           gTGGCGCACAATACCAGTCATCAACCACGCTTTTCAACCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTGAGATATTCCGCCGCAtgt                                                                                                        <  1:256181/149‑1 (MQ=255)
                                                  acaATACCAGTCATCAACCACGCTTTTCAACCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTGAGATATTCCGCCGCATGTGGATTAt                                                                                                 <  7:329443/149‑1 (MQ=255)
                                                      tACCAGTCATCAACCACGCTTTTCAACCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTGAGATATTCCGCCGCATGTGGATTATgcgc                                                                                             <  4:157835/149‑1 (MQ=255)
                                                      tACCAGTCATCAACCACGCTTTTCAACCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTGAGATATTCAGCCGCATGTGGATTATgagc                                                                                             >  1:171364/1‑149 (MQ=255)
                                                      tACCAGTCATCAACCACGCTTTTCAACCAGGCCAGAGTTCCGGCAATAGCTTTATCATGTAGCATACCGATACCCGCCAGCACGCGCCCGGTTTTCGGTAGCGCTTTCATACGCCCGGTGAGATATTCCGCCGCATGTTGATTATgcgc                                                                                             <  7:29721/149‑1 (MQ=255)
                                                        ccAGTCATCAACCACGCTTTTCAACCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTGAGATATTCCGCCGCATGTGTATTATgcgcta                                                                                           >  2:328318/1‑147 (MQ=255)
                                                           gTCATCAACCACGCTTTTCAACCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTGAGATATTCCGCCGCATGTGGATTATGCGCGACAt                                                                                        <  6:350363/149‑1 (MQ=255)
                                                            tcatcaACCACGCTTTTCAACCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCTCTTTCATACGCCCGGTGAGATATTCCGCCGCAtgt                                                                                                        <  1:56547/132‑1 (MQ=255)
                                                            tcatcaACCACGCTTTTCAACCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTGAGATATTCCGCCGCAtgt                                                                                                        >  2:56547/1‑132 (MQ=255)
                                                              atcaACCACGCTTTTCAACCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTGAGATATTCCGCCGCATGTGGATTATGCGCGACa                                                                                         >  8:230895/1‑145 (MQ=255)
                                                              atcaACCACGCTTTTCAACCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTGAGATATTCCGCCGCATGTGGATTATGCGCGACa                                                                                         <  7:230895/145‑1 (MQ=255)
                                                                 aaCCACGCTTTTCAACCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTGAGATATTCCGCCGCATGTGGATTATGCGCGACATCaaaaa                                                                                  >  2:205844/1‑149 (MQ=255)
                                                                   ccACGCTTTTCAACCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTGAGATATTCCGCCGCATGTGGATTATGCGCGACATCAAAAAt                                                                                 >  5:122613/1‑148 (MQ=255)
                                                                    cACGCTTTTCAACCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTGAGATATTCCGCCGCATGTGGATTATGCGCGACATCAAAAATaa                                                                               >  1:178642/1‑149 (MQ=255)
                                                                      cGCTTTTCAACCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTCTCGGTAGCGCTTTCATACGCCCGGTGAGATATTCCGCCGCATGTGGATTATGCGCGACATCAAAAATAAcg                                                                             >  5:112451/1‑149 (MQ=255)
                                                                         ttttCAACCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTGAGATATTCCGCCGCATGTGGATTATGCGCGACATCAAAAATAACGCGt                                                                          >  8:366846/1‑149 (MQ=255)
                                                                           ttCAACCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATACCCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTGAGATATTCCACCGCATGTGGATTATGCCCGACATCAAAAATAACTCGTgg                                                                        >  6:301792/1‑149 (MQ=255)
                                                                           ttCAACCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTGAGATCTTCCGCCGCATGTGGATTATGCGCGACATCAAAAATAACGCGTg                                                                         <  6:127222/148‑1 (MQ=255)
                                                                            tCAACCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTGAGATATTCCGCCGCATGTGGATTATGCGCGACATCAAAACTAACGCGTgg                                                                        <  5:109680/148‑1 (MQ=255)
                                                                              aaCCAGGCCAGAGTTCCGGCAATATCTTTATCATGTCGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTGAGATATTCCGCCGCATGTGGATTATGCGCGACATCAAAAATAACGCGTGGCGa                                                                     <  4:181743/149‑1 (MQ=255)
                                                                              aaCCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTGAGATATTCCGCCGCATGTGGATTATGCGCGACATCAAAAATAACGCGTGGCGa                                                                     >  8:7371/1‑149 (MQ=255)
                                                                              aaCCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTGAGATATTCCGCCGCATGTGGATTATGCGCGACATCAAAAATAACGCGTGGCGa                                                                     >  8:12913/1‑149 (MQ=255)
                                                                                ccAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTGAGATATTCCGCCGCATGTGGATTATGCGCGACATCAAAAATAACGCGTGGCGACt                                                                   >  1:135238/1‑149 (MQ=255)
                                                                                 cAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACAGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTGAGATATTCCGCCGCATGTGGATTATGCGCGACATCAAAAATAACGCGTGGCGACTc                                                                  >  1:238311/1‑149 (MQ=255)
                                                                                     ccAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTGAGATATTCCGCCGCATGTGGATTATGCGCGACATCAAAAATAACGCGTGGCGACTCGCTc                                                              <  7:254194/149‑1 (MQ=255)
                                                                                      cAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTGAGATATTCCGCCGCATGTGGATTATGCGCGACATCAAAAATAACGCGTGGCGACTCGCTca                                                             <  2:135238/149‑1 (MQ=255)
                                                                                       agagTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTGAGATATTCCGCCGCATGTGGATTATGCGCGACATCAAAAATAACGCGTGGCGACTCGCTca                                                             >  2:246195/1‑148 (MQ=255)
                                                                                          gTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTGAGATATTCCGCCGCATGTGGATTATGCGCGACATCAAAAATAACGCGTGGCGACTCGCTCACAAt                                                         <  1:342639/149‑1 (MQ=255)
                                                                                              cGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTGAGATATTCCGCCGCATGTGGATTATGCGCGACATCAAAAATAACGCGTGGCGACTCGCTCACAATCTg                                                      >  5:162240/1‑148 (MQ=255)
                                                                                              cGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTGAGATATTCCGCCGCATGTGGATTATGCGCGACATCAAAAATAACGCGTGGCGACTCGCTCACAATCTg                                                      <  4:299246/148‑1 (MQ=255)
                                                                                              cGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTGAGATATTCCGCCGCATGTGGATTATGCGCGACATCAAAAATAACGCGTGGCGACTCGCTCACAATCTg                                                      <  1:205844/148‑1 (MQ=255)
                                                                                               ggCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTGAGATATTCCGCCGCATGTGGATTATGCGCGCCATCAAAAATAACGCGTGGCGACTCGCTCACAATCTgg                                                     >  2:306272/1‑148 (MQ=255)
                                                                                               ggCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCATTCATACGCCCGGTGAGATATTCCGCCGCATGTGGATTATGCGCGACATCAAAAATAACGCGTGGCGACTCGCTCACAATCTgg                                                     >  2:160993/1‑148 (MQ=255)
                                                                                               ggCAAGATCTTTATCATGTAGCATCCCGATATCCGGCAGCACGCGCGCGTTGTTCGGTAGCGCTTTCATACGCGCGGTGAGATCTTCCGCCGCAGGTGGCTTATGCGCGACCTCAACACTAGCGCGTGGCGACTCGCGCACAATCTg                                                      >  4:260036/1‑147 (MQ=255)
                                                                                                 cAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTGTTTCGGTAGCGCTTTCATACGCCCGGTGAGATATTCCGCCGCCTGTGGATTATGCGCGACATCAAAAATAACGCGTGGCGACTCGCTCACAATCTGGaaa                                                  <  1:334560/149‑1 (MQ=255)
                                                                                                  aaTATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCTGTGAGATATTCCGCCGCATGTGGATTATGCGCGACATCAAAAATAACGCGTGGCGACTCGCTCACAATCTGGAAAc                                                 <  7:366846/149‑1 (MQ=255)
                                                                                                       cTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTGAGATATTCCGCCGCACGTGGATTATGCGCGACATCAAAAATAACGCGTGGCGAc                                                                    <  7:101826/125‑1 (MQ=255)
                                                                                                        tttATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTGAGATATTCCGCCGCATGTGGATTATGCGCGACATCAAAAATAACGCGTGGCGACTCGCTCACAATCTGGAAACGTCCCg                                           >  5:127955/1‑149 (MQ=255)
                                                                                                          aaTCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTGAGATATTCCGCCGCATGTGGATTATGCGCGACATCAAAAATAACGCGTGGCGACTCGCTCACAATCTGGAAACGTccc                                            <  5:81362/145‑1 (MQ=255)
                                                                                                              aTGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTGAGATATTCCGCCGCATGTGGATTATGCGCGACATCAAAAATAACGCGTGGCGACTCGCTCACAATCTGGAAACGTCCCGGCaaaa                                     <  3:29658/149‑1 (MQ=255)
                                                                                                                      tACCGATAACCGCCAGCACGCGCCCGTTTCTCGGTAGCGCTTTCATACGCCCGGTGAGATATTCCGCCGCATGTGGATTATGCGCGACATCAAAAATAACGCGTGGCGACTCGCTCACAATCTGGAAACGTCCCGGCAAAATTGCGCTg                             <  6:112451/149‑1 (MQ=255)
                                                                                                                       aCCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTGAGATATTCCGCCGCATGTGGATTATGCGCGACATCAAAAATAACGCGTGGCGACTCGCTCACAATCTGGAAACGTCCCGGCAAAATTGCGCTgg                            >  1:278194/1‑149 (MQ=255)
                                                                                                                         cGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTGAGATATTCCGCCGCATGTGGATTATGCGCGACATCAAAAATAACGCGTGGCGACTCGCTCACAATCTGGAAACGTCCCGGCAAAATTGCGCTGGc                           >  2:288265/1‑148 (MQ=255)
                                                                                                                              aCCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTGAGATATTCCGCCGCATGTGGATTATGCGCGACATCAAAAATAACGCGTGGCGACTCGCTCACAATCTGGAAACGTCCCGGCAAAATTGCGCTGGCAATccc                     <  1:328318/149‑1 (MQ=255)
                                                                                                                                cGCCAGCACGCGCCCGTTTCTCGGTAGCGCTTTCATACGCCCGGTGAGATATTCCGCCGCATGTGGATTATGCGCGACATCAAAAATAACGCGTGGCGACTCGCTCACAATCTGGAAACGTCCCGGCAAAATTGCGCTGGCAATc                       <  6:308138/145‑1 (MQ=255)
                                                                                                                                cGCCAGCACGCGCCCGTTTCTCGGTAGCGCTTTCATACGCCCGGTGAGATATTCCGCCGCATGTGGATTATGCGCGACATCAAAAATAACGCGTGGCGACTCGCTCACAATCTGGAAACGTCCCGGCAAAATTGCGCTGGCAATc                       >  5:308138/1‑145 (MQ=255)
                                                                                                                                           gcCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTGAGATATTCCGCCGCATGTGGATTATGCGCGACATCAAAAATAACGCGTGGCGACTCGCTCACAATCTGGAAACGTCCCGGCAAAATTGCGCTGGCAATCCCGTCGCGAATGGCa        >  8:200112/1‑149 (MQ=255)
                                                                                                                                             ccGTTTTTCGGTAGCGCTTTCATACGCCCGGTGAGATATTCCGCCGCATGTGGATTATGCGAGACATCAAAAATAACGCGTGGCGACTCGCTCACAATCTGGAAACGTCCCGGCAAAATTGCGCTGGCAATCCCGTCGCGAATGGCAtt      >  6:38855/1‑149 (MQ=255)
                                                                                                                                                  tttCGGTAGCGCTTTCATACGCCCGGTGAGATATTCCGCCGCATGTGGATTATGCGCGACATCAAAAATAACGCGTGGCGACTCGCTCACAATCTGGAAACGTCCCGGCAAAATTGCGCTGGCAATCCCGTCGCGAATGGCATTTTCa  <  1:246195/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TGCTCAAGCAGTTGTTCTGCCGTGGCACCGCGCGGCCCTTCCAGTGGCGCACAATACCAGTCATCAACCACGCTTTTCAACCAGGCCAGAGTTCCGGCAATATCTTTATCATGTAGCATACCGATAACCGCCAGCACGCGCCCGTTTTTCGGTAGCGCTTTCATACGCCCGGTGAGATATTCCGCCGCATGTGGATTATGCGCGACATCAAAAATAACGCGTGGCGACTCGCTCACAATCTGGAAACGTCCCGGCAAAATTGCGCTGGCAATCCCGTCGCGAATGGCATTTTCA  >  NC_000913/2431831‑2432124

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