Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,971,151 T→C 15.4% intergenic (+18/‑120) 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,1510TC15.4% 477.5 / 66.9 199intergenic (+18/‑120)ygdQ/ygdRUPF0053 family inner membrane protein/DUF903 family verified lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (82/86);  new base C (17/14);  total (99/100)
Fisher's exact test for biased strand distribution p-value = 5.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.56e-01

Only 100 of 201 total aligned reads displayed.

GTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTA  >  NC_000913/2971054‑2971250
                                                                                                 |                                                                                                   
gTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTg                                                                                                    >  2:3827469/1‑99 (MQ=255)
gTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGc                                                                                                   >  2:2559113/1‑100 (MQ=255)
 ttACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCg                                                                                                  >  2:2727826/1‑100 (MQ=255)
 ttACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTCAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCg                                                                                                  <  2:678624/100‑1 (MQ=255)
      tCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAgg                                                                                             >  1:3508551/1‑100 (MQ=255)
      tCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAgg                                                                                             <  2:2131439/100‑1 (MQ=255)
        tACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGc                                                                                           >  2:1180987/1‑100 (MQ=255)
          cTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCACAACAAAAAGAATCCGCTCTGATACTCCGTACGCCCTCCCGCAGGAGAGCGt                                                                                         <  2:2496043/100‑1 (MQ=255)
            tCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCGACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCCCTCCCGCAGGAGAGCGTac                                                                                       <  1:1137116/100‑1 (MQ=38)
            tCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTac                                                                                       <  1:2559130/100‑1 (MQ=255)
            tCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTac                                                                                       >  2:2804249/1‑100 (MQ=255)
            tCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTac                                                                                       <  1:563670/100‑1 (MQ=255)
            tCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCCCTCCCGCAGGAGAGCGTac                                                                                       <  1:467353/100‑1 (MQ=38)
              gcgATGTTCTTCTCTATTGCGGTTGAAAGCCTCGACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTactc                                                                                     >  2:2569019/1‑100 (MQ=255)
               cgATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTactca                                                                                    >  2:3104824/1‑100 (MQ=255)
               cgATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTactca                                                                                    <  2:1494042/100‑1 (MQ=255)
                gATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCACTCTGATACTCCGTACGCCCTCCCGCAGGAGAGCGTacgga                                                                                    <  2:943791/99‑4 (MQ=37)
                   gTTCTTCTCTATTGCGGTTGAATGCCTCAACTTGAGGCGCAACAAAAGGAATCCGCTCTGATACTCCGTACGCTCTCCCGCAAGAGAGCGTac                                                                                       <  2:1699464/93‑1 (MQ=255)
                   gTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTcc                                                                                >  1:1301508/1‑100 (MQ=255)
                      cttctCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCttt                                                                             >  2:1378576/1‑100 (MQ=255)
                      cttctCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCttt                                                                             <  1:3780202/100‑1 (MQ=255)
                      cttctCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCttt                                                                             >  2:2348321/1‑100 (MQ=255)
                      cttctCTATTGCGGTTGAAAGCCTCAACGTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCttt                                                                             <  2:2998005/100‑1 (MQ=255)
                       ttctCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCtttt                                                                            <  1:2956303/100‑1 (MQ=255)
                        tctcTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTa                                                                           <  1:1547264/100‑1 (MQ=255)
                           ctATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTAccc                                                                        <  1:1785502/100‑1 (MQ=255)
                           ctATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTAccc                                                                        >  2:3195283/1‑100 (MQ=255)
                           ctATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTAcc                                                                         <  1:1020595/99‑1 (MQ=255)
                            tATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCt                                                                       >  2:710673/1‑100 (MQ=255)
                              ttGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTca                                                                     <  2:1731791/100‑1 (MQ=255)
                               tGCGGTTGAAAGCCTCAACTTGATTCACAACACAAAGAATCCGCTCTGCTACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTcac                                                                    >  2:2607848/1‑100 (MQ=255)
                                 cGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAg                                                                  <  1:1061018/100‑1 (MQ=255)
                                     tGAAAGCCTCAACTTGATTCGCAACAAAGAGAATCCGATCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGt                                                                 >  1:1503229/1‑97 (MQ=255)
                                     tGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTa                                                              >  2:2638939/1‑100 (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)
                                       aaaGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAg                                                            <  1:813452/100‑1 (MQ=255)
                                         aGCGTCAACTTGAGTGGCAACAAAAAGAAGCCGCTCGAATACCCCGTACGCTCTCCTGCGGGAGAGCTTACCAACTCCTGGTACCCTCACAGTGTAAg                                                            <  2:3880551/98‑1 (MQ=255)
                                          gCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGAtt                                                         >  1:3677595/1‑100 (MQ=255)
                                             tCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTc                                                      >  1:1859410/1‑100 (MQ=255)
                                             tCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTc                                                      >  2:175739/1‑100 (MQ=255)
                                               aaCTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCt                                                    >  1:415559/1‑100 (MQ=255)
                                                 cTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGc                                                  <  1:2465462/100‑1 (MQ=255)
                                                      ttCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCa                                             >  1:846146/1‑100 (MQ=255)
                                                      ttCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCa                                             <  2:1866148/100‑1 (MQ=255)
                                                       tCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCaa                                            >  2:856105/1‑100 (MQ=255)
                                                        cGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCTGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAg                                           <  2:3270245/100‑1 (MQ=255)
                                                        cGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAg                                           <  1:3326262/100‑1 (MQ=255)
                                                         gCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGa                                          >  1:1333464/1‑100 (MQ=255)
                                                           aacaaAAAGAATCCGCTCTGATACTCCGTCCGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGata                                        <  2:1610240/100‑1 (MQ=255)
                                                           aacaaAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGata                                        >  1:2107544/1‑100 (MQ=255)
                                                           aacaaAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGata                                        >  2:298081/1‑100 (MQ=255)
                                                            agaaaaaGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGatat                                       >  1:2672264/3‑100 (MQ=255)
                                                             caaAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGatata                                      >  1:2491670/1‑100 (MQ=255)
                                                             caaAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTATTACCCTCACAGTTTAAGATTTTCCTGTTTTCAAGttata                                      >  2:2712536/1‑100 (MQ=255)
                                                              aaaaaGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCTTCACAGTTTAAGATTTTCCTGCTTTCAAGatatat                                     >  1:3418220/1‑100 (MQ=255)
                                                               aaaaGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGAtatata                                    <  1:3885026/100‑1 (MQ=255)
                                                               aaaaGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGAtatata                                    >  1:2387214/1‑100 (MQ=255)
                                                               aaaaGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCATACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGAtatata                                    >  2:1591717/1‑100 (MQ=255)
                                                                 aaGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAAc                                  <  2:1092067/100‑1 (MQ=255)
                                                                    aaTCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTc                               <  1:634583/100‑1 (MQ=255)
                                                                     aTCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCg                              <  2:1347777/100‑1 (MQ=255)
                                                                      tCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCgg                             <  1:1378585/100‑1 (MQ=255)
                                                                        cGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGt                            >  2:3578490/1‑99 (MQ=255)
                                                                          ctctGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTta                         <  1:3041837/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:564820/5‑100 (MQ=37)
                                                                                   tgattACGCTCTCCTGAGTGAGGGCGTACTCACTCTTTATAAAATAACAGTTTAAGATTTACCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGAc                <  2:3710620/96‑1 (MQ=255)
                                                                                   tCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGAc                <  2:2387202/100‑1 (MQ=255)
                                                                                    gagTACGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACt               >  1:444551/3‑100 (MQ=38)
                                                                                    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)
                                                                                      gTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCTGTTTATAAACCGAcgtt             <  2:2985411/100‑4 (MQ=255)
                                                                                      gTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACtat             >  1:92798/1‑100 (MQ=255)
                                                                                      gTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACtat             >  1:2089211/1‑100 (MQ=255)
                                                                                      gTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACtat             <  1:1721469/100‑1 (MQ=255)
                                                                                      gTACGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACtat             >  1:2981664/1‑100 (MQ=38)
                                                                                      gTACGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACtat             >  2:1585847/1‑100 (MQ=38)
                                                                                      gTACGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACtat             >  2:581480/1‑100 (MQ=38)
                                                                                       tACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACtatt            >  1:95047/1‑100 (MQ=255)
                                                                                         cGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACtattat          <  2:2846603/100‑1 (MQ=255)
                                                                                         cGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACtattat          >  1:1793805/1‑100 (MQ=255)
                                                                                          gCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATc         <  2:153542/100‑1 (MQ=255)
                                                                                           ctctCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCa        <  2:3399691/100‑1 (MQ=255)
                                                                                           ctctCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCa        >  1:2287177/1‑100 (MQ=255)
                                                                                            tctcCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAt       <  1:3765704/100‑1 (MQ=255)
                                                                                            tctcCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAt       <  2:2000357/100‑1 (MQ=255)
                                                                                             ctcCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAta      <  1:518860/100‑1 (MQ=255)
                                                                                               ccTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAtagg    >  2:2457562/1‑100 (MQ=255)
                                                                                                cTGCGGGAGGGCGTACTCCCTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAtaggt   <  2:1983183/100‑1 (MQ=255)
                                                                                                cTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAtaggt   >  2:1031623/1‑100 (MQ=255)
                                                                                                 tGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAtaggta  <  2:1827426/100‑1 (MQ=255)
                                                                                                 tGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAtaggta  >  1:478945/1‑100 (MQ=255)
                                                                                                 |                                                                                                   
GTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTA  >  NC_000913/2971054‑2971250

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: