Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 73,013 G→C 9.5% R503R (CGC→CGG)  thiP ← thiamine/thiamine pyrophosphate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091373,0130GC 9.5% 174.6 / 6.6 67R503R (CGC→CGG) thiPthiamine/thiamine pyrophosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (29/31);  new base C (2/5);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGC  >  NC_000913/72866‑73160
                                                                                                                                                   |                                                                                                                                                   
gTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGc                                                                                                                                                    <  1:140791/149‑1 (MQ=255)
  tAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGa                                                                                                                                                  <  2:1089121/149‑1 (MQ=255)
  tAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGa                                                                                                                                                  <  2:1172651/149‑1 (MQ=255)
    aGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATa                                                                                                                                                >  2:463353/1‑149 (MQ=255)
      cTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAgg                                                                                                                                              >  2:590237/1‑149 (MQ=255)
       tAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGa                                                                                                                                             >  1:1160471/1‑149 (MQ=255)
       tAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTATTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCTTCCTGGCTGCGATAGGa                                                                                                                                             <  2:734365/149‑1 (MQ=255)
       tAAAACGCATCGGCAAATGGTGGTAAAGCCAGGAGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACTGTAAACAGCAGAAAACAGCGCCGCAGCAGAATTAACGCGGTGACCGCACCGTCATGCAAGCGATacgc                                                                                                                                             >  1:801509/1‑146 (MQ=255)
       aaataCGCAGCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGa                                                                                                                                             <  2:869220/145‑1 (MQ=255)
         aaaCGCATCGGAAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGc                                                                                                                                           <  2:752945/149‑1 (MQ=255)
           aCGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTCCTATCGCAGCCa                                                                                                                                         <  1:856025/149‑1 (MQ=255)
           aCGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTCCTATCGCAGCCa                                                                                                                                         <  1:590872/149‑1 (MQ=255)
           aCGAATCGTCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTCCTATCGCAGCCa                                                                                                                                         <  2:117199/149‑1 (MQ=255)
            cGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTCCTATCGCAGCCag                                                                                                                                        <  1:559791/149‑2 (MQ=255)
             gCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTACCGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAAt                                                                                                                                       >  1:508040/1‑149 (MQ=255)
             gCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGAGATAGGAGCCAAt                                                                                                                                       >  1:503635/1‑149 (MQ=255)
                  ggCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCt                                                                                                                                  <  2:1160471/149‑1 (MQ=255)
                  ggCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTGTTAACATATCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCt                                                                                                                                  <  1:35562/149‑1 (MQ=255)
                   gCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTg                                                                                                                                 <  2:845756/149‑1 (MQ=255)
                   gCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTCCGATAGGAGCCAATTTGCTg                                                                                                                                 >  2:161670/1‑149 (MQ=255)
                     aaaTGGTGGTAAAGCCAGGTGATATCAGTCAGATTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGGAACCAGCAGACAACAGAGCAGCAGCAGAATTAACGCCGTGCCCGCACCGTCCTGGCACCGATAGGAGCCAAATTGCTGGt                                                                                                                               >  2:947705/1‑149 (MQ=255)
                      aaTGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTa                                                                                                                              <  2:270661/149‑1 (MQ=255)
                         ggtggtAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTagag                                                                                                                           >  1:323099/1‑149 (MQ=255)
                           tggtAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGAt                                                                                                                         >  2:111788/1‑149 (MQ=255)
                             gtAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATaa                                                                                                                       >  1:747299/1‑149 (MQ=255)
                             gtAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTCACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATaa                                                                                                                       >  2:387398/1‑149 (MQ=255)
                                aaGCCAGGTGATATCAGTCAGTTTTAAAATTTCGACCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAAc                                                                                                                    >  1:655326/1‑149 (MQ=255)
                                        tGATATCAGTCAGTTATAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATAAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGt                                                                                                            <  1:103629/149‑1 (MQ=255)
                                         gATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTg                                                                                                           >  2:988401/1‑149 (MQ=255)
                                          atatCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAACCAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGGGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGc                                                                                                          <  1:278088/149‑1 (MQ=255)
                                                   agtTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAACCGGCAGGGTGCGGAAatcat                                                                                                 <  2:963855/149‑1 (MQ=255)
                                                   agtTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAatcat                                                                                                 <  1:144447/149‑1 (MQ=255)
                                                   agtTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAatcat                                                                                                 <  1:321281/149‑1 (MQ=255)
                                                    gtTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAAAAGAGCAGCAGCACAATTAACGCGGTGACCGCACCGTCCTGGCTGCTATAGGAGCCAATTTGCTGGTAGAGAAAAAACGGCAGGGTGCGGATatcatc                                                                                                >  1:1062995/1‑149 (MQ=255)
                                                     ttttAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGCATTAACGCGGTGCCCGCCCCGTCCTGGCTGGGATAGGAGCCATTTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCg                                                                                               >  2:95932/1‑149 (MQ=255)
                                                         aaCATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGCGCAGCAGCTGAATTAACGCGGTGACCGCACCGTCCAGGCTGCGATAGGAGCCAATTTGCGGGGAGAGAGACAACGGCAGGCGGCGGAAATCCTCGtatc                                                                                           >  2:209053/1‑146 (MQ=255)
                                                           cATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCg                                                                                         <  1:955288/149‑1 (MQ=255)
                                                            aTTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGa                                                                                        <  1:812257/149‑1 (MQ=255)
                                                                 gCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGaacaac                                                                                   <  2:617766/149‑1 (MQ=255)
                                                                  ccccGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCTCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCAACGTTACCGAACACCg                                                                                  >  2:1080836/1‑149 (MQ=255)
                                                                   cccGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGc                                                                                 <  2:535766/149‑1 (MQ=255)
                                                                     cGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGcca                                                                               >  2:635110/1‑149 (MQ=255)
                                                                           tttttCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGc                                                                         >  1:169676/1‑149 (MQ=255)
                                                                            ttttCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGTTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGcc                                                                        >  1:461363/1‑149 (MQ=255)
                                                                             tttCAATCACGGGAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGGGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCa                                                                       <  1:607857/149‑1 (MQ=255)
                                                                                   tCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCa                                                                 <  2:921965/149‑1 (MQ=255)
                                                                                   tCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCa                                                                 <  2:627466/149‑1 (MQ=255)
                                                                                        gTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGGGACCGCACCGTCCTGGCTGCGATAGGAGCAAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAAt                                                            <  2:494203/149‑1 (MQ=255)
                                                                                             cagcagAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGCACAACGCCACCACGCCAAAATCTCCAATCGTCa                                                       >  2:22854/1‑149 (MQ=255)
                                                                                              agcagAAAACAGAGCAGTAGCAGAATTAACTCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAg                                                      <  1:801442/149‑1 (MQ=255)
                                                                                              agcagAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAg                                                      >  2:502846/1‑149 (MQ=255)
                                                                                                cagAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCCCCAATCGACAGcc                                                    >  1:575198/1‑148 (MQ=255)
                                                                                                   aaaaCAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTACGGAAATCAACGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGc                                                 >  1:1106817/1‑149 (MQ=255)
                                                                                                    aaaCAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCTGAAATCATCGTTACCGAACAACGCAACCACGCCAAAATCACCACTCGACAGAACCCa                                                >  1:671245/1‑149 (MQ=255)
                                                                                                        agagCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCa                                            <  1:1137934/149‑1 (MQ=255)
                                                                                                             agcagcagAAGAAACGCGGTGACCGCACCGTCCTGGCGGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGc                                       <  2:165228/149‑1 (MQ=255)
                                                                                                                  cagAATTAACGCGGTGACCGCACCGTCCTGCCTGCGATAGGAGCCAATTTGCTGGTATAGATAAAAAGGCAGGGTGCGGAAATCCTCGTTACCGAAAAACGCATCCACGCCAAAATCACGAATCGACAGCCCGCATGAAAAAGACAgcc                                  >  2:616005/1‑148 (MQ=255)
                                                                                                                   agAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACACCGCCACCACGCCAAAATCACCACTCGACAGCACGCATGCAAAGGCCAgcgc                                 >  1:1131465/1‑149 (MQ=255)
                                                                                                                   agAATTAACGCGGTGACCGCACCGTCATGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAgcgc                                 <  1:949944/149‑1 (MQ=255)
                                                                                                                    taaatAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGGCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGcc                                <  1:488954/145‑1 (MQ=255)
                                                                                                                     aaTTAACGCGGTGACCGCACCGTCCTGGCGGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACTCCACGACGCCAAAATCACCACTCGACAGCACGCAGGCAAACGCCAGCGccg                               >  1:780277/1‑148 (MQ=255)
                                                                                                                       ttatctcgGTGAACTCACCGTCCTGGCTGCGATAGTATCCAATTTGCTGGTTGAGATAAAACGGCAGGGAGCGGAAATCATCGTTACCGAACAACGCCACCACGCCTAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTgc                             <  1:865269/143‑1 (MQ=255)
                                                                                                                                tGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAATCGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGa                    <  2:905811/149‑1 (MQ=255)
                                                                                                                                ggACCGCCCCGTCGTTGGTCCGATACGAGCCAATTTGCTGGTAGTTATAAATCGGCAGGGTGCGGAAATCATAGTTAGCGAACTACGCCCCCACGCCAAAATCACAAATCGACAGCACGCATGCTAAGGCCAGCGCCTGCGCCAGTGGa                    <  2:801509/148‑1 (MQ=255)
                                                                                                                                      cACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTc              <  1:463353/149‑1 (MQ=255)
                                                                                                                                      cACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTc              >  2:1031184/1‑149 (MQ=255)
                                                                                                                                                  tGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGc  >  1:1151357/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGC  >  NC_000913/72866‑73160

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