Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,971,159 G→A 16.0% 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.0% 429.5 / 59.7 194intergenic (+26/‑112)ygdQ/ygdRUPF0053 family inner membrane protein/DUF903 family verified lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (87/73);  new base A (15/17);  total (102/92)
Fisher's exact test for biased strand distribution p-value = 4.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.57e-02

Only 100 of 201 total aligned reads displayed.

TCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAAT  >  NC_000913/2971060‑2971258
                                                                                                   |                                                                                                   
tCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAgg                                                                                                     <  2:2561938/100‑1 (MQ=255)
  tACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGc                                                                                                   >  1:40270/1‑100 (MQ=255)
   aCTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCg                                                                                                  <  1:2698294/100‑1 (MQ=255)
    cTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCCCTCCCGCAGGAGAGCGt                                                                                                 <  2:1373818/100‑1 (MQ=255)
     ttCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTa                                                                                                <  1:1316746/100‑1 (MQ=255)
     ttCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCCCTCCCGCAGGAGAGCGTa                                                                                                <  1:3016590/100‑1 (MQ=255)
      tCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTTATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCCCTCCCGCAGGAGAGCGTac                                                                                               <  2:218175/100‑1 (MQ=38)
      tCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTTATACTCCGTACGCCCTCCCGCAGGAGAGCGTac                                                                                               <  2:287073/100‑1 (MQ=38)
      tCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTTTGGGAGGGCGttc                                                                                               <  1:2657885/100‑3 (MQ=255)
      tCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTac                                                                                               <  1:3799787/100‑1 (MQ=255)
       cgcgATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCCCTCCCGCAGGAGAGCGTacg                                                                                              <  1:3749195/100‑2 (MQ=38)
         cgATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTactca                                                                                            >  1:3161926/1‑100 (MQ=255)
         cgATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCCCTCCCGCAGGAGAGCGTacgga                                                                                            <  2:3354597/100‑4 (MQ=38)
           aTGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTactcact                                                                                          <  1:958564/100‑1 (MQ=255)
           aTGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCCCTCCCGCAGGAGAGCGTacggagt                                                                                          <  1:2325711/100‑6 (MQ=37)
             gTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTTCGGGAGGGCGTACTCACTcc                                                                                        <  2:946849/100‑1 (MQ=255)
                cttctCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCttt                                                                                     >  1:1903772/1‑100 (MQ=255)
                cttctCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCttt                                                                                     >  2:1241982/1‑100 (MQ=255)
                 ttctCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCtttt                                                                                    >  1:1032001/1‑100 (MQ=255)
                    tctATTGCTGTTTAAAGCCTCAACTTGATTCGCAATAAAAAGAATCCGCTTTTATACTCCGTACGCTCTCCTGCTGGAGGGCGTACTCACTCCTTTTAcc                                                                                 <  2:366693/100‑1 (MQ=255)
                     ctATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTAGTCACTCCTTTTAccc                                                                                >  1:3875004/1‑100 (MQ=255)
                     ctATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTAccc                                                                                <  2:1975499/100‑1 (MQ=255)
                       aTTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTc                                                                              >  1:2421377/1‑100 (MQ=255)
                         tGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTcac                                                                            >  2:3154043/1‑100 (MQ=255)
                          gCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTcaca                                                                           >  2:1894685/1‑100 (MQ=255)
                          gCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTcaca                                                                           >  1:477564/1‑100 (MQ=255)
                           cGGTTGAAAGCCTCAACTTGATTCGCAGCAAAAAGAATCCGCTCTGATACTCCCTACGCTCTGCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAg                                                                          <  2:2080286/100‑1 (MQ=255)
                             gTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGtt                                                                        >  1:3327859/1‑100 (MQ=255)
                              ttGAATGCCTCAACTTTCTTCGCAACAAAAAGAATCCGCTCTTATACTCCGTACGCTCTCCTGCGGGAGGGCTTACTCACTCCTTTTACCCTCACAGtt                                                                        <  2:3460872/99‑1 (MQ=255)
                                 aaaGCCTCAACTTGATTCGCAACAAAAATAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAg                                                                    <  1:535307/100‑1 (MQ=255)
                                    gCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGAtt                                                                 >  1:2133344/1‑100 (MQ=255)
                                     ccTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGAttt                                                                >  1:1590259/1‑100 (MQ=255)
                                      cTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGAtttt                                                               >  2:3002306/1‑100 (MQ=255)
                                      cTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGAtttt                                                               >  1:1116781/1‑100 (MQ=255)
                                      cTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGAtttt                                                               >  1:1715175/1‑100 (MQ=255)
                                        cAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTcc                                                             >  2:3838351/1‑100 (MQ=255)
                                         aaCTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCt                                                            >  1:149708/1‑100 (MQ=255)
                                         aaCTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCt                                                            >  2:2702201/1‑100 (MQ=255)
                                         aaCTTGATTCGAAAAAAAAAGAATCCGCTCTGTTTCTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTcc                                                             <  1:574779/99‑1 (MQ=255)
                                            ttGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTAAAAGATTAAGATTTTca                                                             >  1:846012/1‑95 (MQ=255)
                                              gATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCttt                                                       <  2:1758482/100‑1 (MQ=255)
                                                  cGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAg                                                   >  2:1601764/1‑100 (MQ=255)
                                                  cGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAg                                                   >  2:192463/1‑100 (MQ=255)
                                                   gCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGa                                                  <  1:385636/100‑1 (MQ=255)
                                                   gCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGa                                                  >  2:2943654/1‑100 (MQ=255)
                                                    caacaaAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGat                                                 <  1:519172/100‑1 (MQ=255)
                                                     aacaaAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGata                                                <  1:3438704/100‑1 (MQ=255)
                                                       caaAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGatata                                              <  2:3878858/100‑1 (MQ=255)
                                                       caaAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGatata                                              >  2:167066/1‑100 (MQ=255)
                                                        aaaaaGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGatatat                                             <  1:1970551/100‑1 (MQ=255)
                                                        aaaaaGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGatatat                                             >  1:2237527/1‑100 (MQ=255)
                                                         aaaaGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGAtatata                                            <  1:2740892/100‑1 (MQ=255)
                                                         aaaaGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGAtatata                                            <  2:206537/100‑1 (MQ=255)
                                                         aaaaGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGAtatata                                            >  1:2074744/1‑100 (MQ=255)
                                                            ccgatctgCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACg                                         <  2:3584495/93‑1 (MQ=255)
                                                             gAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAAc                                          >  2:3804886/1‑98 (MQ=255)
                                                              aaTCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTc                                       >  2:3755345/1‑100 (MQ=255)
                                                                tCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCgg                                     <  1:1926751/100‑1 (MQ=255)
                                                                  cGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGtt                                   <  1:1638979/100‑1 (MQ=255)
                                                                   gCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACgagatcgg                                  >  1:3584534/1‑93 (MQ=255)
                                                                    ctctGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTta                                 <  1:1394669/100‑1 (MQ=255)
                                                                      ctGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTtata                               >  1:2380346/1‑100 (MQ=255)
                                                                      ctGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTtata                               >  2:1748373/1‑100 (MQ=255)
                                                                       tgtttcTCCGTACGCTCTCCTGCGGGAGGGCGTACTCAATCCTTTTACCCTCACAGTTTAAGATCTTGCTGATTTCAAGATATATAACGTCGGTTTATaa                              >  2:2513337/6‑100 (MQ=255)
                                                                       tGATACTCCGTACGCTCTCCTGCGGGAGTGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATaa                              <  1:346530/100‑1 (MQ=255)
                                                                          tACTCCGTACGCTCTCTTGCGGGAGGGCGTACTCACTTCTTTTTCCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACa                           <  1:3376511/100‑1 (MQ=255)
                                                                          tACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACa                           <  1:657802/100‑1 (MQ=255)
                                                                          tACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACa                           <  1:1755323/100‑1 (MQ=255)
                                                                          tACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATAGATAACGTCGGTTTATAAACa                           <  2:2395733/100‑1 (MQ=255)
                                                                          tACTCCGTACGCTCTCCTACGCGAGGGCGTACTCACTCCTTATACCATCCCAGCTTCAGATTTTCCTGCATTCAAAATATATAACGTCGGTTTATAAACa                           >  2:2689091/1‑100 (MQ=255)
                                                                          gagtgagTACGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACa                           >  1:1421750/7‑100 (MQ=37)
                                                                           aCTCCGTACGCTCTCCTGCGGGAGGGCGTACTTACTCCTTTTTCCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAg                          <  1:1207391/100‑1 (MQ=255)
                                                                            cTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGa                         >  2:1452497/1‑100 (MQ=255)
                                                                            cTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGa                         >  2:3136188/1‑100 (MQ=255)
                                                                              gagTACGCCCTCCCGCAGCAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAAATATTCATGCTTTCAAGATATATAACGTCGGTTTATAAACAGAc                        >  1:3225096/3‑99 (MQ=37)
                                                                              ccGTACGCTCTCCTGCGGGAGGGAGTACTCACTCCTTTTAACCTCACAGTTTAAGATTTTCCTGCTTTCTATATATATATCGTCGGTTTATAAACAGAc                        <  2:2759922/99‑1 (MQ=255)
                                                                                gTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACCGTTTAAGATTTTCCTGCTTTAAAGATCTATAACGTCGGTTTATAAACAGACtat                     >  2:3877678/1‑100 (MQ=255)
                                                                                gTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACtat                     >  1:3677825/1‑100 (MQ=255)
                                                                                gTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACta                      >  2:2524053/1‑99 (MQ=255)
                                                                                 tacgttattcTGCGGGAGGGCGTATTCACTACTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGAc                        <  1:2627477/87‑1 (MQ=255)
                                                                                 tACGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACtatt                    >  1:1346601/1‑100 (MQ=255)
                                                                                 tACGCCCTCCCGCAGGAGAGCGTACTCACTACTTTTACCCTCACAGTTTAAGATTTTACTGCTTTCAAGATAAATAACGTCGGTTTATAAACAGACtatt                    >  2:3593120/1‑100 (MQ=255)
                                                                                  aCGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACtatta                   >  1:610959/1‑100 (MQ=255)
                                                                                  aCGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACtatta                   <  2:2430054/100‑1 (MQ=255)
                                                                                  aCGCCCTCCCGCAGGAGAGCGTACTGACTCCTTTTACCCTCCCAGTATAAGATTTTCCTGACTTCAAGATATATAAAGTCGGTTTATACACAGACtatta                   >  2:2129593/1‑100 (MQ=255)
                                                                                  aCGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACtatta                   >  1:2720546/1‑100 (MQ=255)
                                                                                   cGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACtattat                  >  1:3648781/1‑100 (MQ=255)
                                                                                   cGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACtattat                  >  2:709335/1‑100 (MQ=255)
                                                                                       ctcCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACTGACTATTATCaca              <  1:2638570/100‑3 (MQ=255)
                                                                                       ctcCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAta              <  1:710896/100‑1 (MQ=255)
                                                                                       ctcCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAta              >  1:1237801/1‑100 (MQ=255)
                                                                                        tcCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTATCAAGATATATAACAACAAAATATAAACAAACtattataaaa              >  1:335052/1‑95 (MQ=255)
                                                                                          cTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCACGATATATAACGTCGGTTTATAAACAGACTATTATCAtaggt           >  2:1871534/1‑100 (MQ=255)
                                                                                          cTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAtaggt           >  2:104105/1‑100 (MQ=255)
                                                                                           tGCGGGAGGGCGTACTCACTCCTTTTACCCTCACTGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAtaggta          <  1:3202036/100‑1 (MQ=255)
                                                                                           tGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAtaggta          >  2:668263/1‑100 (MQ=255)
                                                                                               ggcgggCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTGTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGATAGGc       >  2:441072/4‑99 (MQ=255)
                                                                                               ggAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGc       <  2:3314178/99‑1 (MQ=255)
                                                                                               ggAGGGCGTACTCACTCCTTTTACCCGCACAGTGTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGcc      <  2:3404777/100‑1 (MQ=255)
                                                                                                   ggCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTGCAAGATATATAACGTCGGTTTATAAACAGGCTATTATCATAGGAAGGCCAAAt  >  2:3739887/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 < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: