Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,466,101 G→C 5.7% noncoding (102/119 nt) PP_mr39 ← P15 ncRNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,466,1010GC 5.7% 283.4 / 5.9 87noncoding (102/119 nt)PP_mr39P15 ncRNA
Reads supporting (aligned to +/- strand):  ref base G (46/36);  new base C (3/2);  total (49/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

GCGGCAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAG  >  NC_002947/3465959‑3466246
                                                                                                                                              |                                                                                                                                                 
gcggcAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCACAGGCCAGCAAAACAAAACCCCCGGACTGGAAGCCGACCg                                                                                                                                            >  2:632850/1‑150 (MQ=255)
  ggcAGGCTGGCCCGAATGGTATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGCTTCCCGACCggg                                                                                                                                          <  1:884757/150‑1 (MQ=255)
  ggcAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGACGCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCggg                                                                                                                                          <  2:1929393/150‑1 (MQ=255)
  ggcAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGCGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAATCCCCCGGTCGGGAAGCCGACCggg                                                                                                                                          <  2:1739534/150‑1 (MQ=255)
   gcAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCgggg                                                                                                                                         <  1:1404565/150‑1 (MQ=255)
         tGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGa                                                                                                                                   <  1:1934858/150‑1 (MQ=255)
         tGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGa                                                                                                                                   >  2:1934858/1‑150 (MQ=255)
             ccGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCgggg                                                                                                                                         <  2:68781/140‑1 (MQ=255)
             ccGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCgggg                                                                                                                                         >  1:68781/1‑140 (MQ=255)
             ccGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTAt                                                                                                                               >  2:178100/1‑150 (MQ=255)
              cGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTAtt                                                                                                                              <  1:226025/150‑1 (MQ=255)
                      aTTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGtt                                                                                                                      <  1:632850/150‑1 (MQ=255)
                       ttGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTGCTAATTGCCCTAGTTCATGAATCTGTCCTCACAGACGTGTGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTc                                                                                                                     <  1:1141081/150‑1 (MQ=255)
                       ttGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGttt                                                                                                                     <  2:1536627/150‑2 (MQ=255)
                          gCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTcggc                                                                                                                  <  2:2117210/150‑1 (MQ=255)
                              aaTCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGcc                                                                                                              >  1:2201159/1‑150 (MQ=255)
                                   ccAGGCGTAGGTGTTTGCGTAGTAATAAGCGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCg                                                                                                         >  2:1963791/1‑150 (MQ=255)
                                     aGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGaa                                                                                                       >  1:184349/1‑150 (MQ=255)
                                        cGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATGCTCATGTTCGGCGGCCCCTCGAAGGt                                                                                                    <  1:1687229/150‑1 (MQ=255)
                                         gTAGGTGTTTGCGTAGTAATAAGTGGCGTCGTTCATCATTCTGTCCTCAATGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTg                                                                                                   <  2:2201159/150‑1 (MQ=255)
                                           aGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGa                                                                                                 >  2:2096480/1‑150 (MQ=255)
                                            ggTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGAc                                                                                                <  1:72076/150‑1 (MQ=255)
                                                tttGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCg                                                                                            <  2:1887656/150‑1 (MQ=255)
                                                         taataaGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCggggg                                                                                                                                        >  1:95125/1‑97 (MQ=255)
                                                         taataaGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCggggg                                                                                                                                        <  2:95125/97‑1 (MQ=255)
                                                               gTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGgcgc                                                                             <  1:715470/150‑1 (MQ=255)
                                                                  gCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCg                                                                                            >  1:1197680/1‑132 (MQ=255)
                                                                  gCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCg                                                                                            <  2:1197680/132‑1 (MQ=255)
                                                                  gCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTa                                                                          >  2:134108/1‑150 (MQ=255)
                                                                   cGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAg                                                                         <  2:184349/150‑1 (MQ=255)
                                                                       gTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGc                                                                     >  1:1280078/1‑150 (MQ=255)
                                                                         tCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTa                                                                   >  1:387565/1‑150 (MQ=255)
                                                                           aTGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGcccc                                                                                                            >  1:1009745/1‑107 (MQ=255)
                                                                           aTGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGcccc                                                                                                            <  2:1009745/107‑1 (MQ=255)
                                                                           aTGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCg                                                                                 <  1:1704358/134‑1 (MQ=255)
                                                                           aTGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCg                                                                                 >  2:1704358/1‑134 (MQ=255)
                                                                           aTGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTaaa                                                                 >  1:1804464/1‑150 (MQ=255)
                                                                           aTGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTaaa                                                                 >  1:534981/1‑150 (MQ=255)
                                                                           aTGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGACGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTaaa                                                                 >  2:2060403/1‑150 (MQ=255)
                                                                           aTGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGCTTCCCGACCGGGGGtttt                                                                                                                                    <  2:63314/83‑2 (MQ=25)
                                                                           aTGAATCTGTCCTCAAAGTGAATCGGTGTAGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGCTTCCCGACCGGGGGtttt                                                                                                                                    >  1:63314/1‑82 (MQ=21)
                                                                            tGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAAc                                                                <  2:1754959/150‑1 (MQ=255)
                                                                                      cTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAGACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCgggt                                                                                         <  2:1180642/115‑1 (MQ=255)
                                                                                      cTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTg                                                                                                   >  2:1579787/1‑105 (MQ=255)
                                                                                      cTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTg                                                                                                   <  1:1579787/105‑1 (MQ=255)
                                                                                      cTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCgggt                                                                                         >  1:1180642/1‑115 (MQ=255)
                                                                                       tCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCaaa                                                     >  2:97434/1‑150 (MQ=255)
                                                                                       tCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCaaa                                                     >  1:1434966/1‑150 (MQ=255)
                                                                                       tCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCaaa                                                     >  1:1100502/1‑150 (MQ=255)
                                                                                            gTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAAt                                                >  2:511223/1‑150 (MQ=255)
                                                                                                 tCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTg                                           <  2:696319/150‑1 (MQ=255)
                                                                                                     tgtTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAg                                       <  1:134108/150‑1 (MQ=255)
                                                                                                          gTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGcc                                  >  1:1715107/1‑150 (MQ=255)
                                                                                                           tCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCa                                 >  1:1436381/1‑150 (MQ=255)
                                                                                                            cAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAg                                                                                                      >  2:692787/1‑80 (MQ=255)
                                                                                                            cAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAg                                                                                                      <  1:692787/80‑1 (MQ=255)
                                                                                                             aaaGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCgggtgg                                                                                       <  1:1042387/94‑1 (MQ=255)
                                                                                                             aaaGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCgggtgg                                                                                       >  2:1042387/1‑94 (MQ=255)
                                                                                                                   cTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCg                         >  2:1471546/1‑150 (MQ=255)
                                                                                                                    tGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGgcgc                                                                             >  2:1551325/1‑97 (MQ=255)
                                                                                                                    tGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGgcgc                                                                             <  1:1551325/97‑1 (MQ=255)
                                                                                                                    tGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAAt                                                >  2:1588638/1‑126 (MQ=255)
                                                                                                                    tGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAAt                                                <  1:1588638/126‑1 (MQ=255)
                                                                                                                     gAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGc                       <  1:178100/150‑1 (MQ=255)
                                                                                                                         actcaaacCCCCGGTCGGCTTCCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATAGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCTACTCAGCCGGCAACg                   >  2:765231/5‑150 (MQ=38)
                                                                                                                           caaaacCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGc                 >  1:669008/1‑150 (MQ=255)
                                                                                                                           caaaacCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGAATTCAACTGATCGGTGCAAGCCCGAAGGTGCGTCCCGATGGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGc                 >  1:656457/1‑150 (MQ=255)
                                                                                                                            aaaacCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGCCGGCGGGTGGGTATCGGCGCCTAGTGGCTACACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGcc                >  2:1379218/1‑150 (MQ=255)
                                                                                                                            aaaacCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGcc                >  1:1889822/1‑150 (MQ=255)
                                                                                                                            aaaacCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGcc                >  2:2188761/1‑150 (MQ=255)
                                                                                                                               acCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTgg                                          >  2:250971/1‑121 (MQ=255)
                                                                                                                               acCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTgg                                          <  1:250971/121‑1 (MQ=255)
                                                                                                                                cccccGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCTTAGTGGGTAAACCTATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGca            >  1:286281/1‑150 (MQ=255)
                                                                                                                                 ccccGGTCGGGAAGTCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGAGCATCTGCGCCTAGTGGCTAAACCAAAAACCACACTAATATCGAGGAGGGGCCACGCAGCCGGCAACGGCCAGcac           >  1:1141128/1‑150 (MQ=255)
                                                                                                                                 ccccGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGcac           >  2:723717/1‑150 (MQ=255)
                                                                                                                                 ccccGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGcac           >  1:2018896/1‑150 (MQ=255)
                                                                                                                                  cccGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGcacc          >  1:435806/1‑149 (MQ=255)
                                                                                                                                  cccGGTCGGCTTCCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGcaca          >  2:1928314/1‑150 (MQ=255)
                                                                                                                                     ggTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAggg       <  1:707211/150‑1 (MQ=255)
                                                                                                                                        cGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTgg                                                                      <  2:1096021/84‑1 (MQ=255)
                                                                                                                                        cGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTgg                                                                      >  1:1096021/1‑84 (MQ=255)
                                                                                                                                         gggAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCAc                                >  2:1852781/1‑121 (MQ=255)
                                                                                                                                         gggAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCAc                                <  1:1852781/121‑1 (MQ=255)
                                                                                                                                         gggAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCa   >  1:1016347/1‑150 (MQ=255)
                                                                                                                                          ggAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAg  <  1:2096480/150‑1 (MQ=255)
                                                                                                                                          ggAAGCCGACCGGGGGTTTGAGTAATCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCgggt                                                                                         >  2:2233223/1‑63 (MQ=255)
                                                                                                                                          ggAAGCCGACCGGGGGTTTGAGTAATCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCgggt                                                                                         <  1:2233223/63‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
GCGGCAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAG  >  NC_002947/3465959‑3466246

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