Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,663,482 A→C 6.5% K14T (AAG→ACG)  suhB → inositol monophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,663,4820AC 6.5% 329.3 / 12.3 92K14T (AAG→ACG) suhBinositol monophosphatase
Reads supporting (aligned to +/- strand):  ref base A (44/42);  new base C (3/3);  total (47/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01

TCTAAAGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTAAATCTTACCCACAGCAC  >  NC_000913/2663336‑2663627
                                                                                                                                                  |                                                                                                                                                 
tCTAAAGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAgg                                                                                                                                                 >  4:1327026/1‑149 (MQ=255)
          attCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCg                                                                                                                                               <  3:1257268/141‑1 (MQ=255)
          attCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCg                                                                                                                                               >  4:1257268/1‑141 (MQ=255)
            tCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGtaatttaa                                                                                                                                     <  3:215009/149‑1 (MQ=255)
              aCGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCTTGCGCGCTGCGCGCACgg                                                                                                                                                 <  1:152355/135‑1 (MQ=255)
              aCGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCTTGCGCGCTGCGCGCACgg                                                                                                                                                 >  2:152355/1‑135 (MQ=255)
               cGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTg                                                                                                                                  <  3:1327026/149‑1 (MQ=255)
                 cATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGtaattt                                                                                                                                       <  1:1203896/142‑1 (MQ=255)
                 cATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGtaattt                                                                                                                                       >  2:1203896/1‑142 (MQ=255)
                    cTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAgg                                                                                                                                                 <  4:320633/129‑1 (MQ=255)
                    cTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAgg                                                                                                                                                 >  3:320633/1‑129 (MQ=255)
                     ttATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCgg                                                                                                                                              <  2:655249/131‑1 (MQ=255)
                     ttATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCgg                                                                                                                                              >  1:655249/1‑131 (MQ=255)
                          aTAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTAt                                                                                                                       <  1:989740/149‑1 (MQ=255)
                          aTAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTAt                                                                                                                       >  2:989740/1‑149 (MQ=255)
                          aTAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTAt                                                                                                                       >  3:1288592/1‑149 (MQ=255)
                          aTAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTAt                                                                                                                       >  1:883394/1‑149 (MQ=255)
                            aaaaCGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCATTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGa                                                                                                                     <  2:883394/149‑1 (MQ=255)
                               aCGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGtaatttaatt                                                                                                                                   <  3:131628/132‑1 (MQ=255)
                               aCGAAGACAGAAGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAAAATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCTGTGCGCGCAGCGCGCAAGGCGGGtaatttaatt                                                                                                                                   >  4:131628/1‑132 (MQ=255)
                                           gCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCaaaaa                                                                                                                           >  3:3906/1‑128 (MQ=255)
                                           gCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCaaaaa                                                                                                                           <  4:3906/128‑1 (MQ=255)
                                              gATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGtaatttaat                                                                                                                                    <  2:1332117/116‑1 (MQ=255)
                                              gATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGtaatttaat                                                                                                                                    >  1:1332117/1‑116 (MQ=255)
                                              gATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGa                                                                                                   >  3:910733/1‑149 (MQ=255)
                                                   cGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGa                                                                                              <  4:910733/149‑1 (MQ=255)
                                                     ctgctATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTa                                                                                                                        >  4:1165752/1‑121 (MQ=255)
                                                     ctgctATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTa                                                                                                                        <  3:1165752/121‑1 (MQ=255)
                                                        ctATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGc                                                                                                                                 >  2:370305/1‑109 (MQ=255)
                                                        ctATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGc                                                                                                                                 <  1:370305/109‑1 (MQ=255)
                                                        ctATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAg                                                                                         >  2:102326/1‑149 (MQ=255)
                                                         tataCTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGa                                                                                        >  3:1401291/1‑149 (MQ=255)
                                                           taCTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGtaa                                                                                                                                          >  3:120405/1‑97 (MQ=255)
                                                           taCTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGtaa                                                                                                                                          <  4:120405/97‑1 (MQ=255)
                                                             ctctGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTAt                                                                                                                       >  4:677270/1‑114 (MQ=255)
                                                             ctctGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTAt                                                                                                                       <  3:677270/114‑1 (MQ=255)
                                                             ctctGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAgg                                                                                    <  4:1288592/149‑1 (MQ=255)
                                                             ctctGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAgg                                                                                    <  4:1401291/149‑1 (MQ=255)
                                                             ctctGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAgg                                                                                    >  2:841334/1‑149 (MQ=255)
                                                              tctGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGaa                                                                                       >  3:590153/1‑145 (MQ=255)
                                                              tctGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGaa                                                                                       <  4:590153/145‑1 (MQ=255)
                                                              tctGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGc                                                                                   <  3:1204615/149‑1 (MQ=255)
                                                              tctGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGc                                                                                   >  4:1204615/1‑149 (MQ=255)
                                                               ctGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCa                                                                                  <  1:102326/149‑1 (MQ=255)
                                                                tGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGcc                                                                                                                                >  2:196070/1‑102 (MQ=255)
                                                                tGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGcc                                                                                                                                <  1:196070/102‑1 (MQ=255)
                                                                  cgcCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGc                                                                                                                                 >  2:441764/1‑99 (MQ=255)
                                                                  cgcCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGc                                                                                                                                 <  1:441764/99‑1 (MQ=255)
                                                                       ttttCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGtaatt                                                                                                                                        <  4:464957/87‑1 (MQ=255)
                                                                       ttttCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGtaatt                                                                                                                                        >  3:464957/1‑87 (MQ=255)
                                                                           cccGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAAccc                                                                                                                >  4:561599/1‑107 (MQ=255)
                                                                           cccGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAAccc                                                                                                                <  3:561599/107‑1 (MQ=255)
                                                                               ttCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGAt                                                                          <  2:235545/141‑1 (MQ=255)
                                                                               ttCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGAt                                                                          >  1:235545/1‑141 (MQ=255)
                                                                                  tttAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGa                                                                   >  2:600233/1‑145 (MQ=255)
                                                                                  tttAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGa                                                                   <  1:600233/145‑1 (MQ=255)
                                                                                    tAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGa                                                                                              <  2:714900/116‑1 (MQ=255)
                                                                                    tAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGa                                                                                              >  1:714900/1‑116 (MQ=255)
                                                                                     aaCATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCa                                                                                  >  3:606045/1‑127 (MQ=255)
                                                                                     aaCATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCa                                                                                  <  4:606045/127‑1 (MQ=255)
                                                                                     aaCATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGa                                                                           >  3:388990/1‑134 (MQ=255)
                                                                                     aaCATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGa                                                                           <  4:388990/134‑1 (MQ=255)
                                                                                      aCATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTa                                                           <  1:841334/149‑1 (MQ=255)
                                                                                        aTCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGtaatttaatt                                                                                                                                   <  4:319077/75‑1 (MQ=255)
                                                                                        aTCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGtaatttaatt                                                                                                                                   >  3:319077/1‑75 (MQ=255)
                                                                                                          aTGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGaaa                                                                                      <  4:1328631/102‑1 (MQ=255)
                                                                                                          aTGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGaaa                                                                                      >  3:1328631/1‑102 (MQ=255)
                                                                                                             cATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGa                                                                                              >  4:1205788/1‑91 (MQ=255)
                                                                                                             cATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGCAGAAGCGa                                                                                              <  3:1205788/91‑1 (MQ=255)
                                                                                                               tCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAg                                                                                             <  4:723438/90‑1 (MQ=255)
                                                                                                               tCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAg                                                                                             >  3:723438/1‑90 (MQ=255)
                                                                                                               tCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTa                                  >  4:1460946/1‑149 (MQ=255)
                                                                                                                ccGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTAt                                 <  3:1460946/149‑1 (MQ=255)
                                                                                                                 cGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAg                                                                                             <  1:1319823/88‑1 (MQ=255)
                                                                                                                 cGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAg                                                                                             >  2:1319823/1‑88 (MQ=255)
                                                                                                                  gATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCg                               >  4:663100/1‑149 (MQ=255)
                                                                                                                       tGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGt                                       >  1:697563/1‑136 (MQ=255)
                                                                                                                       tGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGt                                       <  2:697563/136‑1 (MQ=255)
                                                                                                                               gccttgcgcgctgcgcgcACGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCa                                                                <  3:502228/99‑1 (MQ=255)
                                                                                                                               gCCTTGCGCGCTGCGCGCACGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCa                                                                <  4:371395/103‑1 (MQ=255)
                                                                                                                               gCCTTGCGCGCTGCGCGCACGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCa                                                                >  3:371395/1‑103 (MQ=255)
                                                                                                                               gCCTTGCGCGCTGCGCGCACGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCa                                                                >  4:502228/1‑103 (MQ=255)
                                                                                                                                ccGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTaaa                 <  3:663100/149‑1 (MQ=255)
                                                                                                                                       cgcagcgcgcaAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGt                                                                     >  4:1315616/1‑90 (MQ=255)
                                                                                                                                       cgcagcgcgcaAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGt                                                                     <  3:1315616/90‑1 (MQ=255)
                                                                                                                                       cgcagcgcgcaAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTAAATCTTAcc          >  3:574802/1‑149 (MQ=255)
                                                                                                                                          agcgcgcaAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTAAATCTTa            >  1:489756/1‑144 (MQ=255)
                                                                                                                                          agcgcgcaAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTAAATCTTa            <  2:489756/144‑1 (MQ=255)
                                                                                                                                          agcgcgcaAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGACCCAGAACGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGaca                            <  4:1178161/128‑1 (MQ=255)
                                                                                                                                          agcgcgcaAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCCGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGaca                            >  3:1178161/1‑128 (MQ=255)
                                                                                                                                               gcaAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTTGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTAAATCTTACCCACAGcac  >  2:690767/1‑149 (MQ=255)
                                                                                                                                               gcaAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTAAATCTTACCCACAGcac  >  4:954972/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TCTAAAGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTAAATCTTACCCACAGCAC  >  NC_000913/2663336‑2663627

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: