Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,712,133 G→C 8.0% P211R (CCG→CGG)  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,1330GC 8.0% 294.8 / 13.2 98P211R (CCG→CGG) trmNtRNA1(Val) (adenine(37)‑N6)‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (45/44);  new base C (5/3);  total (51/47)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.12e-01

ACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATCCCTCT  >  NC_000913/2711989‑2712282
                                                                                                                                                 |                                                                                                                                                     
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                                                                                                                                                <  1:582219/150‑1 (MQ=255)
  ccATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGaa                                                                                                                                                >  2:495998/1‑150 (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                                                                                                                                         >  2:1038085/1‑150 (MQ=255)
         ggTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCcagc                                                                                                                                         <  1:965446/150‑1 (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                                                                                                                             >  1:1156318/1‑145 (MQ=255)
                          ccccGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGc                                                                                                                             <  2:1156318/145‑1 (MQ=255)
                           cccGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGt                                                                                                                       >  2:787875/1‑150 (MQ=255)
                           cccGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGt                                                                                                                       >  2:794653/1‑150 (MQ=255)
                           cccGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGt                                                                                                                       >  2:117502/1‑150 (MQ=255)
                               gCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCTGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTcgcg                                                                                                                   <  1:434973/150‑1 (MQ=255)
                                  aTCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCtt                                                                                                                >  2:1596682/1‑150 (MQ=255)
                                        aCATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGtttt                                                                                                          >  2:1330650/1‑150 (MQ=255)
                                        aCATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGtttt                                                                                                          <  1:1719/150‑1 (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:794653/150‑1 (MQ=255)
                                                                   aGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTaaa                                                                                <  1:787875/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                       >  2:1728130/1‑150 (MQ=255)
                                                                                                                            cTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGgcagc                       <  1:1394667/150‑1 (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)
                                                                                                                                             cTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATcc      <  2:690817/150‑1 (MQ=255)
                                                                                                                                                 ggggAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATCCctct  >  2:1695869/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
ACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGATGCGGCAGTCGCGCTTCGTTTTCCGCCACATCCGTCCGCAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATCCCTCT  >  NC_000913/2711989‑2712282

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: