Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,306,788 T→C 17.7% intergenic (+20/+33) oppF → / ← yciU oligopeptide ABC transporter ATPase/UPF0263 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,306,7880TC17.7% 342.8 / 69.7 158intergenic (+20/+33)oppF/yciUoligopeptide ABC transporter ATPase/UPF0263 family protein
Reads supporting (aligned to +/- strand):  ref base T (61/69);  new base C (15/13);  total (76/82)
Fisher's exact test for biased strand distribution p-value = 5.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.07e-02

Only 100 of 159 total aligned reads displayed.

CAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCA  >  NC_000913/1306689‑1306884
                                                                                                   |                                                                                                
   agTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTg                                                                                                >  1:1885104/1‑99 (MQ=255)
    gTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTc                                                                                              <  1:1800679/100‑1 (MQ=255)
     tGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCa                                                                                             <  2:564287/100‑1 (MQ=255)
       cgcCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAg                                                                                            >  1:2099352/1‑99 (MQ=255)
       cgcCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAg                                                                                            >  1:133092/1‑99 (MQ=255)
       cgcCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAg                                                                                            <  2:3761614/99‑1 (MQ=255)
       cgcCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGc                                                                                           >  1:23146/1‑100 (MQ=255)
        gcCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTCAAAACAATAAGGGCTGACAGTTGTCAGcc                                                                                          >  1:2998840/1‑100 (MQ=255)
        gcCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGCCAGTTGTCAGcc                                                                                          >  1:3292247/1‑100 (MQ=255)
        gcCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGcc                                                                                          >  2:3135214/1‑100 (MQ=255)
        gcCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAACTGTCAGcc                                                                                          <  2:216182/100‑1 (MQ=255)
         ccAAAACACGTCCTGTTCTGGAGGGGCGTTGCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAACTGTCAGcc                                                                                          <  2:1799331/99‑1 (MQ=255)
         ccAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGccc                                                                                         >  2:2230335/1‑100 (MQ=255)
          cAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAACTGTCAGCCCt                                                                                        <  2:1654906/100‑1 (MQ=255)
           aaaaCACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCtt                                                                                       >  2:3204273/1‑100 (MQ=255)
            aaaCACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCttt                                                                                      <  2:2091425/100‑1 (MQ=255)
             aaCACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCtttt                                                                                     <  2:2816388/100‑1 (MQ=255)
             aaCACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCtttt                                                                                     <  1:617218/100‑1 (MQ=255)
             aaCACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAACTGTCAGCCCttat                                                                                     <  2:1283845/100‑3 (MQ=255)
             aaCACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAACTGTCAGCCCttat                                                                                     <  1:862148/100‑3 (MQ=255)
             aaCACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAACTGTCAGCCCttat                                                                                     <  1:312355/100‑3 (MQ=255)
              acacGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGCTCCGCTTTAACAACAATAAGGGCTGACAGTTGTCAGCCCtttt                                                                                     <  2:3292225/99‑1 (MQ=255)
              acacGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTCAAAAACAATAAGGGCTGACAGTTGTCAGCCCttttt                                                                                    <  1:2788496/100‑1 (MQ=255)
               cacGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTc                                                                                   >  1:972082/1‑100 (MQ=255)
               cacGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCAATCCGATTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTc                                                                                   <  2:354680/100‑1 (MQ=255)
                acGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCa                                                                                  <  2:242265/100‑1 (MQ=255)
                 cGTCCTGTTGTCGAGGGGAGTTTCCGCCACACCGTTTCTTGCTTGACACTCGATCCGCTTTAAAAACAATAAGGGCTGACAACTGTCAGCCCttattgt                                                                                  <  2:175421/99‑3 (MQ=21)
                 cGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCAc                                                                                 <  2:75457/100‑1 (MQ=255)
                 cGTCCTGTTCTGGAAGGTAGTTTCAGACACTCCGTTTCTTGCCTTAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAACTGTCAGCCCttattg                                                                                   <  2:277187/98‑2 (MQ=37)
                   tCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGc                                                                               >  1:35671/1‑100 (MQ=255)
                   tCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGc                                                                               <  1:3824716/100‑1 (MQ=255)
                   tCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGc                                                                               >  1:1474634/1‑100 (MQ=255)
                   tCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAACTGTCAGCCCttattgtttt                                                                               <  1:418850/100‑9 (MQ=38)
                    ccTGTTCTGGAGGGGAGTTTCAGACCCTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCt                                                                              <  2:2027798/100‑1 (MQ=255)
                    ccTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAACTGTCAGCCCttattgttttt                                                                              <  2:601564/100‑10 (MQ=37)
                      tGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTaa                                                                            >  2:739177/1‑100 (MQ=255)
                       gTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTaaa                                                                           <  1:1075992/100‑1 (MQ=255)
                         tCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAg                                                                         >  2:1025225/1‑100 (MQ=255)
                            ggAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGa                                                                      >  1:2347408/1‑100 (MQ=255)
                            ggAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGa                                                                      >  1:44846/1‑100 (MQ=255)
                             gAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGat                                                                     <  1:64964/100‑1 (MQ=255)
                              aGGGGAGTTTCAGACACTCCGTTTCTTGCATGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGatt                                                                    <  1:778866/100‑1 (MQ=255)
                               ggggAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGatta                                                                   <  1:2428475/100‑1 (MQ=255)
                               ggggAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGatta                                                                   >  1:2243651/1‑100 (MQ=255)
                                 ggAGTTTCAGACACACCATTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGattatt                                                                 <  2:2688497/100‑1 (MQ=255)
                                   attttCAGACACGCCGTTTATTGCCGGAAAGTCGATCCGCTTTAAAAACTATAAGGGCTGACAGTTGGCAGGCCGTTTTCACGCTAAAAGCGattatt                                                                 <  2:3583251/96‑1 (MQ=255)
                                    gTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGAttatttat                                                              <  2:611955/100‑1 (MQ=255)
                                      ttCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTc                                                            >  2:919853/1‑100 (MQ=255)
                                        cAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTccc                                                          <  1:1821602/100‑1 (MQ=255)
                                        cAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTGTTCACGCTAAAAGCGATTATTTATTccc                                                          >  2:133130/1‑100 (MQ=255)
                                         aGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCg                                                         <  2:30544/100‑1 (MQ=255)
                                           acacTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGcc                                                       >  1:1364912/1‑100 (MQ=255)
                                            cacTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCa                                                      <  1:985402/100‑1 (MQ=255)
                                            cacTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCa                                                      <  2:1652741/100‑1 (MQ=255)
                                              cTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAga                                                    >  1:1188198/1‑100 (MQ=255)
                                              cTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAga                                                    >  1:3372584/1‑100 (MQ=255)
                                                ccGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAgatg                                                  >  2:2902340/1‑100 (MQ=255)
                                                 cGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAgatga                                                 >  1:1100946/1‑100 (MQ=255)
                                                     tCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATg                                             <  2:3060677/100‑1 (MQ=255)
                                                     tCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATg                                             >  2:442159/1‑100 (MQ=255)
                                                     tCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATg                                             <  1:3128751/100‑1 (MQ=255)
                                                      cTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATgg                                            >  2:2557995/1‑100 (MQ=255)
                                                        tGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCa                                          >  1:2400025/1‑100 (MQ=255)
                                                         gCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCaa                                         >  2:2500273/1‑100 (MQ=255)
                                                          ccTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCaaa                                        >  1:3859369/1‑100 (MQ=255)
                                                           cTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAg                                       <  2:3337407/100‑1 (MQ=255)
                                                               aaGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTtg                                   <  1:1577826/100‑1 (MQ=255)
                                                                 gTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTtgtg                                 <  1:3204295/100‑1 (MQ=255)
                                                                      tCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCtt                            <  1:2202783/100‑1 (MQ=255)
                                                                            ttAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTcgc                      <  1:3631968/100‑1 (MQ=255)
                                                                            ttAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTcgc                      >  1:3917354/1‑100 (MQ=255)
                                                                             tAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTcgcg                     <  2:3759594/100‑1 (MQ=255)
                                                                              aaaaaCAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCgg                    <  2:3912417/100‑1 (MQ=255)
                                                                               aaaaCAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGc                   <  1:2048529/100‑1 (MQ=255)
                                                                               aaaaCAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGc                   >  1:53065/1‑100 (MQ=255)
                                                                                 aaCAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCAt                 <  2:2787158/100‑1 (MQ=255)
                                                                                    aaTAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAAc              <  1:1583337/100‑1 (MQ=255)
                                                                                     aaaaGGGCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACa             <  2:3944719/98‑1 (MQ=255)
                                                                                     aaaaGGGCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACa             >  1:2751593/3‑100 (MQ=255)
                                                                                     aTAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACa             >  1:760286/1‑100 (MQ=255)
                                                                                     aTAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACa             >  1:1773806/1‑100 (MQ=255)
                                                                                      tAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAg            <  2:1456817/100‑1 (MQ=255)
                                                                                      aaaGGGCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAg            >  2:1002936/2‑100 (MQ=255)
                                                                                      aaaGGGCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAg            >  2:3945159/2‑100 (MQ=255)
                                                                                       aaGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAGAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTCTGTGATCTTTTTCGCGGCATAACAGa           >  2:574107/1‑100 (MQ=255)
                                                                                       aaGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGa           <  1:3029856/100‑1 (MQ=255)
                                                                                       aaGGGCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGa           >  1:187683/1‑100 (MQ=255)
                                                                                       aaGGGCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGa           >  1:3543411/1‑100 (MQ=255)
                                                                                        aGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGaa          >  1:2641134/1‑100 (MQ=255)
                                                                                        aGGGCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGaa          >  1:2150318/1‑100 (MQ=255)
                                                                                        aGGGCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGaa          >  1:1375919/1‑100 (MQ=255)
                                                                                        aGGGCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGaa          >  2:2959733/1‑100 (MQ=255)
                                                                                        aGGGCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGaa          >  1:3920375/1‑100 (MQ=255)
                                                                                             tGACAGTTGTCAGCCCTTTTTTTCTTTTTTTTCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCgg     >  2:1830574/1‑100 (MQ=255)
                                                                                             tGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCgg     >  1:2462391/1‑100 (MQ=255)
                                                                                              gACAGTTGTCAGCTCTTTTTCACGCTGAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCggg    <  1:2702881/100‑1 (MQ=255)
                                                                                              gACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCggg    <  2:3705454/100‑1 (MQ=255)
                                                                                              gACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCggg    <  1:2698625/100‑1 (MQ=255)
                                                                                               aCAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGc   <  2:2400013/100‑1 (MQ=255)
                                                                                                 aGTTGTCAGACCTTTTTCACGCTAAAAGCGATTATTTATTCACGCCAGATGATATGGCAACGATAGTGATCTTTTTCGCGGCATAACAGAATGCGGGCa  >  2:2134488/1‑99 (MQ=255)
                                                                                                   |                                                                                                
CAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCA  >  NC_000913/1306689‑1306884

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