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

GAACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAA  >  NC_000913/2711980‑2712273
                                                                                                                                                     |                                                                                                                                                 
gAACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGcc                                                                                                                                                  <  1:1601591/150‑1 (MQ=255)
         ccccTTTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCTGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGgag                                                                                                                                         <  2:1170479/149‑1 (MQ=255)
         aCCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACTGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGgag                                                                                                                                         <  1:339936/150‑1 (MQ=255)
         aCCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGgag                                                                                                                                         >  1:266605/1‑150 (MQ=255)
          cccATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGTATAGTTCTGGTCTGGCCCACGAATGACTAAGCGCTCGCTAACGCATTCTCCCGCCGGCGGGgaga                                                                                                                                        >  2:643995/1‑150 (MQ=255)
           ccATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGaa                                                                                                                                       >  2:495998/1‑150 (MQ=255)
           ccATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGaa                                                                                                                                       <  1:582219/150‑1 (MQ=255)
            cATTCCGGTCGGTCGATCCTTTCCCGCGGCAATCATTACATAAACAGATAAAATTCCTGGGTCAGCGCCGTATTCGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAATGCATTCTCCCGCCTGCGGGGAGAAt                                                                                                                                      <  2:833454/150‑1 (MQ=255)
            cATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAAt                                                                                                                                      <  1:1580438/150‑1 (MQ=255)
            cATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCACACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAAt                                                                                                                                      >  2:1332765/1‑150 (MQ=255)
            cATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGGCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATGCTCCCGCCTGCGGGGAGAAt                                                                                                                                      <  1:1032832/150‑1 (MQ=255)
             aTTCCTGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCATCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATg                                                                                                                                     <  2:1640356/150‑1 (MQ=255)
             aTTCCGGTCCGTCGATCCTTACCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATg                                                                                                                                     <  2:393237/150‑1 (MQ=255)
              ttCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAAACAGATAAAAAGCCTGGGGCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCCGCAGGCGGGAGAATg                                                                                                                                     <  1:1196655/150‑1 (MQ=255)
               tCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGCAATAGTGCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCCGCAGGCGGGAGAATGct                                                                                                                                   <  1:1111880/150‑2 (MQ=255)
                ccGGTCCGGCGATCCTTTCCCCCGGCAATCATTACCTAAACAGATAAAACGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCca                                                                                                                                  <  1:655881/150‑1 (MQ=255)
                 cGGTCCGTCGCTCCTTTCCCCCGGCAATCATTACATAAACAGATAAACAGCCTTGGGCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCCGCAGGCGGGAGAATGcttt                                                                                                                                 <  1:645678/150‑4 (MQ=255)
                 cGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCcag                                                                                                                                 <  2:642944/150‑1 (MQ=255)
                  ggTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGTGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCcagc                                                                                                                                <  2:471380/150‑1 (MQ=255)
                  ggTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCcagc                                                                                                                                <  1:965446/150‑1 (MQ=255)
                  ggTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCcagc                                                                                                                                >  2:1038085/1‑150 (MQ=255)
                   gTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCTGCGCCGTATACGCTACGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCcagca                                                                                                                               <  2:1549721/150‑1 (MQ=255)
                   gTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCcagca                                                                                                                               <  1:1011735/150‑1 (MQ=255)
                      cGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGAGGGGAGAATTCCagcagcc                                                                                                                            >  2:1532572/1‑149 (MQ=255)
                                 tCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGTGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGc                                                                                                                 <  2:1545337/150‑1 (MQ=255)
                                   ccccGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGc                                                                                                                    <  2:1156318/145‑1 (MQ=255)
                                   ccccGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGc                                                                                                                    >  1:1156318/1‑145 (MQ=255)
                                    cccGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGt                                                                                                              >  2:794653/1‑150 (MQ=255)
                                    cccGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGt                                                                                                              >  2:117502/1‑150 (MQ=255)
                                    cccGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGt                                                                                                              >  2:787875/1‑150 (MQ=255)
                                        gCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCTGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTcgcg                                                                                                          <  1:434973/150‑1 (MQ=255)
                                           aTCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCtt                                                                                                       >  2:1596682/1‑150 (MQ=255)
                                                 aCATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGtttt                                                                                                 <  1:1719/150‑1 (MQ=255)
                                                 aCATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGtttt                                                                                                 >  2:1330650/1‑150 (MQ=255)
                                                  cATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGTGAATGCCAGCAGCACCCGATCCGGCAGACGCGCTTCGTTTTc                                                                                                >  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)
                                                                                                                                                     |                                                                                                                                                 
GAACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAA  >  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: