Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000731 81 N→G 100% intergenic (–/‑151)  / → repA –/replication protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000731810NG100.0% 5137.3 / ‑6.5 1364intergenic (–/‑151)–/repA–/replication protein A
Reads supporting (aligned to +/- strand):  ref base N (0/0);  major base G (793/567);  minor base A (2/0);  total (797/567)
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 1364 total aligned reads displayed.

‑‑‑‑‑‑‑‑ATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTANAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAA  >  CP000731/1‑170
                                                                                        |                                                                                         
ataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGaa                                                                                         >  1:2717175‑M1/9‑91 (MQ=255)
 taagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGaaa                                                                                        >  2:2054417‑M1/8‑91 (MQ=255)
 taagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGaaa                                                                                        >  1:863327‑M1/8‑91 (MQ=255)
 taagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGaaa                                                                                        >  2:2221242‑M1/8‑91 (MQ=255)
  aagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAAc                                                                                       >  1:3180261‑M1/7‑91 (MQ=255)
   agtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACt                                                                                      <  1:2792517‑M1/86‑1 (MQ=255)
    gtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTg                                                                                     <  1:2505242‑M1/87‑1 (MQ=255)
    gtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTg                                                                                     <  2:1507030‑M1/87‑1 (MQ=255)
          ttAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGtttatt                                                                               >  1:1640479/1‑91 (MQ=255)
           tAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGtttattt                                                                              >  2:2717155/1‑91 (MQ=255)
           tAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGtttattt                                                                              >  1:1994701/1‑91 (MQ=255)
            aaaGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGtttattta                                                                             <  1:1766230/91‑1 (MQ=255)
             aaGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTttatttat                                                                            <  1:3420520/91‑1 (MQ=255)
              aGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTtata                                                                           <  2:2548643/91‑1 (MQ=255)
               gTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTtatat                                                                          >  1:2231962/1‑91 (MQ=255)
               gTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTtatat                                                                          <  1:1713795/91‑1 (MQ=255)
                 gCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATaa                                                                        >  2:3727501/1‑91 (MQ=255)
                   tataATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAAtt                                                                      >  2:1601679/1‑91 (MQ=255)
                    ataataAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAAttt                                                                     >  1:930919/1‑91 (MQ=255)
                     taataaACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAAtttt                                                                    <  1:1132092/91‑1 (MQ=255)
                        taaACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGa                                                                 <  1:1955083/91‑1 (MQ=255)
                        taaACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGa                                                                 >  2:397760/1‑91 (MQ=255)
                           aCTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATAc                                                              >  1:2573558/1‑91 (MQ=255)
                            gttAATGTAGTAAATACCCGAGCTGCTGGTAATAGTCGGGGTTTGGTAATAAGTCGTTTAGAAACTGTTTATTTATATAGTTTTCGATACt                                                             <  2:915265/90‑1 (MQ=255)
                            cTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTACTAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGAtctt                                                             >  2:46705/1‑88 (MQ=255)
                            cTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACt                                                             <  2:3661755/91‑1 (MQ=255)
                            cTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACt                                                             <  1:3636518/91‑1 (MQ=255)
                              tAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACtcg                                                           >  2:3416183/1‑89 (MQ=255)
                                  ataACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGtatt                                                       <  1:267812/91‑1 (MQ=255)
                                  ataACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGtatt                                                       <  2:1835782/91‑1 (MQ=255)
                                   taACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGtatta                                                      >  1:2053414/1‑91 (MQ=255)
                                      cAAAAACCCCATCAGCTGTTAATAGTCGGTGTTTTGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATtatat                                                   <  2:3345989/91‑1 (MQ=255)
                                      cAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATtatat                                                   >  1:3496938/1‑91 (MQ=255)
                                      cAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATtatat                                                   >  1:3852503/1‑91 (MQ=255)
                                      cAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATtatat                                                   <  1:529467/91‑1 (MQ=255)
                                       aaaaaCCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATtatata                                                  <  2:2837276/91‑1 (MQ=255)
                                          aaccccgactattaccagtagttggggtttGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATg                                               >  1:3796586‑M2/31‑91 (MQ=255)
                                          AACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATg                                               <  1:1007725‑M1/91‑1 (MQ=255)
                                           ACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGa                                              <  1:184796‑M1/91‑1 (MQ=255)
                                           ACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGa                                              >  1:1391868‑M1/1‑91 (MQ=255)
                                           ACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACCTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGa                                              <  2:3251465‑M1/91‑1 (MQ=255)
                                            ccccgactattaccagtagttggggtttGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGaa                                             >  1:789465‑M2/29‑91 (MQ=255)
                                            CCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGaa                                             <  1:3311220‑M1/91‑1 (MQ=255)
                                            CCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGaa                                             <  1:2277080‑M1/91‑1 (MQ=255)
                                             CCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGaat                                            >  1:1455352‑M1/1‑91 (MQ=255)
                                             CCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGaat                                            <  2:2420410‑M1/91‑1 (MQ=255)
                                             CCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGaat                                            <  1:1230872‑M1/91‑1 (MQ=255)
                                             CCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGaat                                            >  1:628471‑M1/1‑91 (MQ=255)
                                              CCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGaata                                           <  1:3708524‑M1/91‑1 (MQ=255)
                                               CAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGaataa                                          >  2:1028952‑M1/1‑91 (MQ=255)
                                               CAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGaataa                                          >  2:1248312‑M1/1‑91 (MQ=255)
                                               CAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGaataa                                          <  2:1928612‑M1/91‑1 (MQ=255)
                                               CAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGaataa                                          >  2:2490737‑M1/1‑91 (MQ=255)
                                                aaCTACTGGTAATAGTCGGGGTTTGGTAATAACACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGaataat                                         >  2:3885048/1‑91 (MQ=255)
                                                   tattaccagtagttggggtttGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATc                                      >  1:2864361‑M2/22‑91 (MQ=255)
                                                   tattaccagtagttggggtttGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATc                                      >  2:570693‑M2/22‑91 (MQ=255)
                                                   tattaccagtagttggggtttGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATATTCGATACTTGTATTATATAATGAATAATATc                                      >  1:3365813‑M2/22‑91 (MQ=255)
                                                   tactGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATc                                      >  2:3103758/1‑91 (MQ=255)
                                                    attaccagtagttggggtttGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCt                                     >  1:1063777‑M2/21‑91 (MQ=255)
                                                    attaccagtagttggggtttGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCt                                     >  2:1331862‑M2/21‑91 (MQ=255)
                                                    attaccagtagttggggtttGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCt                                     >  1:1509980‑M2/21‑91 (MQ=255)
                                                    actGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCt                                     <  1:487300/91‑1 (MQ=255)
                                                    actGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGATTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCt                                     <  1:3943935/91‑1 (MQ=255)
                                                     ctagTAATCGACGCGGTATGGCAATATAACGTCTAGCATCAATTTATTTATATAATTTCCAATACTTCTATTATATAATTAATAATATCtt                                    <  2:115070/88‑1 (MQ=255)
                                                     ctGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCtt                                    >  2:1076734/1‑91 (MQ=255)
                                                      taccagtagttggggtttGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTa                                   >  1:154289‑M2/19‑91 (MQ=255)
                                                      tGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTa                                   <  1:1909479/91‑1 (MQ=255)
                                                      tGGAAATAGTCGGGGTTTGGTAATAAAACGTTTAGTAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATAGCTTa                                   <  2:3813563/91‑1 (MQ=255)
                                                        gTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGt                                 <  2:648623/91‑1 (MQ=255)
                                                              gTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGc                           <  1:741871/91‑1 (MQ=255)
                                                               tCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCa                          >  2:2136467/1‑91 (MQ=255)
                                                               tCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCa                          <  2:3485009/91‑1 (MQ=255)
                                                                cGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCAc                         <  1:1305046/91‑1 (MQ=255)
                                                                 ggggtttGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCAcc                        <  1:2143973‑M2/84‑1 (MQ=255)
                                                                 ggggtttGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCAcc                        >  2:3585225‑M2/8‑91 (MQ=255)
                                                                 ggggtttGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCAcc                        <  1:1367202‑M2/84‑1 (MQ=255)
                                                                 ggggtttGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCAcc                        >  2:3012060‑M2/8‑91 (MQ=255)
                                                                 ggggtttGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCAcc                        <  1:3001826‑M2/84‑1 (MQ=255)
                                                                 ggggtttGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATCTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCAcc                        <  1:3129904‑M2/84‑1 (MQ=255)
                                                                  gggtttGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCt                       >  2:2765868‑M2/7‑91 (MQ=255)
                                                                   ggtttGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTa                      >  2:199134‑M2/6‑91 (MQ=255)
                                                                   ggtttGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTa                      <  1:3570315‑M2/86‑1 (MQ=255)
                                                                    gtttGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTaa                     <  1:3235923‑M2/87‑1 (MQ=255)
                                                                        ggTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCa                 <  2:855788/91‑1 (MQ=255)
                                                                         gTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCAt                <  2:2349138/91‑1 (MQ=255)
                                                                         gTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCAt                >  1:97249/1‑91 (MQ=255)
                                                                          taataaAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCAtt               >  2:1670284/1‑91 (MQ=255)
                                                                             taaAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATtaaa            >  2:2951895/1‑91 (MQ=255)
                                                                             taaAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATtaaa            <  1:3080846/91‑1 (MQ=255)
                                                                              aaaaCGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATtaaat           >  1:2432801/1‑91 (MQ=255)
                                                                              aaaaCGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATtaaat           >  1:2351997/1‑91 (MQ=255)
                                                                               aaaCGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATtaaata          <  1:2934156/91‑1 (MQ=255)
                                                                                aaCGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATtaaataa         <  2:3148952/91‑1 (MQ=255)
                                                                                 aCGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATtaaataaa        >  1:42864/1‑91 (MQ=255)
                                                                                 aCGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATtaaataaa        >  2:745207/1‑91 (MQ=255)
                                                                                    tttAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGt     >  1:1606487/1‑91 (MQ=255)
                                                                                      gaGAAGCTGTTTATTTATATAATGGTGGATGCTGGTATTATATAGTGAAGAATATGGTAGGGCAAGCACCTAAATCAGTAAATAAATGTaa   <  2:1416596/90‑1 (MQ=255)
                                                                                       aGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTaaa  >  2:2888528/1‑91 (MQ=255)
                                                                                       aGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTaaa  >  2:2624218/1‑91 (MQ=255)
                                                                                       aGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTaaa  >  2:173823/1‑91 (MQ=255)
                                                                                        |                                                                                         
‑‑‑‑‑‑‑‑ATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTANAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAA  >  CP000731/1‑170

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