Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,306,788 T→C 16.5% 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,7880TC16.5% 195.4 / 31.6 93intergenic (+20/+33)oppF/yciUoligopeptide ABC transporter ATPase/UPF0263 family protein
Reads supporting (aligned to +/- strand):  ref base T (42/35);  new base C (9/6);  total (52/41)
Fisher's exact test for biased strand distribution p-value = 7.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

AGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGC  >  NC_000913/1306690‑1306883
                                                                                                  |                                                                                               
agagTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGtt                                                                                                <  1:1201330/100‑1 (MQ=255)
 gagTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTg                                                                                               <  2:5606/100‑1 (MQ=255)
      cgcCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACAGTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGc                                                                                          >  1:844026/1‑100 (MQ=255)
      cgcCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGc                                                                                          >  1:547379/1‑100 (MQ=255)
      cgcCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGc                                                                                          >  1:1650546/1‑100 (MQ=255)
      cgcCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGc                                                                                          >  2:2402256/1‑100 (MQ=255)
       gcCAAAGCACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAACTGTCAGcc                                                                                         <  2:2367011/100‑1 (MQ=255)
       gcCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGcc                                                                                         >  1:259551/1‑100 (MQ=255)
       gcCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGcc                                                                                         >  1:1010674/1‑100 (MQ=255)
        ccAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGccc                                                                                        >  2:1186794/1‑100 (MQ=255)
        ccAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGccc                                                                                        >  1:1244703/1‑100 (MQ=255)
         cAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAACTGTCAGCCCt                                                                                       <  2:2131134/100‑1 (MQ=255)
         cAAAACACGTCCTGTTCTGGAGGGGAGTTCCAGACACAACGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCt                                                                                       <  2:1236926/100‑1 (MQ=255)
          aaaaCACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCt                                                                                       >  2:1211335/1‑99 (MQ=255)
           aaaCACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAACTGTCAGCCCtta                                                                                     <  1:847267/100‑2 (MQ=255)
            aaCACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCttt                                                                                     <  1:2385587/99‑1 (MQ=255)
              ctcGTGCTGTTCCGGAGGGGAGTTTCATACACGCCTTTTCTGGCCTGAAAGTCGCTCCGCTTTAAAAACAATAAGGGCTGCCAACTGTCAGCCCttat                                                                                    <  2:2500451/96‑3 (MQ=255)
              cacTTCCTGTTCTGGATGGGATTTTAACACCCTCCGTTTCTTGCCTGAAAGTCGCTCCGCTTTAAAAACAATAAGGGCTTACAACTGTCAGCCCt                                                                                       <  2:353022/95‑1 (MQ=255)
              cacGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTc                                                                                  <  1:2427931/100‑1 (MQ=255)
               acGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAACTGTCAGCCCttattg                                                                                  <  2:928648/99‑2 (MQ=38)
                cGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCAc                                                                                >  2:1960304/1‑100 (MQ=255)
                 gTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCAc                                                                                >  2:830433/1‑99 (MQ=255)
                 gTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACg                                                                               <  2:1354524/100‑1 (MQ=255)
                   ccTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAGCAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCt                                                                             <  1:2444151/100‑1 (MQ=255)
                   ccTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGc                                                                              >  2:1173432/1‑99 (MQ=255)
                    cTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTa                                                                            >  2:796361/1‑100 (MQ=255)
                      gTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAGAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTaa                                                                           <  2:187378/99‑1 (MQ=255)
                         cTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGc                                                                       <  1:421749/100‑1 (MQ=255)
                           ggAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGa                                                                     >  2:1544611/1‑100 (MQ=255)
                           ggAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGa                                                                     >  2:387062/1‑100 (MQ=255)
                           ggAGGGCAGTTTCAGACACTCCGTTTCTTGCCTGAAATACCATCGGTTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGc                                                                       <  2:935280/98‑1 (MQ=255)
                             aGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGatt                                                                   >  1:1075072/1‑100 (MQ=255)
                             aGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGAACCGCTTTAAAAACAATAAGGGATGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGatt                                                                   >  2:336804/1‑100 (MQ=255)
                              ggggAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGatta                                                                  >  2:2545074/1‑100 (MQ=255)
                              ggggAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGatta                                                                  >  2:970327/1‑100 (MQ=255)
                                ggAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGattatt                                                                >  1:1078004/1‑100 (MQ=255)
                                ggAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGatta                                                                  >  2:959198/1‑98 (MQ=255)
                                 gAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGAttattt                                                               <  2:9546/100‑1 (MQ=255)
                                  aGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGAttattta                                                              <  1:988789/100‑1 (MQ=255)
                                   gTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGAttattta                                                              <  2:827833/99‑1 (MQ=255)
                                     ttCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGTCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAACGATTATTTATTc                                                           >  2:1174378/1‑100 (MQ=255)
                                     ttCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTc                                                           >  1:1764073/1‑100 (MQ=255)
                                      tCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTcc                                                          >  1:2069602/1‑100 (MQ=255)
                                       cAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTccc                                                         >  2:1203870/1‑100 (MQ=255)
                                          acacTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGcc                                                      <  1:1173395/100‑1 (MQ=255)
                                            acTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGcc                                                      <  2:2157632/98‑1 (MQ=255)
                                               ccGTTCCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGCGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCAAga                                                   >  2:2233311/1‑98 (MQ=255)
                                                cGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAgatga                                                <  2:652743/100‑1 (MQ=255)
                                                   ttCTTGCCTGAATGTCGATCCGCTTTAAAAACAATAATGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGatat                                             <  2:420227/100‑1 (MQ=255)
                                                   ttCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGatat                                             >  2:1328902/1‑100 (MQ=255)
                                                    tCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATg                                            <  1:1186974/100‑1 (MQ=255)
                                                     cTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATgg                                           <  2:352023/100‑1 (MQ=255)
                                                     cTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATgg                                           >  1:1689715/1‑100 (MQ=255)
                                                      ttGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGc                                          >  2:102892/1‑100 (MQ=255)
                                                        gCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCa                                         <  1:129505/99‑1 (MQ=255)
                                                        gCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGATAAAAGCGAATATTTATTCCCGCTAGATGATATGGaaa                                        >  2:2316598/1‑100 (MQ=255)
                                                         ccTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCaaa                                       >  2:514402/1‑100 (MQ=255)
                                                         ccTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCaaa                                       >  1:2073650/1‑100 (MQ=255)
                                                          cTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAg                                      <  2:1515514/100‑1 (MQ=255)
                                                            gAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGtt                                    >  2:2262507/1‑100 (MQ=255)
                                                              aaGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTtg                                  <  2:2485541/100‑1 (MQ=255)
                                                                 tCGATCCGCTTTAAAAACAATAAGGACTGACAGTTGTCAGCCCTTTTCGACGCTAAAAGCGATCATTTATTCCCGCAATATGATATGGCAATGttttt                                 >  2:646454/1‑96 (MQ=255)
                                                                  cGATCCGCTTTAAAAACAATAAGGGCTGACAGATGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGAt                              >  1:2426211/1‑100 (MQ=255)
                                                                   gATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATc                             <  1:387066/100‑1 (MQ=255)
                                                                   gATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATc                             <  2:2517981/100‑1 (MQ=255)
                                                                   gATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATc                             >  2:1653769/1‑100 (MQ=255)
                                                                       cGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCtttt                         <  1:2043902/100‑1 (MQ=255)
                                                                         cTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTc                       >  2:1276498/1‑100 (MQ=255)
                                                                          tttAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTcg                      <  2:1760153/100‑1 (MQ=255)
                                                                            tAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCCGATGATATGGCAAAGTTTGTGATCTTTTTcgcg                    <  1:640250/100‑1 (MQ=255)
                                                                            tAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTcgcg                    >  1:1303252/1‑100 (MQ=255)
                                                                             aaaaaCAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCgg                   <  2:52883/100‑1 (MQ=255)
                                                                             aaaaaCAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCgg                   <  1:2506600/100‑1 (MQ=255)
                                                                             aaaaaCAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCgg                   >  1:775089/1‑100 (MQ=255)
                                                                              aaaaCAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGTGGc                  <  1:1207643/100‑1 (MQ=255)
                                                                              aaaaCAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGc                  <  1:227337/100‑1 (MQ=255)
                                                                               cgtgaaaaaGGGCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCa                 >  1:1346868/8‑100 (MQ=38)
                                                                               aaaCAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCa                 <  1:1545790/100‑1 (MQ=255)
                                                                               aaaCAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCa                 >  1:784592/1‑100 (MQ=255)
                                                                                gtgaaaaaGGGCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCAt                >  1:1979127/7‑100 (MQ=38)
                                                                                gtgaaaaaGGGCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCAt                >  1:1979008/7‑100 (MQ=38)
                                                                                 tgaaaaaGGGCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGCTATGGCAACGTTTGGGATCTTTTTCGCGGCATa               >  1:1159977/6‑100 (MQ=37)
                                                                                  cAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATaa              >  2:2253557/1‑100 (MQ=255)
                                                                                    aaaaGGGCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACa            >  1:903847/3‑100 (MQ=255)
                                                                                       aGGGCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGaa         >  1:1244144/1‑100 (MQ=255)
                                                                                       aGGGCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGa          >  1:1201591/1‑99 (MQ=255)
                                                                                        gggCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAAt        >  2:1112562/1‑100 (MQ=255)
                                                                                         ggCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATg       >  1:2473484/1‑100 (MQ=255)
                                                                                          gCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGc      >  1:677096/1‑100 (MQ=255)
                                                                                            tGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCgg    <  1:1265364/100‑1 (MQ=255)
                                                                                            tGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCgg    >  1:2269283/1‑100 (MQ=255)
                                                                                             gACAGTTGTCAGCCCTTTTTCACGCTAAATGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCggg   <  1:796368/100‑1 (MQ=255)
                                                                                              aCAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGc  <  1:1703461/100‑1 (MQ=255)
                                                                                              aCAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATAGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCgg    <  2:239354/98‑1 (MQ=255)
                                                                                                  |                                                                                               
AGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGC  >  NC_000913/1306690‑1306883

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: