Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_pae_pgi 3,439,589 G→C 6.1% intergenic (‑80/+27) yhdN ← / ← rplQ DUF1992 family protein/50S ribosomal subunit protein L17

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi3,439,5890GC 6.1% 173.1 / 5.4 68intergenic (‑80/+27)yhdN/rplQDUF1992 family protein/50S ribosomal subunit protein L17
Reads supporting (aligned to +/- strand):  ref base G (33/30);  new base C (2/2);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.27e-01

TTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACG  >  NC_000913_3_pae_pgi/3439442‑3439732
                                                                                                                                                   |                                                                                                                                                
ttATCAAACTCACCTTTGTCCTGTCCTTCTTCTATATTGCGCTCTGCCCACTGGTTAAGTAACCACATAAATACTCCAGGGGAAAAACAAAAAGTTTTATGCTATTGTAGATAAGTGGGGCTACTACGTGTATAAAAAAACCAGCCggg                                                                                                                                                 <  2:213093/149‑1 (MQ=255)
  aTCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGc                                                                                                                                               >  1:395672/1‑149 (MQ=255)
      aaCTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATACACAGAAGGAGTACAGCTATTGTAGATCAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGt                                                                                                                                           <  1:443353/149‑1 (MQ=255)
       aCTCTCCTTACTCTTGCACTTTTTGTATTTGGCCCTCTGCCCACGGGTCAATTAACCACATAAATAGCCCAGGGGATGAACCAACAGAGTTTACCTTATGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGtt                                                                                                                                          <  2:418502/149‑1 (MQ=255)
         tCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTGAAGTAACCACATAAATACTTGAGGGGATGAACAAAAATAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTAGAAAAAAACCCGCCGGGGCGGGtttt                                                                                                                                        >  1:907167/1‑149 (MQ=255)
          cACCTTTCGCTTTCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGttttt                                                                                                                                       <  1:520961/149‑1 (MQ=255)
           aCCTTTCGCTTGCGCTTGTGTTATATGGCGCTCTGTCCACTGGTCAAGTAACCACATAAATACTCCAGTGGATGAACAAAAAGAGGACAGCTATTGTAGATAAGTGGGGATACTACGGCTTCAAATTAACCCGCCGGGGCGGGtttttt                                                                                                                                      <  2:972396/149‑1 (MQ=255)
           aCCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGtttttt                                                                                                                                      <  1:622989/149‑1 (MQ=255)
              tttCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGAGAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGAGGCGGGgttttttgc                                                                                                                                   >  2:243831/1‑147 (MQ=255)
                      gcgcTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTc                                                                                                                           <  1:96680/149‑1 (MQ=255)
                        gcTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAg                                                                                                                         >  1:811082/1‑149 (MQ=255)
                            cTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTa                                                                                                                     >  1:469339/1‑149 (MQ=255)
                             tGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGTTTTAAAAAAACCCGCCGGGGCGGGTTTTTTTACATTGCTTCAgcctac                                                                                                                    >  1:439607/1‑144 (MQ=255)
                              gCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTAct                                                                                                                   >  1:952048/1‑149 (MQ=255)
                              gCTATATGGCGCTCTGACAACTGGTCAAGAAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTAct                                                                                                                   <  1:289262/149‑1 (MQ=255)
                               cTATATGGCGCCCTGCCCACTGGTCAAGTAACCACATAAATAGTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTGTTTTACGTCGCTGCAGACTAgtc                                                                                                                  >  2:605253/1‑149 (MQ=255)
                                tataTGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTActct                                                                                                                 <  1:738104/149‑1 (MQ=255)
                                tataTGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGAATACTACGGGTATAACAAAACCCCCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTActct                                                                                                                 >  2:950622/1‑149 (MQ=255)
                                  tatGGCGCTGTGCCCATGGTTCCAGTTCCCAGCGAAATACACCACGGGACGTACAAAAAGTGTCCATCTATAGTTGTTAAGGTGGGTTATGTCTGGTTTAAAAAGACCCTCTGGGTCAGGTTTTTTTACGTTGCTTCCGATTACTCTgc                                                                                                               <  2:950793/149‑1 (MQ=255)
                                   atGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGCGTAAAGCTATTGTAGATAAGTGGGGATCCTACGGGTATAAAAAAACGCGCCGGGGCGGGTTTGTTTTCGTTGCTTGTGCTTACTCTgca                                                                                                              >  1:410363/1‑149 (MQ=255)
                                    tGGCGCTCCGCCAACTGGTCAAATAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTgcag                                                                                                             <  2:717998/149‑1 (MQ=255)
                                     ggCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTgcagc                                                                                                            <  1:857518/149‑1 (MQ=255)
                                       cgcgcggCCCAGTGGTCAAGTAGCCACATAAATACTCCAGGGGATGAACAAAAAGAGGACAGTTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGcagcag                                                                                                          <  2:426169/143‑1 (MQ=255)
                                              cccACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATATAAAAACCCGCCGGGGCGGGGTTTTTTACGTTGCTTCAGATTAGTCTGCAGCAGATTCTgt                                                                                                   >  1:435625/1‑148 (MQ=255)
                                              cccACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGc                                                                                                   >  2:38147/1‑149 (MQ=255)
                                                    ggTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCgg                                                                                                                                             >  2:185960/1‑101 (MQ=255)
                                                    ggTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCgg                                                                                                                                             <  1:185960/101‑1 (MQ=255)
                                                        aaGTAACCACATAAATACTTCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGaa                                                                                         <  1:734474/149‑1 (MQ=255)
                                                        aaGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGAACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAACCCCGCCGGGGCGGGTGTTTTTACGTTGCTTCAGATTACTCTGCACCAGCTTCTACTTTCTATGaa                                                                                         >  2:831020/1‑149 (MQ=255)
                                                        aaGTAACCACATAAATACTCAAGGGGATGAACACAACGAGTACCGCAATAAAAGATAGGAGGGGATACTAGGGGTATAAAAAAACCCGCCGGGGCGCGTTTTTTTGCTTTACTGCACTTTCATCTGCAGCACCTCGCGCTATCTCTgca                                                                                         >  2:369571/1‑147 (MQ=255)
                                                           tAACCACATAAATACTCCTGGGGATGTATAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGGATAAAAAAACCCGCCGGGGCGGGTTTTTTTAGGTTGCTTGAGATTGGTCTGCAGCAGCTTTTGGTTTCTGTGataga                                                                                      >  2:703134/1‑145 (MQ=255)
                                                           tAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAACCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGa                                                                                      <  2:811082/149‑1 (MQ=255)
                                                               cacaTAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGCTACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGATCCTGCTTTCTCTGATCGATCaa                                                                                  >  1:229772/1‑149 (MQ=255)
                                                                 caTAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGGTTTTTTACGTTGCTTCAGATTACTCTGCCGCAGCTTCTGCTTTCTCTGAACGATCAAcc                                                                                >  2:742385/1‑149 (MQ=255)
                                                                      aTACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTc                                                                           >  1:242820/1‑149 (MQ=255)
                                                                            cAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTACTTTCTCTGAACGATCAACCAGCTCGATGTa                                                                     >  1:838568/1‑149 (MQ=255)
                                                                               gggATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGAATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCAGAACGATCAACCAGCTCGATGTCAGc                                                                  >  2:976829/1‑149 (MQ=255)
                                                                                  aTGAACAAAAAGATCACAGCTATTGGAGATAAGTGGGTATACTACGGGTATAAAAAATCCCGCCAGGGCGGGTTCTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCAt                                                               <  2:935609/149‑1 (MQ=255)
                                                                                     aaCAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCgg                                                           >  1:766574/1‑149 (MQ=255)
                                                                                     aaCAAAAAGAGTACAGCTATTGTAGATAAGTGAGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCATATTCTGCTTTCTCTTAACGATCAACCAGCTCGATGTGAGCCATCgc                                                           >  2:265125/1‑148 (MQ=255)
                                                                                        aaaaaGAGTACGGGTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTTCTGCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGgcgc                                                        <  1:641580/149‑1 (MQ=255)
                                                                                         aaaaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTGTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAAACAGCTCAATGTAAACCATCTgcgca                                                       >  1:164344/1‑148 (MQ=255)
                                                                                         aaaaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCCCGGCGGGTTTTTTTa                                                                                                                                    <  1:672994/73‑1 (MQ=37)
                                                                                         aaaaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCCCGGCGGGTTTTTTTa                                                                                                                                    >  2:672994/1‑73 (MQ=37)
                                                                                             gagTACAGATATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGt                                                   <  2:333246/149‑1 (MQ=255)
                                                                                                tatctcTATTGTAGATAAGTGGGGATACTATGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGAAGCTTATGCTTTCGCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGc                                                <  2:5198/144‑1 (MQ=255)
                                                                                                          gTAGATAAATGGGGAAACTACGGTTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGaa                                      <  1:997044/149‑1 (MQ=255)
                                                                                                           gaGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGATTCTGCTTTCTCTGAACGATCAACCAGCTCGAAGTAAGCCATCGACGCGTTGTCGCCGGCACGGAAg                                     >  1:801651/2‑149 (MQ=255)
                                                                                                             gATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTGGGCTGCACGGAAggc                                   >  1:632248/1‑147 (MQ=255)
                                                                                                              aTAAGTGGGGATACTACGGGTATAAAAAAACCCCCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATTTAAGCCATCGGCGCGTTGTCGCGGCCACAGAAGcaa                                  >  1:344122/1‑147 (MQ=255)
                                                                                                                     gggATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCCAGCAGGTTTCGCCTGCACGGAAGCCACACTTCa                           <  2:469339/149‑1 (MQ=255)
                                                                                                                        aTACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACTATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGTACGCAACACTTCAGaa                        <  2:907167/149‑1 (MQ=255)
                                                                                                                        aTACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGaa                        <  1:38147/149‑1 (MQ=255)
                                                                                                                        aTACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGaa                        <  2:661145/149‑1 (MQ=255)
                                                                                                                            tacGGTTATAAAAAAACCCGCCGGGGCGCGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACg                    <  1:183939/149‑1 (MQ=255)
                                                                                                                             gggggTACAATAGAACCCGGCGGGGCGGGTTTTTTTACGTTCGTTCGGATTACTCTGTAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGa                   <  2:77700/147‑1 (MQ=255)
                                                                                                                             atgggTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGa                   <  1:220920/147‑1 (MQ=255)
                                                                                                                             acGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGa                   <  2:242820/149‑1 (MQ=255)
                                                                                                                             aagtgTATAAAAAAACCCGCCGGGGCCGGTTTTTTTACTTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGa                   <  2:164344/145‑1 (MQ=255)
                                                                                                                                 gTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAgtgt               >  1:357458/1‑149 (MQ=255)
                                                                                                                                       aaaaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAaccac         >  2:365162/1‑149 (MQ=255)
                                                                                                                                       aaaaaCCCGCCGGGGCGGGTTGTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAaccac         <  1:950622/149‑1 (MQ=255)
                                                                                                                                       aaaaaCCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGt                                                                      <  1:943799/89‑1 (MQ=255)
                                                                                                                                       aaaaaCCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCACCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGt                                                                      >  2:943799/1‑89 (MQ=38)
                                                                                                                                        aaaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTAAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCATGCACGCAAGCCAAACTTCAGAATCCGAGTGTAaccacc        >  2:98221/1‑149 (MQ=255)
                                                                                                                                        aaaaCCCGCCGGGGCGGGTTTTTTGACGTTGCTTCAGATTGCTCTGGGGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTCGTCACCTGCACGGAAGCTACACTTCGGATTACGACCGTAaccccc        >  2:168195/1‑149 (MQ=255)
                                                                                                                                             ccgccgGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCAc   >  2:17368/1‑149 (MQ=255)
                                                                                                                                             ccgccgGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCAc   >  2:969001/1‑149 (MQ=255)
                                                                                                                                              cgccgGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTGAGTAAAAGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCaag  >  1:737114/1‑147 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
TTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACG  >  NC_000913_3_pae_pgi/3439442‑3439732

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