Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,971,154 G→A 16.0% intergenic (+21/‑117) ygdQ → / → ygdR UPF0053 family inner membrane protein/DUF903 family verified lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,971,1540GA16.0% 471.3 / 63.8 193intergenic (+21/‑117)ygdQ/ygdRUPF0053 family inner membrane protein/DUF903 family verified lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (79/82);  new base A (17/14);  total (96/97)
Fisher's exact test for biased strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 202 total aligned reads displayed.

TTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGC  >  NC_000913/2971055‑2971253
                                                                                                   |                                                                                                   
ttACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCg                                                                                                     <  1:484263/100‑1 (MQ=255)
ttACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCg                                                                                                     <  1:851449/100‑1 (MQ=255)
ttACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTCAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCg                                                                                                     <  2:678624/100‑1 (MQ=255)
 tACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCgg                                                                                                    <  1:1398296/100‑1 (MQ=255)
 tACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCgg                                                                                                    <  2:2981991/100‑1 (MQ=255)
     tCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAgg                                                                                                >  1:815766/1‑100 (MQ=255)
     tCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAgg                                                                                                <  2:2131439/100‑1 (MQ=255)
     tCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAgg                                                                                                <  2:2592856/100‑1 (MQ=255)
      cTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAggg                                                                                               >  1:3778068/1‑100 (MQ=255)
      cTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAggg                                                                                               <  1:546576/100‑1 (MQ=255)
        aCTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCg                                                                                             <  1:3920122/100‑1 (MQ=255)
        aCTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCCCTCCCGCAGGAGAGCg                                                                                             <  1:1297752/100‑1 (MQ=255)
         cTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCACAACAAAAAGAATCCGCTCTGATACTCCGTACGCCCTCCCGCAGGAGAGCGt                                                                                            <  2:2496043/100‑1 (MQ=255)
           tCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCGACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCCCTCCCGCAGGAGAGCGTac                                                                                          <  1:1137116/100‑1 (MQ=38)
           tCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTac                                                                                          <  1:2559130/100‑1 (MQ=255)
              cgATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTactca                                                                                       <  2:1494042/100‑1 (MQ=255)
               gATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCACTCTGATACTCCGTACGCCCTCCCGCAGGAGAGCGTacgga                                                                                       <  2:943791/99‑4 (MQ=37)
                  gTTCTTCTCTATTGCGGTTGAATGCCTCAACTTGAGGCGCAACAAAAGGAATCCGCTCTGATACTCCGTACGCTCTCCCGCAAGAGAGCGTac                                                                                          <  2:1699464/93‑1 (MQ=255)
                     cttctCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCttt                                                                                >  2:2348321/1‑100 (MQ=255)
                     cttctCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCttt                                                                                >  2:1378576/1‑100 (MQ=255)
                      ttctCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCtttt                                                                               <  1:2956303/100‑1 (MQ=255)
                         tctATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTAcc                                                                            <  2:1136452/100‑1 (MQ=255)
                          ctATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTAccc                                                                           <  2:2786825/100‑1 (MQ=255)
                          ctATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTAccc                                                                           <  1:563708/100‑1 (MQ=255)
                          ctATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTAccc                                                                           <  1:3302422/100‑1 (MQ=255)
                          ctATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTAccc                                                                           <  1:1785502/100‑1 (MQ=255)
                          ctATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTAcc                                                                            <  1:1020595/99‑1 (MQ=255)
                           tATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCt                                                                          >  2:710673/1‑100 (MQ=255)
                             ttGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTca                                                                        <  2:1731791/100‑1 (MQ=255)
                               gCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTcaca                                                                      >  1:3828285/1‑100 (MQ=255)
                                cGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAg                                                                     <  1:1061018/100‑1 (MQ=255)
                                     naaaGCCTCAACTTGATTCGCAACAAAAAGAATCCGATCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTaa                                                                >  1:2779182/2‑100 (MQ=255)
                                     gAAAGCCTCAACTTGATTCGCAACAAAGAGAATCCGCTCTGATACTCCGTACACTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTaa                                                                <  2:2219457/100‑1 (MQ=255)
                                     gAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTaa                                                                >  2:1271351/1‑100 (MQ=255)
                                     gAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTaa                                                                <  1:3143647/100‑1 (MQ=255)
                                     gAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTaa                                                                <  1:1527511/100‑1 (MQ=255)
                                          ccTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGAttt                                                           <  2:336726/100‑1 (MQ=255)
                                            tCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTc                                                         >  1:1859410/1‑100 (MQ=255)
                                            tCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTc                                                         >  2:175739/1‑100 (MQ=255)
                                                cTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGc                                                     <  1:2465462/100‑1 (MQ=255)
                                                    aTTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTc                                                 >  1:437970/1‑100 (MQ=255)
                                                     ttCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCa                                                <  2:1866148/100‑1 (MQ=255)
                                                      tCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCaa                                               >  2:856105/1‑100 (MQ=255)
                                                        gCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGa                                             >  1:1333376/1‑100 (MQ=255)
                                                        gCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGa                                             >  1:3533769/1‑100 (MQ=255)
                                                          aacaaAAAGAATCCGCTCTGATACTCCGTCCGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGata                                           <  2:1610240/100‑1 (MQ=255)
                                                          aacaaAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGata                                           <  2:2619121/100‑1 (MQ=255)
                                                          aacaaAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGata                                           >  1:2107544/1‑100 (MQ=255)
                                                          aacaaAAAGAATCCGCTCTGATACGCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGat                                            <  1:3718275/99‑1 (MQ=255)
                                                           agaaaaaGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGatat                                          >  1:2672264/3‑100 (MQ=255)
                                                            caaAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGatata                                         <  2:1025066/100‑1 (MQ=255)
                                                            caaAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGatata                                         >  1:2491670/1‑100 (MQ=255)
                                                            caaAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGatata                                         <  1:633185/100‑1 (MQ=255)
                                                            caaAAAGAATCCGCTCTGATACTCCGTACGCCCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCATCACAGTTTAAGATTTTCCTGCTTTCATGatata                                         >  2:267581/1‑100 (MQ=255)
                                                              aaaaGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGAtatata                                       >  1:2387214/1‑100 (MQ=255)
                                                              aaaaGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCATACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGAtatata                                       >  2:1591717/1‑100 (MQ=255)
                                                               cgctcttCCGATCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATaa                                      <  2:1040052/94‑1 (MQ=255)
                                                               aaaGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATaa                                      <  2:3113497/100‑1 (MQ=255)
                                                                aaGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAAc                                     >  2:2717590/1‑100 (MQ=255)
                                                                    aTCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCg                                 <  2:1347777/100‑1 (MQ=255)
                                                                    aTCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGAGTACTCACTCGTTTTACCCTCACAGTTTAAGATTGTCCTGCTTTCGAGATATATAACGt                                   >  2:2971696/1‑98 (MQ=255)
                                                                       cGCTCTGCTACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGtt                              <  2:281969/100‑1 (MQ=255)
                                                                       cGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGtt                              <  2:2567179/100‑1 (MQ=255)
                                                                       cGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGtt                              >  1:477303/1‑100 (MQ=255)
                                                                       cGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGt                               >  2:3578490/1‑99 (MQ=255)
                                                                        gCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGttt                             <  2:3828256/100‑1 (MQ=255)
                                                                         ctctGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTta                            <  1:3041837/100‑1 (MQ=255)
                                                                           ctGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTtata                          >  1:457788/1‑100 (MQ=255)
                                                                               tACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACa                      <  2:2990653/100‑1 (MQ=255)
                                                                               tACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACa                      <  1:3491591/100‑1 (MQ=255)
                                                                               tACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCGCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACa                      <  1:160443/100‑1 (MQ=255)
                                                                                aCTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAg                     >  2:20633/1‑100 (MQ=255)
                                                                                 gtgagTACGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCCCAGTCCAAGATTTTCCTGCTTTCAAGATATATAACGTCCGTCTATAAACAGa                    >  2:483258/5‑100 (MQ=21)
                                                                                 gtgagTACGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGa                    >  1:665286/5‑100 (MQ=38)
                                                                                 cTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGa                    >  1:3399714/1‑100 (MQ=255)
                                                                                  tgattACGCTCTCCTGAGTGAGGGCGTACTCACTCTTTATAAAATAACAGTTTAAGATTTACCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGAc                   <  2:3710620/96‑1 (MQ=255)
                                                                                  tgagTACGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCCCCGTTTACGCTTTTGCTGCGTTCAAGATATATAACGTGCGTTTACAAAc                       >  2:3428878/4‑96 (MQ=16)
                                                                                  tCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGAc                   <  2:2387202/100‑1 (MQ=255)
                                                                                   ccGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACt                  <  2:3098036/100‑1 (MQ=255)
                                                                                   ccGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACt                  >  1:2184033/1‑100 (MQ=255)
                                                                                    cGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACta                 >  2:455284/1‑100 (MQ=255)
                                                                                    cGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATACATAACCTCGAGTTATAAACAGACta                 >  2:2266989/1‑100 (MQ=255)
                                                                                    agTACGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACta                 >  1:269774/2‑100 (MQ=38)
                                                                                    agTACGCCCTCCCGCAGGAAAGCGAACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTGTATAAACAGACta                 >  2:3478993/2‑100 (MQ=255)
                                                                                     gTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACtat                <  1:1721469/100‑1 (MQ=255)
                                                                                     gTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACtat                >  1:2089211/1‑100 (MQ=255)
                                                                                     gTACGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACtat                >  2:1585847/1‑100 (MQ=38)
                                                                                     gTACGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACtat                >  2:581480/1‑100 (MQ=38)
                                                                                     gTACGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACtat                >  1:2981664/1‑100 (MQ=38)
                                                                                      tACGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACtatt               >  2:3229077/1‑100 (MQ=255)
                                                                                       aCGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACtatta              >  1:1377095/1‑100 (MQ=255)
                                                                                        cGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACtattat             >  1:1793805/1‑100 (MQ=255)
                                                                                         gCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATc            <  2:153542/100‑1 (MQ=255)
                                                                                          ctctCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCa           >  1:2287177/1‑100 (MQ=255)
                                                                                           tctcCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAt          <  2:2000357/100‑1 (MQ=255)
                                                                                              ccTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAtagg       >  1:933962/1‑100 (MQ=255)
                                                                                              ccTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAtagg       >  2:2457562/1‑100 (MQ=255)
                                                                                               cTGCGGGAGGGCGTACTCCCTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAtaggt      <  2:1983183/100‑1 (MQ=255)
                                                                                                tGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAtaggta     <  2:1827426/100‑1 (MQ=255)
                                                                                                   gggAGGGCGTACTGACTCCTTGTACCCTCACAGTTTAAGATTTTCCTGGTTTCAAGATATAGAACGTCGGTTTATAAACAGACTATTATCAtaggtagg   <  2:3931607/99‑1 (MQ=255)
                                                                                                   |                                                                                                   
TTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGC  >  NC_000913/2971055‑2971253

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: