Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,971,159 G→A 16.2% intergenic (+26/‑112) 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,1590GA16.2% 473.5 / 57.1 192intergenic (+26/‑112)ygdQ/ygdRUPF0053 family inner membrane protein/DUF903 family verified lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (75/85);  new base A (17/15);  total (92/100)
Fisher's exact test for biased strand distribution p-value = 5.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 200 total aligned reads displayed.

TCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAAT  >  NC_000913/2971060‑2971258
                                                                                                   |                                                                                                   
 cTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAggg                                                                                                    >  1:3778068/1‑100 (MQ=255)
   aCTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCg                                                                                                  <  1:3920122/100‑1 (MQ=255)
   aCTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCCCTCCCGCAGGAGAGCg                                                                                                  <  1:1297752/100‑1 (MQ=255)
    cTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCCCTCCCGCAGGAGAGCGt                                                                                                 <  1:209962/100‑1 (MQ=255)
      tCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCGACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCCCTCCCGCAGGAGAGCGTac                                                                                               <  1:1137116/100‑1 (MQ=38)
      tCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTac                                                                                               <  1:563670/100‑1 (MQ=255)
      tCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTac                                                                                               <  1:2559130/100‑1 (MQ=255)
      tCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGAGTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCCCTCCCGCAGGAGAGCGTac                                                                                               <  2:724315/100‑1 (MQ=38)
        gcgTTGTTCTTCTCTATTGCGGTTGAATGCCTCACTTTGATTCGCAACAAACAGAATCCGCTCTGATATTTCTTACGCCCTCTCGCAGGAGAGCGTac                                                                                               <  2:1355254/98‑1 (MQ=255)
        gcgATGTTCTTCTCTATTGCGGTTGAAAGCCTCGACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTactc                                                                                             >  2:2569019/1‑100 (MQ=255)
         cgATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTactca                                                                                            >  2:3104824/1‑100 (MQ=255)
         cgATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTactca                                                                                            <  2:1494042/100‑1 (MQ=255)
               tcttctCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCCCTCCCGCAGGAGAGCGTacggagtatca                                                                                      <  1:3138135/100‑10 (MQ=25)
                cttctCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCCCTCCCGCAGGAGAGCGTacggagtatcag                                                                                     <  2:2377644/100‑11 (MQ=25)
                cttctCTATTGCGGTTGAAAGCCTCAACGTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCttt                                                                                     <  2:2998005/100‑1 (MQ=255)
                  tctcTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTa                                                                                   <  1:1547264/100‑1 (MQ=255)
                    tctATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTAcc                                                                                 <  2:2486281/100‑1 (MQ=255)
                    tctATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTAcc                                                                                 <  2:1136452/100‑1 (MQ=255)
                     ctATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTAccc                                                                                <  2:2786825/100‑1 (MQ=255)
                     ctATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTAccc                                                                                <  1:1785501/100‑1 (MQ=255)
                     ctATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTAccc                                                                                >  2:3195283/1‑100 (MQ=255)
                     ctATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTAcc                                                                                 <  1:1020595/99‑1 (MQ=255)
                       aTTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTc                                                                              <  2:1939521/100‑1 (MQ=255)
                        ttGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTca                                                                             <  2:1740346/100‑1 (MQ=255)
                        ttGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTca                                                                             >  2:69928/1‑100 (MQ=255)
                         tGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTcac                                                                            >  2:938065/1‑100 (MQ=255)
                         tGCGGTTGAAAGCCTCAACTTGATTCACAACACAAAGAATCCGCTCTGCTACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTcac                                                                            >  2:2607848/1‑100 (MQ=255)
                          gCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTcaca                                                                           >  1:3828285/1‑100 (MQ=255)
                           cGGTTGAAAGCCTCAACTTGATTCGCAACAAGAAGCAGCCGCTCTGATACTCCGTACGCTCTCGTGCGGGAGGGCGTACTTACTGCTTTTACCCTCACAg                                                                          <  2:227374/100‑1 (MQ=255)
                           cGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAg                                                                          <  1:1061018/100‑1 (MQ=255)
                               tGAAAGCCTCAACTTGATTCGCAACAAAGAGAATCCGATCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGt                                                                         >  1:1503229/1‑97 (MQ=255)
                               tGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTa                                                                      >  2:2638939/1‑100 (MQ=255)
                                gAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTaa                                                                     <  1:3143647/100‑1 (MQ=255)
                                 aaaGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAg                                                                    <  1:813452/100‑1 (MQ=255)
                                     ccTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGAttt                                                                >  2:1989083/1‑100 (MQ=255)
                                     ccTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGAttt                                                                <  2:336726/100‑1 (MQ=255)
                                           cTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGc                                                          <  1:2465462/100‑1 (MQ=255)
                                           cTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGc                                                          >  1:2383147/1‑100 (MQ=255)
                                               aTTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTc                                                      <  2:950708/100‑1 (MQ=255)
                                               aTTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTc                                                      <  1:2234468/100‑1 (MQ=255)
                                               aTTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTc                                                      >  1:437970/1‑100 (MQ=255)
                                                ttCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCa                                                     >  1:846146/1‑100 (MQ=255)
                                                  cGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAg                                                   >  2:2223449/1‑100 (MQ=255)
                                                  cGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAg                                                   <  2:1588667/100‑1 (MQ=255)
                                                   gCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGa                                                  >  1:3533769/1‑100 (MQ=255)
                                                   gCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGa                                                  >  1:1333376/1‑100 (MQ=255)
                                                    caacaaAAAGAATCCGCTCTGATACTCCGCACGCTCTCCTGCGGCAGGGCGTACTCACTCCTTTTACCCTCCCCGCTTAACATTTTCCTCCTTTCa                                                     >  2:2058843/1‑96 (MQ=255)
                                                     aacaaAAAGAATCCGCTCTGATACTCCGTCCGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGata                                                <  2:1610240/100‑1 (MQ=255)
                                                     aacaaAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGata                                                <  2:1702640/100‑1 (MQ=255)
                                                     aacaaAAAGAATCCGCTCTGATACGCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGat                                                 <  1:3718275/99‑1 (MQ=255)
                                                       caaAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGatata                                              <  2:1025066/100‑1 (MQ=255)
                                                       caaAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGatata                                              >  1:2685545/1‑100 (MQ=255)
                                                       caaAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGatata                                              >  1:2491670/1‑100 (MQ=255)
                                                       caaAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTATTACCCTCACAGTTTAAGATTTTCCTGTTTTCAAGttata                                              >  2:2712536/1‑100 (MQ=255)
                                                          cgctcttCCGATCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATaa                                           <  2:1040052/94‑1 (MQ=255)
                                                            aGACTCCTCTCTGTTACTCCGTACGCTCTCCTTCGGGCGGGCTTACTTACTTCTTTTTCCCACACAGTTTGAGATGTTCCTGCTTTCAAGatatat                                             <  2:2142650/96‑1 (MQ=255)
                                                              cttCCGATCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTc                                       <  1:1854115/98‑1 (MQ=255)
                                                              aaTCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTc                                       <  1:634583/100‑1 (MQ=255)
                                                               aTCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGAGTACTCACTCGTTTTACCCTCACAGTTTAAGATTGTCCTGCTTTCGAGATATATAACGt                                        >  2:2971696/1‑98 (MQ=255)
                                                                 ccGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGt                                    >  1:3038243/1‑100 (MQ=255)
                                                                  cGCTCTGCTACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGtt                                   <  2:281969/100‑1 (MQ=255)
                                                                  cGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGtt                                   <  1:217961/100‑1 (MQ=255)
                                                                  cGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGtt                                   >  1:477303/1‑100 (MQ=255)
                                                                  cGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGtt                                   >  2:1780614/1‑100 (MQ=255)
                                                                  cGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGt                                    >  2:3578490/1‑99 (MQ=255)
                                                                   gCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGttt                                  <  2:3828256/100‑1 (MQ=255)
                                                                      ctGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTtata                               >  1:457788/1‑100 (MQ=255)
                                                                        gATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTcagatcggaag                             >  2:1854106/1‑90 (MQ=255)
                                                                          tACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACa                           <  1:3491591/100‑1 (MQ=255)
                                                                           aCTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAg                          >  1:179082/1‑100 (MQ=255)
                                                                            gtgagTACGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCCCAGTCCAAGATTTTCCTGCTTTCAAGATATATAACGTCCGTCTATAAACAGa                         >  2:483258/5‑100 (MQ=21)
                                                                            gtgagTACGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGa                         >  1:564820/5‑100 (MQ=37)
                                                                            cTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGa                         >  1:3399714/1‑100 (MQ=255)
                                                                            cTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGa                         >  1:2710526/1‑100 (MQ=255)
                                                                             tgattACGCTCTCCTGAGTGAGGGCGTACTCACTCTTTATAAAATAACAGTTTAAGATTTACCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGAc                        <  2:3710620/96‑1 (MQ=255)
                                                                             tgagTACGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCCCCGTTTACGCTTTTGCTGCGTTCAAGATATATAACGTGCGTTTACAAAc                            >  2:3428878/4‑96 (MQ=16)
                                                                              gagTACGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACt                       >  1:1964738/3‑100 (MQ=38)
                                                                               cGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACta                      >  1:2960266/1‑100 (MQ=255)
                                                                               cGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACta                      >  2:455284/1‑100 (MQ=255)
                                                                               agTACGCCCTCCCGCAGGAAAGCGAACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTGTATAAACAGACta                      >  2:3478993/2‑100 (MQ=255)
                                                                                gTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCTGTTTATAAACCGAcgtt                     <  2:2985411/100‑4 (MQ=255)
                                                                                gTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACtat                     <  1:1721469/100‑1 (MQ=255)
                                                                                gTACGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACtat                     >  2:1585847/1‑100 (MQ=38)
                                                                                gTACGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACtat                     >  2:581480/1‑100 (MQ=38)
                                                                                gTACGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACtat                     >  1:2797894/1‑100 (MQ=38)
                                                                                gTACGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATCTTCCTGCCTTCAAGATATATAACGTCGGTTTATAAACAGACtat                     >  2:211193/1‑100 (MQ=38)
                                                                                  aCGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACtatta                   >  1:1377095/1‑100 (MQ=255)
                                                                                    gCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATc                 <  2:153542/100‑1 (MQ=255)
                                                                                     ctctCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTGTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCa                >  2:2447388/1‑100 (MQ=255)
                                                                                      tctcCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCa                <  2:1294015/99‑1 (MQ=255)
                                                                                       ctcCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAta              <  1:518860/100‑1 (MQ=255)
                                                                                        tcCTGCGGGAGGGCGTACTCTCTGCTTTTGTCCTCACAGTTTGAGATTTTCCTGCTTTCAAGATATATACCGTCGGTTTATAAACAGACTATTATCAt               <  2:2547758/98‑1 (MQ=255)
                                                                                         ccTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAtagg            >  2:1270803/1‑100 (MQ=255)
                                                                                         ccTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAtagg            >  1:933962/1‑100 (MQ=255)
                                                                                              gggAGGGCGTACTGACTCCTTGTACCCTCACAGTTTAAGATTTTCCTGGTTTCAAGATATAGAACGTCGGTTTATAAACAGACTATTATCAtaggtagg        <  2:3931607/99‑1 (MQ=255)
                                                                                               ggAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGcc      <  2:3234962/100‑1 (MQ=255)
                                                                                               ggAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGcc      <  1:1623704/100‑1 (MQ=255)
                                                                                               ggAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGcc      <  1:3304749/100‑1 (MQ=255)
                                                                                               ggAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGcc      <  1:1222913/100‑1 (MQ=255)
                                                                                                   ggCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAAt  >  1:239380/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
TCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAAT  >  NC_000913/2971060‑2971258

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