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.

TTATTGGCCATTGCATCCAGTACCAATACAACAACAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTT  >  NC_000913/558521‑558770
                                                                                                                                      |                                                                                                                    
ttATTGGGCATTGCATCCAGTACCAATACAACAACAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCAtt                                                                                                        <  4:381399/148‑1 (MQ=255)
       ccATTGCATCCAGTACCAATACAACAACAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGttctt                                                                                                            <  4:551165/137‑1 (MQ=255)
       ccATTGCATCCAGTACCAATACAACAACAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGttctt                                                                                                            >  3:551165/1‑137 (MQ=255)
         aTTGCATCCAGTACCAATACAACAACAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCAc                                                                                               >  2:553766/1‑148 (MQ=255)
                ccAGTACCAATACAACAACAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATaa                                                                                        >  4:1040215/1‑148 (MQ=255)
                     aCCAATACAACATCAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCa                                                                                   <  3:1040215/148‑1 (MQ=255)
                        aaTACAACAACAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCg                                                                                >  4:488919/1‑148 (MQ=255)
                          tACAACAACAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTta                                                                                           <  3:674508/135‑1 (MQ=255)
                          tACAACAACAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTta                                                                                           >  4:674508/1‑135 (MQ=255)
                           acaacaacaGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATaa                                                                                        <  3:1138070/137‑1 (MQ=255)
                           acaacaacaGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATaa                                                                                        >  4:1138070/1‑137 (MQ=255)
                            caacaacaGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATaa                                                                                        >  4:900023/1‑136 (MQ=255)
                            caacaacaGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATaa                                                                                        <  3:900023/136‑1 (MQ=255)
                               caacaGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAg                                                                         >  2:44152/1‑148 (MQ=255)
                               caacaGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAg                                                                         >  2:481792/1‑148 (MQ=255)
                               caacaGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAg                                                                         <  3:488919/148‑1 (MQ=255)
                                 acaGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCa                                                                                                                    >  3:1021613/1‑103 (MQ=255)
                                 acaGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCa                                                                                                                    <  4:1021613/103‑1 (MQ=255)
                                 acaGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGt                                                                                             >  1:602124/1‑126 (MQ=255)
                                 acaGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGt                                                                                             <  2:602124/126‑1 (MQ=255)
                                 acaGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTg                                                                       >  4:917079/1‑148 (MQ=255)
                                    gCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCa                                                                                   <  3:389293/133‑1 (MQ=255)
                                    gCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCa                                                                                   >  4:389293/1‑133 (MQ=255)
                                    gCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCAt                                                                    <  1:481792/148‑1 (MQ=255)
                                     cAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGAGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATc                                                                   <  1:553766/148‑1 (MQ=255)
                                       aCGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACGTTGGTCCCGGGGATCAAGttctg                                                                                                            >  2:984168/1‑104 (MQ=38)
                                       aCGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACGTTGGTCCCGGGGATCAAGttctg                                                                                                            <  1:902305/105‑2 (MQ=38)
                                       aCGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACGTTGGTCCCGGGGATCAAGttctg                                                                                                            <  1:984168/105‑2 (MQ=37)
                                       aCGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACGTTGGTCCCGGGGATCAAGttctg                                                                                                            >  2:902305/1‑104 (MQ=38)
                                          ggTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATcagcag                                                              <  1:44152/148‑1 (MQ=255)
                                                   aTTGAAATACTTGATAATACATCCGCAATTCTCAAAACTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGgcaagcaa                                                     <  3:917079/148‑1 (MQ=255)
                                                       aaaTACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAATCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGgca                                                          >  1:740655/1‑139 (MQ=255)
                                                       aaaTACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGgca                                                          <  2:740655/139‑1 (MQ=255)
                                                        aaTACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAg                                                                         >  3:236819/1‑123 (MQ=255)
                                                        aaTACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAg                                                                         <  4:236819/123‑1 (MQ=255)
                                                         aTACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGt                                                                                             <  3:234364/102‑1 (MQ=255)
                                                         aTACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGt                                                                                             >  4:234364/1‑102 (MQ=255)
                                                            cTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCAt                                                                    <  4:83450/124‑1 (MQ=255)
                                                            cTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCAt                                                                    >  3:83450/1‑124 (MQ=255)
                                                                 taataCATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACt                                       >  1:472874/1‑148 (MQ=255)
                                                                   ataCATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACttt                                     <  2:472874/148‑1 (MQ=255)
                                                                         ccGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTtattat                               >  2:727104/1‑148 (MQ=255)
                                                                           gCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCAttttt                                                                                                     <  4:185969/76‑1 (MQ=255)
                                                                           gCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCAttttt                                                                                                     >  3:185969/1‑76 (MQ=255)
                                                                            cAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGAcgcg                                          >  2:684768/1‑134 (MQ=255)
                                                                            cAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGAcgcg                                          <  1:684768/134‑1 (MQ=255)
                                                                                tctcAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACa                                                                           >  2:820791/1‑97 (MQ=255)
                                                                                tctcAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACa                                                                           <  1:820791/97‑1 (MQ=255)
                                                                                   cAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGATCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGaa                     >  2:860294/1‑148 (MQ=255)
                                                                                     aaCCTGATGGGAATAGCTTCTCAACCAACCAGACCTTGATCCCCGGGACCAACGttcttc                                                                                                           <  3:760930/60‑1 (MQ=255)
                                                                                     aaCCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGttcttc                                                                                                           >  4:760930/1‑60 (MQ=255)
                                                                                     aaCCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAg                                                                         <  3:843850/94‑1 (MQ=255)
                                                                                     aaCCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAg                                                                         >  4:843850/1‑94 (MQ=255)
                                                                                         tGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAg                                                                                      <  2:895208/77‑1 (MQ=255)
                                                                                         tGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAg                                                                                      >  1:895208/1‑77 (MQ=255)
                                                                                          gATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACa                                                                           >  2:648276/1‑87 (MQ=255)
                                                                                          gATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACa                                                                           <  1:648276/87‑1 (MQ=255)
                                                                                           aTGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCGGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCaaa             <  1:727104/148‑1 (MQ=255)
                                                                                              ggAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGt                                                                        <  2:165161/86‑1 (MQ=255)
                                                                                              ggAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGt                                                                        >  1:165161/1‑86 (MQ=255)
                                                                                              ggAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAAt                                                   <  2:462463/107‑1 (MQ=255)
                                                                                              ggAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGATATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAAt                                                   >  1:462463/1‑107 (MQ=255)
                                                                                                aaTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGt                                                                                             >  2:110610/1‑63 (MQ=255)
                                                                                                aaTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGt                                                                                             <  1:110610/63‑1 (MQ=255)
                                                                                                aaTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTtta                                   >  2:1017830/1‑121 (MQ=255)
                                                                                                aaTAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTtta                                   <  1:1017830/121‑1 (MQ=255)
                                                                                                  tAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCa                                                                     <  1:1096006/85‑1 (MQ=255)
                                                                                                  tAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCa                                                                     >  2:1096006/1‑85 (MQ=255)
                                                                                                  tAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTtta                                   >  4:915312/1‑119 (MQ=255)
                                                                                                  tAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTtta                                   <  3:915312/119‑1 (MQ=255)
                                                                                                   aGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGc                                                                                    <  3:535504/69‑1 (MQ=255)
                                                                                                   aGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGc                                                                                    >  4:535504/1‑69 (MQ=255)
                                                                                                   aGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGc                                                                                    <  1:871370/69‑1 (MQ=255)
                                                                                                   aGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGc                                                                                    >  2:871370/1‑69 (MQ=255)
                                                                                                     cTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAgg                                                             >  1:514594/1‑90 (MQ=255)
                                                                                                     cTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAgg                                                             <  2:514594/90‑1 (MQ=255)
                                                                                                      ttCTCAACCAACCAGAACTTGATCCCCGGGATCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACgttgtt  <  1:860294/148‑1 (MQ=255)
                                                                                                        ctcAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACgttgtt  >  1:238436/1‑146 (MQ=255)
                                                                                                        ctcAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACgttgtt  <  2:238436/146‑1 (MQ=255)
                                                                                                         tcAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAg                                                                         >  4:1053669/1‑74 (MQ=255)
                                                                                                         tcAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAg                                                                         <  3:1053669/74‑1 (MQ=255)
                                                                                                          caaccaacCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAgg                                                                                     >  1:521212/1‑61 (MQ=255)
                                                                                                          caaccaacCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAgg                                                                                     <  2:521212/61‑1 (MQ=255)
                                                                                                            accaaccaGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATgag                            >  3:948595/1‑116 (MQ=255)
                                                                                                              caaccaGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGc                                                                      >  4:306512/1‑72 (MQ=255)
                                                                                                              caaccaGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGc                                                                      <  3:306512/72‑1 (MQ=255)
                                                                                                                  caGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGc                                                                      >  4:254772/1‑68 (MQ=255)
                                                                                                                  caGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGc                                                                      <  3:254772/68‑1 (MQ=255)
                                                                                                                   aGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTttatta                                <  2:611552/105‑1 (MQ=255)
                                                                                                                   aGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTttatta                                >  1:611552/1‑105 (MQ=255)
                                                                                                                    gAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCAt                                                                    >  4:842549/1‑68 (MQ=255)
                                                                                                                    gAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCAt                                                                    <  3:842549/68‑1 (MQ=255)
                                                                                                                     aaCTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCAt                                                                    <  3:709975/67‑1 (MQ=255)
                                                                                                                     aaCTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGAACCGGTACAAGTGCAt                                                                    >  4:709975/1‑67 (MQ=255)
                                                                                                                      aCTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGgcaa                                                         <  3:738435/77‑1 (MQ=255)
                                                                                                                      aCTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGgcaa                                                         >  4:738435/1‑77 (MQ=255)
                                                                                                                      aCTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCa         >  4:25090/1‑125 (MQ=255)
                                                                                                                      aCTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCa         <  3:25090/125‑1 (MQ=255)
                                                                                                                         ttGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGa                      <  4:610446/110‑1 (MQ=255)
                                                                                                                         ttGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGa                      >  3:610446/1‑110 (MQ=255)
                                                                                                                                      |                                                                                                                    
TTATTGGCCATTGCATCCAGTACCAATACAACAACAGCAACGGGTGTCGGTATTGAAATACTTGATAATACATCCGCAATTCTCAAACCTGATGGGAATAGCTTCTCAACCAACCAGAACTTGATCCCCGGGACCAACGTTCTTCATTTTTCTGCACGTTATAAAGGCACCGGTACAAGTGCATCAGCAGGGCAAGCAAATGCTGACGCGACTTTTATTATGAGATATGAATAATCAAAACCACGTTGTT  >  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: