Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,712,128 C→G 8.1% A213P (GCG→CCG)  trmN ← tRNA1(Val) (adenine(37)‑N6)‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,712,1280CG 8.1% 300.3 / 13.2 97A213P (GCG→CCG) trmNtRNA1(Val) (adenine(37)‑N6)‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (46/43);  new base G (5/3);  total (51/46)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.79e-01

GAACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACAT‑AAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAA  >  NC_000913/2711980‑2712273
                                                                                                                                                     |                                                                                                                                                 
gAACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACAT‑AAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGcc                                                                                                                                                  <  1:1601591/150‑1 (MQ=255)
         ccccTTTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACAT‑AAACAGATAAAAAGCCTGGGTCTGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGgag                                                                                                                                         <  2:1170479/149‑1 (MQ=255)
         aCCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACAT‑AAACTGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGgag                                                                                                                                         <  1:339936/150‑1 (MQ=255)
         aCCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACAT‑AAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGgag                                                                                                                                         >  1:266605/1‑150 (MQ=255)
          cccATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACAT‑AAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGTATAGTTCTGGTCTGGCCCACGAATGACTAAGCGCTCGCTAACGCATTCTCCCGCCGGCGGGgaga                                                                                                                                        >  2:643995/1‑150 (MQ=255)
           ccATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACAT‑AAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGaa                                                                                                                                       >  2:495998/1‑150 (MQ=255)
           ccATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACAT‑AAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGaa                                                                                                                                       <  1:582219/150‑1 (MQ=255)
            cATTCCGGTCGGTCGATCCTTTCCCGCGGCAATCATTACAT‑AAACAGATAAAATTCCTGGGTCAGCGCCGTATTCGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAATGCATTCTCCCGCCTGCGGGGAGAAt                                                                                                                                      <  2:833454/150‑1 (MQ=255)
            cATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACAT‑AAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAAt                                                                                                                                      <  1:1580438/150‑1 (MQ=255)
            cATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACAT‑AAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCACACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAAt                                                                                                                                      >  2:1332765/1‑150 (MQ=255)
            cATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACAT‑AAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGGCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATGCTCCCGCCTGCGGGGAGAAt                                                                                                                                      <  1:1032832/150‑1 (MQ=255)
             aTTCCTGTCCGTCGATCCTTTCCCCCGGCAATCATTACAT‑AAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCATCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATg                                                                                                                                     <  2:1640356/150‑1 (MQ=255)
             aTTCCGGTCCGTCGATCCTTACCCCCGGCAATCATTACAT‑AAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATg                                                                                                                                     <  2:393237/150‑1 (MQ=255)
              ttCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAAACAGATAAAAAGCCTGGGGCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCCGCAGGCGGGAGAATg                                                                                                                                     <  1:1196655/150‑1 (MQ=255)
               tCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACAT‑AAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGCAATAGTGCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCCGCAGGCGGGAGAATGct                                                                                                                                   <  1:1111880/150‑2 (MQ=255)
                ccGGTCCGGCGATCCTTTCCCCCGGCAATCATTACCT‑AAACAGATAAAACGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCca                                                                                                                                  <  1:655881/150‑1 (MQ=255)
                 cGGTCCGTCGCTCCTTTCCCCCGGCAATCATTACAT‑AAACAGATAAACAGCCTTGGGCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCCGCAGGCGGGAGAATGcttt                                                                                                                                 <  1:645678/150‑4 (MQ=255)
                 cGGTCCGTCGATCCTTTCCCCCGGCAATCATTACAT‑AAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCcag                                                                                                                                 <  2:642944/150‑1 (MQ=255)
                  ggTCCGTCGATCCTTTCCCCCGGCAATCATTACAT‑AAACAGATAAAAAGCCTGTGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCcagc                                                                                                                                <  2:471380/150‑1 (MQ=255)
                  ggTCCGTCGATCCTTTCCCCCGGCAATCATTACAT‑AAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCcagc                                                                                                                                <  1:965446/150‑1 (MQ=255)
                  ggTCCGTCGATCCTTTCCCCCGGCAATCATTACAT‑AAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCcagc                                                                                                                                >  2:1038085/1‑150 (MQ=255)
                   gTCCGTCGATCCTTTCCCCCGGCAATCATTACAT‑AAACAGATAAAAAGCCTGGGTCTGCGCCGTATACGCTACGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCcagca                                                                                                                               <  2:1549721/150‑1 (MQ=255)
                   gTCCGTCGATCCTTTCCCCCGGCAATCATTACAT‑AAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCcagca                                                                                                                               <  1:1011735/150‑1 (MQ=255)
                      cGTCGATCCTTTCCCCCGGCAATCATTACAT‑AAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGAGGGGAGAATTCCagcagcc                                                                                                                            >  2:1532572/1‑149 (MQ=255)
                                 tCCCCCGGCAATCATTACAT‑AAACAGATAAAAAGCCTGTGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGc                                                                                                                 <  2:1545337/150‑1 (MQ=255)
                                   ccccGGCAATCATTACAT‑AAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGc                                                                                                                    <  2:1156318/145‑1 (MQ=255)
                                   ccccGGCAATCATTACAT‑AAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGc                                                                                                                    >  1:1156318/1‑145 (MQ=255)
                                    cccGGCAATCATTACAT‑AAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGt                                                                                                              >  2:794653/1‑150 (MQ=255)
                                    cccGGCAATCATTACAT‑AAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGt                                                                                                              >  2:117502/1‑150 (MQ=255)
                                    cccGGCAATCATTACAT‑AAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGt                                                                                                              >  2:787875/1‑150 (MQ=255)
                                        gCAATCATTACAT‑AAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCTGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTcgcg                                                                                                          <  1:434973/150‑1 (MQ=255)
                                           aTCATTACAT‑AAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCtt                                                                                                       >  2:1596682/1‑150 (MQ=255)
                                                 aCAT‑AAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGtttt                                                                                                 <  1:1719/150‑1 (MQ=255)
                                                 aCAT‑AAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGtttt                                                                                                 >  2:1330650/1‑150 (MQ=255)
                                                  cAT‑AAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGTGAATGCCAGCAGCACCCGATCCGGCAGACGCGCTTCGTTTTc                                                                                                >  2:232625/1‑150 (MQ=255)
                                                           gATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACAt                                                                                        >  1:1572739/1‑150 (MQ=255)
                                                            aTAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACAtc                                                                                       <  2:949319/150‑1 (MQ=255)
                                                            aTAAAAAGCCTGGGTCAGCGCCGAATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAACCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACAtc                                                                                       <  1:155061/150‑1 (MQ=255)
                                                              aaaaaGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGTACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCcaccaccg                                                                                     >  2:118279/1‑145 (MQ=255)
                                                              aaaaaGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGAAGCACCAGATGCGGCATTCGCGCTTCGGTTTCCGCCCCAtccg                                                                                     >  1:129126/1‑150 (MQ=255)
                                                              aaaaaGCCTGGGTCAGCGCCGTATAAGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGAAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACAtccg                                                                                     <  1:935019/150‑1 (MQ=255)
                                                                     cTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCaa                                                                              >  1:1541056/1‑150 (MQ=255)
                                                                      tGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCTGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCaaa                                                                             >  1:223909/1‑150 (MQ=255)
                                                                      tGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCaaa                                                                             >  1:422437/1‑150 (MQ=255)
                                                                       gggTCAGCGCCGTATACGCTTCGGAATAGTTATGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCaacc                                                                            >  1:1184507/1‑148 (MQ=255)
                                                                        ggTCAGCGCCGTATACGCTTCGGAATAGTTCTGGGCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACg                                                                           <  2:851872/150‑1 (MQ=255)
                                                                          gaaGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGTATGACTAAGCGATGGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATTCGGCAGTCGCGCTGCGTTTTCCGCCACATCCGTCCGCAAACGTa                                                                         <  2:876029/148‑1 (MQ=255)
                                                                           cAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTaa                                                                        >  1:1503100/1‑150 (MQ=255)
                                                                           cAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTaa                                                                        >  2:1087752/1‑150 (MQ=255)
                                                                            aGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTaaa                                                                       <  1:787875/150‑1 (MQ=255)
                                                                            aGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTaaa                                                                       <  1:794653/150‑1 (MQ=255)
                                                                            aGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTaaa                                                                       <  1:117502/150‑1 (MQ=255)
                                                                              cgcCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATg                                                                     >  2:739988/1‑150 (MQ=255)
                                                                                ccGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGcc                                                                   >  2:446183/1‑150 (MQ=255)
                                                                                  gTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAg                                                                 >  1:559961/1‑150 (MQ=255)
                                                                                    ataCGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGcc                                                               <  2:1424607/150‑1 (MQ=255)
                                                                                       cGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCAt                                                            <  1:979816/150‑1 (MQ=255)
                                                                                         cTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGc                                                          >  2:1357453/1‑150 (MQ=255)
                                                                                             ggAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACAGCCGTCCGCAAACGTAAATGCCAGCCCATGCTTaa                                                      >  1:1283138/1‑150 (MQ=255)
                                                                                             ggAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCACCCGCCTGCGGGGAGAATGCCAGCAGCACACGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCCAACGTAAATGCCAGCACATGATTaa                                                      >  1:627628/1‑150 (MQ=255)
                                                                                                  aGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCa                                                 <  1:923075/150‑1 (MQ=255)
                                                                                                  aGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGAGGCAGTAGCGCTTCGTTTTCCGCCACTTCCGTCCGCAAACGTAAATGCCAGCCAATGCTTAATGCCa                                                 >  1:1069659/1‑150 (MQ=255)
                                                                                                       tggtctggCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTcc                                            <  1:890459/150‑1 (MQ=255)
                                                                                                       tggtctggCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTcc                                            <  1:1087752/150‑1 (MQ=255)
                                                                                                        ggtctggCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCCCCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCg                                           >  1:1645853/1‑150 (MQ=255)
                                                                                                        ggtctggCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCCGCCCATGCTTAATGCCAGCTCCg                                           >  2:1385384/1‑150 (MQ=255)
                                                                                                          tctggCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGTACCCTTTTCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTa                                         <  1:118279/150‑1 (MQ=255)
                                                                                                           ctggCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTaa                                        >  2:1733689/1‑150 (MQ=255)
                                                                                                              gCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCGGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAAc                                     <  1:584763/150‑1 (MQ=255)
                                                                                                               cccACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCAATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAAcc                                    >  2:1624338/1‑150 (MQ=255)
                                                                                                                     aaTGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCAGCCGATGCGGGAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTAcc                              <  2:1416797/150‑1 (MQ=255)
                                                                                                                     aaTGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCTGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTAcc                              <  1:299910/150‑1 (MQ=255)
                                                                                                                          cTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCt                         >  1:1182280/1‑150 (MQ=255)
                                                                                                                          cTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCAAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCt                         >  2:232336/1‑150 (MQ=255)
                                                                                                                           tAAGCGCTCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTg                        <  2:1170509/150‑1 (MQ=255)
                                                                                                                            aaGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGc                       >  2:654178/1‑150 (MQ=255)
                                                                                                                                 aTCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTgg                  <  1:186901/150‑1 (MQ=255)
                                                                                                                                   cGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGCCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCCGCTCCGTAAACCCATTACCAATCTGCTCTGgca                >  1:988661/1‑150 (MQ=255)
                                                                                                                                    gCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGgcag               >  2:398785/1‑150 (MQ=255)
                                                                                                                                    ctgctggCATTCTCCCCGCAGGCGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGgcag               >  2:609735/7‑150 (MQ=255)
                                                                                                                                     cTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGgcagc              <  1:1394667/150‑1 (MQ=255)
                                                                                                                                     cTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGgcagc              >  2:1728130/1‑150 (MQ=255)
                                                                                                                                      tAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCAATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGgcagca             >  2:965561/1‑150 (MQ=255)
                                                                                                                                      tAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCAAATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTTCAACCATTGCCAATCTGCTCTGgcagca             >  2:1525750/1‑150 (MQ=255)
                                                                                                                                       aaaGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCTATGCGGCAGTCGCGCTTCGTTTTCAGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGATTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGgcagcaa            >  1:507645/1‑149 (MQ=255)
                                                                                                                                       aaaGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTCCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCAc            <  1:495998/150‑1 (MQ=255)
                                                                                                                                        tggCATTCTCCCCGCAGGCGGGAGAATGCTAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACt           >  2:607438/3‑150 (MQ=255)
                                                                                                                                        aaGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCCCt           >  1:1679009/1‑150 (MQ=255)
                                                                                                                                        aaGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACt           <  2:266605/150‑1 (MQ=255)
                                                                                                                                         ggCATTCTCCCCGCAGGCGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGAGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTa          >  1:263432/2‑150 (MQ=255)
                                                                                                                                           cATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACg        >  2:1532810/1‑150 (MQ=255)
                                                                                                                                           cATTCTCCCCGCAGGCGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCCGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACg        >  1:1308924/1‑150 (MQ=255)
                                                                                                                                           cATTCTCCCCGCAGGCGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCAGGCAGCACTACg        >  1:517722/1‑150 (MQ=255)
                                                                                                                                               ctcCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCCGCTCCGTAAAACCATTTCCAATCTGCTCTGGCAGCCCTACGCAGa    >  2:1497756/1‑150 (MQ=255)
                                                                                                                                                 cccGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAATTGCCGGCCCATGCTTATTGCCAGCGCGGTAAACCCATTACCAAACTGCTCTGGCGGCCCTAGCCAGaaa  >  2:410402/1‑150 (MQ=255)
                                                                                                                                                 cccGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGaaa  <  1:1115030/150‑1 (MQ=255)
                                                                                                                                                 cccGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGaaa  <  1:805524/150‑1 (MQ=255)
                                                                                                                                                 cccGCCAGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGaaa  <  1:232625/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GAACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACAT‑AAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAA  >  NC_000913/2711980‑2712273

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑