Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,541,559 G→A *201* (TAG→TAA fimB → tyrosine recombinase/inversion of on/off regulator of fimA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,541,5590GA98.9% 291.9 / ‑4.5 92*201* (TAG→TAAfimBtyrosine recombinase/inversion of on/off regulator of fimA
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (41/50);  total (42/50)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

TAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGT  >  NC_000913/4541423‑4541687
                                                                                                                                        |                                                                                                                                
tAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATaaag                                                                                                                                >  1:1250824/1‑139 (MQ=255)
 aGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAaaga                                                                                                                               >  2:976223/1‑139 (MQ=255)
  gATACGCGACTTATCCAGGATTATCTTGGTCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGAt                                                                                                                              >  2:982326/1‑139 (MQ=255)
  gATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGAt                                                                                                                              <  2:1448086/139‑1 (MQ=255)
      cgcgACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGtt                                                                                                                          >  2:40279/1‑139 (MQ=255)
       gcgACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTg                                                                                                                         <  2:627775/139‑1 (MQ=255)
        cgACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTga                                                                                                                        <  2:795504/139‑1 (MQ=255)
        cgACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTga                                                                                                                        <  1:220559/139‑1 (MQ=255)
           cTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGAt                                                                                                                     <  1:176184/139‑1 (MQ=255)
                   ggATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGa                                                                                                             <  1:523164/139‑1 (MQ=255)
                   ggATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGa                                                                                                             <  1:985087/139‑1 (MQ=255)
                     aTTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATg                                                                                                           <  2:1476883/139‑1 (MQ=255)
                        aTCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAAACACCAGGGATGgtg                                                                                                        <  1:982326/139‑1 (MQ=255)
                         tCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTgaga                                                                                                                      <  2:365776/124‑1 (MQ=255)
                         tCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTgaga                                                                                                                      >  1:365776/1‑124 (MQ=255)
                          cTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGtt                                                                                                      <  1:234179/139‑1 (MQ=255)
                              ggCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTaaaa                                                                                                  >  2:156284/1‑139 (MQ=255)
                               gCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTaaaaa                                                                                                 <  1:401189/139‑1 (MQ=255)
                                   cGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTAc                                                                                             >  2:1174420/1‑139 (MQ=255)
                                         atTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACattatt                                                                                       >  2:656217/1‑139 (MQ=255)
                                           tCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCt                                                                                     <  1:976223/139‑1 (MQ=255)
                                             gTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTa                                                                                   >  2:422777/1‑139 (MQ=255)
                                                  aCTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGata                                                                              >  1:590225/1‑139 (MQ=255)
                                                    tGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGa                                                                                >  1:1272107/1‑135 (MQ=255)
                                                    tGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGa                                                                                <  2:1272107/135‑1 (MQ=255)
                                                     gTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTCATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGatatat                                                                           >  1:890857/1‑139 (MQ=255)
                                                      tCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGGGGACGACAGCGTCACGCTGTTTTATAAAGAt                                                                                                                              <  1:304711/87‑1 (MQ=255)
                                                      tCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGGGGACGACAGCGTCACGCTGTTTTATAAAGAt                                                                                                                              >  2:304711/1‑87 (MQ=255)
                                                      tCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATc                                                                          >  2:226538/1‑139 (MQ=255)
                                                        tGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGTTATATCaa                                                                        >  2:1475586/1‑139 (MQ=255)
                                                          gTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGGGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCaaaa                                                                      <  2:120264/139‑1 (MQ=255)
                                                           tataCCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCaaaaa                                                                     >  1:788124/1‑139 (MQ=255)
                                                                cGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTgag                                                                                                                       <  1:1474817/84‑1 (MQ=255)
                                                                cGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTgag                                                                                                                       >  2:1474817/1‑84 (MQ=255)
                                                                cGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATg                                                                                                           >  1:1149374/1‑96 (MQ=255)
                                                                cGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATg                                                                                                           <  2:1149374/96‑1 (MQ=255)
                                                                cGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTaaaaa                                                                                                 >  2:455529/1‑106 (MQ=255)
                                                                cGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTaaaaa                                                                                                 <  1:455529/106‑1 (MQ=255)
                                                                cGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACa                                                                                            <  1:594415/111‑1 (MQ=255)
                                                                cGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACa                                                                                            >  2:594415/1‑111 (MQ=255)
                                                                cGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGAtt                                                                >  2:1190240/1‑139 (MQ=255)
                                                                cGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGAtt                                                                >  1:675167/1‑139 (MQ=255)
                                                                cGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGAtt                                                                >  1:680749/1‑139 (MQ=255)
                                                                cGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGAtt                                                                >  1:76606/1‑139 (MQ=255)
                                                                  ccAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATg                                                                                                           >  2:453924/1‑94 (MQ=255)
                                                                  ccAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATg                                                                                                           <  1:453924/94‑1 (MQ=255)
                                                                   cAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAAc                                                             <  2:590225/139‑1 (MQ=255)
                                                                    aGCAATGCAGGGCGATTTTACGGCATCTGGGATAGAGCCAGAGGAAGACAGCCACACGCTGTGTTATAACGATAGTTAAGATGCCAGAGCTGGTGTTAAAAGCTACTTTATTCTTACGATATATCAAAGATGATTAAca                                                            >  1:1409120/1‑138 (MQ=255)
                                                                     gCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCa                                                           >  1:585675/1‑139 (MQ=255)
                                                                              ggCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCaa                                                                        <  1:328615/117‑1 (MQ=255)
                                                                              ggCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCaa                                                                        >  2:328615/1‑117 (MQ=255)
                                                                               gCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAAtt                                                 >  1:1375897/1‑139 (MQ=255)
                                                                                cGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTg                                                >  1:87723/1‑139 (MQ=255)
                                                                                cGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTg                                                >  2:898520/1‑139 (MQ=255)
                                                                                  tttttACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGAt                                              <  1:40279/139‑1 (MQ=255)
                                                                                   ttttACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATg                                             >  1:1128945/1‑139 (MQ=255)
                                                                                    tttACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGt                                            <  1:898520/139‑1 (MQ=255)
                                                                                    tttACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGt                                            <  1:422777/139‑1 (MQ=255)
                                                                                       aCGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTAt                                         <  2:788124/139‑1 (MQ=255)
                                                                                               tGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAgtg                                 <  2:87723/139‑1 (MQ=255)
                                                                                                 ggATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGtgtg                               <  2:1128945/139‑1 (MQ=255)
                                                                                                   aTAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGAt                             <  1:226538/139‑1 (MQ=255)
                                                                                                     agagCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTga                                                                                                                        >  1:1449689/1‑46 (MQ=255)
                                                                                                     agagCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTga                                                                                                                        <  2:1449689/46‑1 (MQ=255)
                                                                                                     agagCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTc                           <  1:1190240/139‑1 (MQ=255)
                                                                                                       agCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTgag                                                                                                                       >  2:498707/1‑45 (MQ=255)
                                                                                                       agCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTgag                                                                                                                       <  1:498707/45‑1 (MQ=255)
                                                                                                       agCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTga                                                                                                                        <  2:255282/44‑1 (MQ=255)
                                                                                                       agCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTga                                                                                                                        >  1:255282/1‑44 (MQ=255)
                                                                                                       agCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATa                                                                                                                    <  1:1312552/48‑1 (MQ=255)
                                                                                                       agCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATa                                                                                                                    >  2:1312552/1‑48 (MQ=255)
                                                                                                       agCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACa                                                                                            >  2:352767/1‑72 (MQ=255)
                                                                                                       agCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACa                                                                                            <  1:352767/72‑1 (MQ=255)
                                                                                                       agCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAAt                                                  <  2:1078806/114‑1 (MQ=255)
                                                                                                       agCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAAt                                                  >  1:1078806/1‑114 (MQ=255)
                                                                                                       agCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTg                                                >  2:618310/1‑116 (MQ=255)
                                                                                                       agCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTg                                                <  1:618310/116‑1 (MQ=255)
                                                                                                       agCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGAt                                              <  2:1399565/118‑1 (MQ=255)
                                                                                                       agCCAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGAt                                              >  1:1399565/1‑118 (MQ=255)
                                                                                                         ccAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGa                                                                                                             >  1:1397534/1‑53 (MQ=255)
                                                                                                         ccAGAGGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGa                                                                                                             <  2:1397534/53‑1 (MQ=255)
                                                                                                           agagGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGAtt                            >  2:1017911/1‑132 (MQ=255)
                                                                                                           agagGACGACAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGAtt                            <  1:1017911/132‑1 (MQ=255)
                                                                                                             aagacgacAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATAAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTAt                   <  1:41605/137‑1 (MQ=255)
                                                                                                                   acAGCGTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCataat             <  1:174295/139‑1 (MQ=255)
                                                                                                                        gTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACa        <  2:1375897/139‑1 (MQ=255)
                                                                                                                        gTCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACa        <  2:675167/139‑1 (MQ=255)
                                                                                                                         tCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACa                                                                                            <  2:1048203/54‑1 (MQ=255)
                                                                                                                         tCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACa                                                                                            >  1:1048203/1‑54 (MQ=255)
                                                                                                                         tCACGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACAt       >  2:736290/1‑139 (MQ=255)
                                                                                                                            cGCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTaa    <  2:680749/139‑1 (MQ=255)
                                                                                                                             gCTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAg   <  1:656217/139‑1 (MQ=255)
                                                                                                                              cTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGt  <  1:1174420/139‑1 (MQ=255)
                                                                                                                              cTGTTTTATAAAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGt  <  1:736290/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                
TAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGT  >  NC_000913/4541423‑4541687

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: