Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,584,031 A→G 31.5% intergenic (‑72/+176) safA ← / ← ydeP two‑component system connector membrane protein, EvgSA to PhoQP/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,584,0310AG31.5% 183.3 / 118.9 148intergenic (‑72/+176)safA/ydePtwo‑component system connector membrane protein, EvgSA to PhoQP/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (45/55);  new base G (25/21);  total (71/77)
Fisher's exact test for biased strand distribution p-value = 3.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 153 total aligned reads displayed.

TTTTGTTTTTCACTGTGGTCGCATGCATTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATGTAGTTTAATTTAATCCATGT  >  NC_000913/1583932‑1584130
                                                                                                   |                                                                                                   
ttttGTTTTTCACTGTGGTCGCATGCATTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGaa                                                                                                     >  1:641469/1‑100 (MQ=255)
ttttGTTTTTCACTGTGGTCGCATGCATTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGaa                                                                                                     >  1:3725776/1‑100 (MQ=255)
 tttGTTTTTCACTGTGGTCGCATGCATTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGaaa                                                                                                    <  1:784751/100‑1 (MQ=255)
   tGTTTTTCACTGTGGTCGCATGCATTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAag                                                                                                  <  1:3007237/100‑1 (MQ=255)
   tGTTTTTCACTGTGGTCGCATGCATTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAag                                                                                                  >  2:347950/1‑100 (MQ=255)
      ttttCACTGTGGTCGCATGCATTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAgg                                                                                               <  2:1359075/100‑1 (MQ=255)
      ttttCACTGTGGTCGCATGCATTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGATTAATATAACAGGCTGGAGCGAAGAgg                                                                                               >  1:4179650/1‑100 (MQ=255)
       tttCACTGTGGTCGCATGCATTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGa                                                                                              <  2:2910464/100‑1 (MQ=255)
        ttCACTGTGGTCGCATGCATTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGaa                                                                                             >  1:2930156/1‑100 (MQ=255)
        ttCACTGTGGTCGCATGCATTTCAAATATGTTTAATTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGaa                                                                                             <  2:20489/100‑1 (MQ=255)
         tCACTGTGGTCGCATGCATTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGaaa                                                                                            >  2:185801/1‑100 (MQ=255)
          cACTGTGGTCGCATGCATTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAAt                                                                                           <  2:2930168/100‑1 (MQ=255)
           aCTGTGGTCGCATGCATTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATa                                                                                          >  2:1230228/1‑100 (MQ=255)
           aCTGTGGTCGCATGCATTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATa                                                                                          <  2:421905/100‑1 (MQ=255)
           aCTGTGGTCGCATGCATTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATa                                                                                          <  2:1452368/100‑1 (MQ=255)
            cTGTGGTCGCATGCATTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAg                                                                                         <  2:2205441/100‑1 (MQ=255)
                 gTCGCATGCATTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTg                                                                                    <  2:3291816/100‑1 (MQ=255)
                  tCGCATGCATTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTgg                                                                                   <  1:4117900/100‑1 (MQ=255)
                  tCGCATACATTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTgg                                                                                   <  1:856831/100‑1 (MQ=255)
                   cGCATGCATTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGt                                                                                  >  1:453907/1‑100 (MQ=255)
                     catgcatTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGc                                                                                >  2:1768265/1‑100 (MQ=255)
                      atgcatTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCt                                                                               >  2:3793172/1‑100 (MQ=255)
                       tgcatTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAACATAACATGCTGGAGAAAAGAGGAAATAGGACTGGTGCtt                                                                              >  1:2998816/1‑100 (MQ=255)
                        gcatTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTc                                                                             <  1:2317698/100‑1 (MQ=255)
                        gcatTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCttn                                                                             <  2:2726427/100‑2 (MQ=255)
                        gcatTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTc                                                                             <  2:2992968/100‑1 (MQ=255)
                         catTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCa                                                                            >  2:2387059/1‑100 (MQ=255)
                          atTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAg                                                                           >  1:2388384/1‑100 (MQ=255)
                          atTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAg                                                                           >  1:1810311/1‑100 (MQ=255)
                            ttCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTa                                                                         <  2:1796257/100‑1 (MQ=255)
                             tCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTAc                                                                        >  2:1326422/1‑100 (MQ=255)
                             tCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTAc                                                                        >  2:2052344/1‑100 (MQ=255)
                                 atatGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAg                                                                    <  2:2803988/100‑1 (MQ=255)
                                 atatGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATACTTAATATAACAGGCTGGAGAAAAGAGGAAATTGGACTGGTGCTTCAGTACTGAg                                                                    <  1:933161/100‑1 (MQ=255)
                                    tGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGTGAAGAGGAAATAGGAATGGTGCTTCAGTACTGCGCgg                                                                 >  1:4012536/1‑100 (MQ=255)
                                     gTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGa                                                                >  1:2109548/1‑100 (MQ=255)
                                     gTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGa                                                                >  2:174218/1‑100 (MQ=255)
                                     gTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGa                                                                >  2:657242/1‑100 (MQ=255)
                                     gTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGa                                                                >  2:3090827/1‑100 (MQ=255)
                                        tatttaGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGtt                                                             <  2:2355665/100‑1 (MQ=255)
                                         atttaGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGttt                                                            <  2:3655943/100‑1 (MQ=255)
                                           ttaGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCt                                                          <  2:1890134/100‑1 (MQ=255)
                                           ttaGCGGATAACGTGAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCt                                                          <  1:2010907/100‑1 (MQ=255)
                                            taGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCtt                                                         >  1:4268527/1‑100 (MQ=255)
                                            taGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCtt                                                         <  1:4123068/100‑1 (MQ=255)
                                            taGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCtt                                                         <  1:347948/100‑1 (MQ=255)
                                             aGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTa                                                        <  2:3542324/100‑1 (MQ=255)
                                              gCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTTCTTCAGTACTGAGCGGAGTTTCTTAc                                                       <  1:469634/100‑1 (MQ=255)
                                              gCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTAc                                                       >  2:2725785/1‑100 (MQ=255)
                                                 gATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGc                                                    <  2:4104258/100‑1 (MQ=255)
                                                 gATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGc                                                    >  2:199623/1‑100 (MQ=255)
                                                 gATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGc                                                    <  1:4249936/100‑1 (MQ=255)
                                                 gATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGc                                                    <  1:3140723/100‑1 (MQ=255)
                                                 gATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGc                                                    <  2:1466192/100‑1 (MQ=255)
                                                  aTAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCt                                                   <  1:1954137/100‑1 (MQ=255)
                                                    aaCGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGt                                                 <  2:2087294/100‑1 (MQ=255)
                                                      cGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAg                                               >  1:3702673/1‑100 (MQ=255)
                                                       gTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAgg                                              >  1:2248070/1‑100 (MQ=255)
                                                       gTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAgg                                              <  2:252271/100‑1 (MQ=255)
                                                        ttAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGc                                             <  2:2272755/100‑1 (MQ=255)
                                                        ttAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGc                                             <  2:1092359/100‑1 (MQ=255)
                                                         tAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGc                                             <  1:733208/99‑1 (MQ=255)
                                                         tAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCa                                            >  1:3332741/1‑100 (MQ=255)
                                                          aaaaaTCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGc                                             <  1:3262073/98‑1 (MQ=255)
                                                           aaaaTCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGa                                          >  2:311204/1‑100 (MQ=255)
                                                           aaaaTCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGa                                          >  2:2750974/1‑100 (MQ=255)
                                                              aTCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGt                                       <  2:1180158/100‑1 (MQ=255)
                                                                 gTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAGATAGGCCTGGTGCTTCAGTACTGAGCGGAGGTTCTTACAGCTGTAGGCAGAAGTTta                                    >  1:1304490/1‑100 (MQ=255)
                                                                  ttAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTtat                                   <  2:3023057/100‑1 (MQ=255)
                                                                   tAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTtata                                  <  2:830848/100‑1 (MQ=255)
                                                                    aaTCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTatat                                 <  1:1004861/100‑1 (MQ=255)
                                                                        aGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTta                             <  1:3115817/100‑1 (MQ=255)
                                                                        aGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATATGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTta                             <  2:307250/100‑1 (MQ=255)
                                                                           aTGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAg                          <  2:1786328/100‑1 (MQ=255)
                                                                           aTGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAg                          >  1:3197802/1‑100 (MQ=255)
                                                                            tGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGt                         >  2:3432287/1‑100 (MQ=255)
                                                                             gCTTACTATAACAGGCCGGGGGCCAGCGTGAATAGGACTGGTGCTTCAGTACTCAGCGGAGTTTCTTACAGCTGTAGGCAGAGGTTTATGTCTTATAg                          >  1:2593921/1‑98 (MQ=255)
                                                                             gCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAg                          <  1:1193678/98‑1 (MQ=255)
                                                                              cTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTAt                       >  2:3173873/1‑100 (MQ=255)
                                                                               ttAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATg                      >  1:3503881/1‑100 (MQ=255)
                                                                                tAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATGt                     >  2:1408681/1‑100 (MQ=255)
                                                                                tAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATGt                     <  2:3359209/100‑1 (MQ=255)
                                                                                tAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATGt                     >  2:438434/1‑100 (MQ=255)
                                                                                tAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGATTCAGTACTGAGCGGAGTGTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATGt                     >  1:1032177/1‑100 (MQ=255)
                                                                                  atatAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATGTAg                   <  2:132156/100‑1 (MQ=255)
                                                                                   tataACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATGTAGt                  >  2:1442621/1‑100 (MQ=255)
                                                                                   tataACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATGTAGt                  >  1:3380964/1‑100 (MQ=255)
                                                                                      aaCAGGCTGGTGAAAAGAGGAAATAGGTCTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATGTAGttta               >  1:3935346/1‑100 (MQ=255)
                                                                                      aaCAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATGTAGttta               <  1:903356/100‑1 (MQ=255)
                                                                                       aCAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATGTAGtttaa              <  2:3960421/100‑1 (MQ=255)
                                                                                         aGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATGTAGtttaatt            <  2:684297/100‑1 (MQ=255)
                                                                                         aGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTATTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATGTAGtttaatt            <  2:4066848/100‑1 (MQ=255)
                                                                                          ggCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATGTAGtttaattt           >  2:1283146/1‑100 (MQ=255)
                                                                                          ggCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATGTAGtttaattt           >  2:1996103/1‑100 (MQ=255)
                                                                                          ggCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATGTAGtttaattt           >  1:1394457/1‑100 (MQ=255)
                                                                                              ggAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATGTAGTTTAATTTAATc       >  1:1222027/1‑100 (MQ=255)
                                                                                              ggAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATGTAGTTTAATTTAATc       >  2:3941810/1‑100 (MQ=255)
                                                                                               gagaAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATGTAGTTTAATTTAATcc      >  1:2073853/1‑100 (MQ=255)
                                                                                                   aaaGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATGTAGTTTAATTTAATCCAtgt  <  2:2456229/100‑1 (MQ=255)
                                                                                                   aaaGAGGAAATAGGACTAGTGCTTCAGTACTGAGCGGATTTTCTTAAAGCTGTAGGCAGAAGTTTATATCTTATAGTATGTAGTTTAATTTAATCCacg   <  1:167891/99‑3 (MQ=255)
                                                                                                   |                                                                                                   
TTTTGTTTTTCACTGTGGTCGCATGCATTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATGTAGTTTAATTTAATCCATGT  >  NC_000913/1583932‑1584130

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: