Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,306,787 G→A 15.8% intergenic (+19/+34) 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,7870GA15.8% 225.2 / 28.2 94intergenic (+19/+34)oppF/yciUoligopeptide ABC transporter ATPase/UPF0263 family protein
Reads supporting (aligned to +/- strand):  ref base G (43/36);  new base A (9/6);  total (52/42)
Fisher's exact test for biased strand distribution p-value = 7.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.48e-01

CCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGC  >  NC_000913/1306688‑1306883
                                                                                                   |                                                                                                
ccAGAGTGCGCCAAAACACGTCCTGTTCTTGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAg                                                                                                  >  2:1240695/1‑100 (MQ=255)
  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                                           >  1:1689715/1‑100 (MQ=255)
                                                       cTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATgg                                           <  2:352023/100‑1 (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                             >  2:1653769/1‑100 (MQ=255)
                                                                     gATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATc                             <  2:2517981/100‑1 (MQ=255)
                                                                     gATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATc                             <  1:387066/100‑1 (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                   >  1:775089/1‑100 (MQ=255)
                                                                               aaaaaCAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCgg                   <  1:2506600/100‑1 (MQ=255)
                                                                               aaaaaCAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCgg                   <  2:52883/100‑1 (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)
                                                                                                   |                                                                                                
CCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGC  >  NC_000913/1306688‑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: