Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 558,654 C→T 5.9% T148I (ACC→ATC)  sfmA → FimA homolog, function unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913558,6540CT 5.9% 339.0 / 8.5 101T148I (ACC→ATC) sfmAFimA homolog, function unknown
Reads supporting (aligned to +/- strand):  ref base C (47/48);  new base T (3/3);  total (50/51)
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.56e-01

Only 100 of 101 total aligned reads displayed.

TTATTGGCCATTGCATCCAGTACCAATACAACAACAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTT  >  NC_000913/558521‑558770
                                                                                                                                      |                                                                                                                    
ttATTGGGCATTGCATCCAGTACCAATACAACAACAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCAtt                                                                                                        <  4:381399/148‑1 (MQ=255)
       ccATTGCATCCAGTACCAATACAACAACAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGttctt                                                                                                            <  4:551165/137‑1 (MQ=255)
       ccATTGCATCCAGTACCAATACAACAACAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGttctt                                                                                                            >  3:551165/1‑137 (MQ=255)
         aTTGCATCCAGTACCAATACAACAACAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCAc                                                                                               >  2:553766/1‑148 (MQ=255)
                ccAGTACCAATACAACAACAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATaa                                                                                        >  4:1040215/1‑148 (MQ=255)
                     aCCAATACAACATCAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCa                                                                                   <  3:1040215/148‑1 (MQ=255)
                        aaTACAACAACAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCg                                                                                >  4:488919/1‑148 (MQ=255)
                          tACAACAACAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTta                                                                                           <  3:674508/135‑1 (MQ=255)
                          tACAACAACAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTta                                                                                           >  4:674508/1‑135 (MQ=255)
                           acaacaacaGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATaa                                                                                        <  3:1138070/137‑1 (MQ=255)
                           acaacaacaGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATaa                                                                                        >  4:1138070/1‑137 (MQ=255)
                            caacaacaGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATaa                                                                                        >  4:900023/1‑136 (MQ=255)
                            caacaacaGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATaa                                                                                        <  3:900023/136‑1 (MQ=255)
                               caacaGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAg                                                                         >  2:44152/1‑148 (MQ=255)
                               caacaGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAg                                                                         >  2:481792/1‑148 (MQ=255)
                               caacaGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAg                                                                         <  3:488919/148‑1 (MQ=255)
                                 acaGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCa                                                                                                                    >  3:1021613/1‑103 (MQ=255)
                                 acaGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCa                                                                                                                    <  4:1021613/103‑1 (MQ=255)
                                 acaGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGt                                                                                             >  1:602124/1‑126 (MQ=255)
                                 acaGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGt                                                                                             <  2:602124/126‑1 (MQ=255)
                                 acaGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTg                                                                       >  4:917079/1‑148 (MQ=255)
                                    gCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCa                                                                                   <  3:389293/133‑1 (MQ=255)
                                    gCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCa                                                                                   >  4:389293/1‑133 (MQ=255)
                                    gCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCAt                                                                    <  1:481792/148‑1 (MQ=255)
                                     cAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGAGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATc                                                                   <  1:553766/148‑1 (MQ=255)
                                       aCGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACGTTGGTCCCGGGGATCAA‑Gttctg                                                                                                            >  2:984168/1‑104 (MQ=38)
                                       aCGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACGTTGGTCCCGGGGATCAA‑Gttctg                                                                                                            <  1:902305/105‑2 (MQ=38)
                                       aCGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACGTTGGTCCCGGGGATCAA‑Gttctg                                                                                                            <  1:984168/105‑2 (MQ=37)
                                       aCGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACGTTGGTCCCGGGGATCAA‑Gttctg                                                                                                            >  2:902305/1‑104 (MQ=38)
                                          ggTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATcagcag                                                              <  1:44152/148‑1 (MQ=255)
                                                   aTTGAAATACTTGATAATACATCCGCAATTCTCAAAACTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGgcaagcaa                                                     <  3:917079/148‑1 (MQ=255)
                                                       aaaTACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAATCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGgca                                                          >  1:740655/1‑139 (MQ=255)
                                                       aaaTACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGgca                                                          <  2:740655/139‑1 (MQ=255)
                                                        aaTACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAg                                                                         >  3:236819/1‑123 (MQ=255)
                                                        aaTACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAg                                                                         <  4:236819/123‑1 (MQ=255)
                                                         aTACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGt                                                                                             <  3:234364/102‑1 (MQ=255)
                                                         aTACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGt                                                                                             >  4:234364/1‑102 (MQ=255)
                                                            cTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCAt                                                                    <  4:83450/124‑1 (MQ=255)
                                                            cTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCAt                                                                    >  3:83450/1‑124 (MQ=255)
                                                                 taataCATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACt                                       >  1:472874/1‑148 (MQ=255)
                                                                   ataCATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACttt                                     <  2:472874/148‑1 (MQ=255)
                                                                         ccGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTtattat                               >  2:727104/1‑148 (MQ=255)
                                                                           gCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCAttttt                                                                                                     <  4:185969/76‑1 (MQ=255)
                                                                           gCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCAttttt                                                                                                     >  3:185969/1‑76 (MQ=255)
                                                                            cAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGAcgcg                                          >  2:684768/1‑134 (MQ=255)
                                                                            cAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGAcgcg                                          <  1:684768/134‑1 (MQ=255)
                                                                                tctcAAACCTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACa                                                                           >  2:820791/1‑97 (MQ=255)
                                                                                tctcAAACCTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACa                                                                           <  1:820791/97‑1 (MQ=255)
                                                                                   cAAACCTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGATCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGaa                     >  2:860294/1‑148 (MQ=255)
                                                                                     aaCCTGATGGGAATAGCTTCTCAACCAACCAGACC‑TTGATCCCCGGGACCAACGttcttc                                                                                                           <  3:760930/60‑1 (MQ=255)
                                                                                     aaCCTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGttcttc                                                                                                           >  4:760930/1‑60 (MQ=255)
                                                                                     aaCCTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAg                                                                         <  3:843850/94‑1 (MQ=255)
                                                                                     aaCCTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAg                                                                         >  4:843850/1‑94 (MQ=255)
                                                                                         tGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAg                                                                                      <  2:895208/77‑1 (MQ=255)
                                                                                         tGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAg                                                                                      >  1:895208/1‑77 (MQ=255)
                                                                                          gATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACa                                                                           >  2:648276/1‑87 (MQ=255)
                                                                                          gATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACa                                                                           <  1:648276/87‑1 (MQ=255)
                                                                                           aTGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCGGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCaaa             <  1:727104/148‑1 (MQ=255)
                                                                                              ggAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGt                                                                        <  2:165161/86‑1 (MQ=255)
                                                                                              ggAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGt                                                                        >  1:165161/1‑86 (MQ=255)
                                                                                              ggAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAAt                                                   <  2:462463/107‑1 (MQ=255)
                                                                                              ggAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGATATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAAt                                                   >  1:462463/1‑107 (MQ=255)
                                                                                                aaTAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGt                                                                                             >  2:110610/1‑63 (MQ=255)
                                                                                                aaTAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGt                                                                                             <  1:110610/63‑1 (MQ=255)
                                                                                                aaTAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTtta                                   >  2:1017830/1‑121 (MQ=255)
                                                                                                aaTAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTtta                                   <  1:1017830/121‑1 (MQ=255)
                                                                                                  tAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCa                                                                     <  1:1096006/85‑1 (MQ=255)
                                                                                                  tAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCa                                                                     >  2:1096006/1‑85 (MQ=255)
                                                                                                  tAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTtta                                   >  4:915312/1‑119 (MQ=255)
                                                                                                  tAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTtta                                   <  3:915312/119‑1 (MQ=255)
                                                                                                   aGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGc                                                                                    <  3:535504/69‑1 (MQ=255)
                                                                                                   aGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGc                                                                                    >  4:535504/1‑69 (MQ=255)
                                                                                                   aGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGc                                                                                    <  1:871370/69‑1 (MQ=255)
                                                                                                   aGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGc                                                                                    >  2:871370/1‑69 (MQ=255)
                                                                                                     cTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAgg                                                             >  1:514594/1‑90 (MQ=255)
                                                                                                     cTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAgg                                                             <  2:514594/90‑1 (MQ=255)
                                                                                                      ttCTCAACCAACCAGAAC‑TTGATCCCCGGGATCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACgttgtt  <  1:860294/148‑1 (MQ=255)
                                                                                                        ctcAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACgttgtt  >  1:238436/1‑146 (MQ=255)
                                                                                                        ctcAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACgttgtt  <  2:238436/146‑1 (MQ=255)
                                                                                                         tcAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAg                                                                         >  4:1053669/1‑74 (MQ=255)
                                                                                                         tcAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAg                                                                         <  3:1053669/74‑1 (MQ=255)
                                                                                                          caaccaacCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAgg                                                                                     >  1:521212/1‑61 (MQ=255)
                                                                                                          caaccaacCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAgg                                                                                     <  2:521212/61‑1 (MQ=255)
                                                                                                            accaaccaGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATgag                            >  3:948595/1‑116 (MQ=255)
                                                                                                              caaccaGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGc                                                                      >  4:306512/1‑72 (MQ=255)
                                                                                                              caaccaGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGc                                                                      <  3:306512/72‑1 (MQ=255)
                                                                                                                  caGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGc                                                                      >  4:254772/1‑68 (MQ=255)
                                                                                                                  caGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGc                                                                      <  3:254772/68‑1 (MQ=255)
                                                                                                                   aGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTttatta                                <  2:611552/105‑1 (MQ=255)
                                                                                                                   aGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTttatta                                >  1:611552/1‑105 (MQ=255)
                                                                                                                    gAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCAt                                                                    >  4:842549/1‑68 (MQ=255)
                                                                                                                    gAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCAt                                                                    <  3:842549/68‑1 (MQ=255)
                                                                                                                     aaC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCAt                                                                    <  3:709975/67‑1 (MQ=255)
                                                                                                                     aaC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGAACCGGTACAAGTGCAt                                                                    >  4:709975/1‑67 (MQ=255)
                                                                                                                      aC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGgcaa                                                         <  3:738435/77‑1 (MQ=255)
                                                                                                                      aC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGgcaa                                                         >  4:738435/1‑77 (MQ=255)
                                                                                                                      aC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCa         >  4:25090/1‑125 (MQ=255)
                                                                                                                      aC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCa         <  3:25090/125‑1 (MQ=255)
                                                                                                                         ttGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGa                      <  4:610446/110‑1 (MQ=255)
                                                                                                                         ttGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGa                      >  3:610446/1‑110 (MQ=255)
                                                                                                                                      |                                                                                                                    
TTATTGGCCATTGCATCCAGTACCAATACAACAACAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAAC‑TTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTT  >  NC_000913/558521‑558770

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