Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_pgi 1,677,492 Δ1 bp coding (412/1533 nt) pntA ← pyridine nucleotide transhydrogenase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi1,677,4910C.96.2% 362.1 / 7.5 79coding (413/1533 nt)pntApyridine nucleotide transhydrogenase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base . (37/39);  total (38/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.74e-01

GACTGTTTGCTGCGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCA  >  NC_000913_3_pae_pgi/1677345‑1677635
                                                                                                                                                  |                                                                                                                                                
gaCTGTTTGCTGCGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAAcggc                                                                                                                                               <  1:923211/150‑4 (MQ=255)
         ctgcGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGc                                                                                                                                     <  2:336750/150‑1 (MQ=255)
            cgcCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCGAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCAt                                                                                                                                  >  3:256499/1‑150 (MQ=255)
             gcCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATc                                                                                                                                 >  1:6199/1‑150 (MQ=255)
              gcAATGGCTGCTCGACTTGCAACACCCGCACCAATCGCCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCg                                                                                                                                <  4:406068/149‑1 (MQ=255)
              ccAATGGCGGCCAGACCTGCAACACCCGCGCCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTGAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCg                                                                                                                                <  2:746301/150‑1 (MQ=255)
              ccAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCg                                                                                                                                >  2:32711/1‑150 (MQ=255)
               cAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGa                                                                                                                               >  1:1105664/1‑150 (MQ=255)
               cAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGa                                                                                                                               >  4:9056/1‑150 (MQ=255)
                aaTGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAg                                                                                                                              >  3:495061/1‑150 (MQ=255)
                  tGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCt                                                                                                                            <  2:1105664/150‑1 (MQ=255)
                    gcggcCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTa                                                                                                                          >  3:434328/1‑150 (MQ=255)
                      ggcCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGt                                                                                                                        >  1:1097752/1‑150 (MQ=255)
                      ggcCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGt                                                                                                                        >  1:506198/1‑150 (MQ=255)
                          aGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGt                                                                                                                    >  1:248225/1‑150 (MQ=255)
                           ggccTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTc                                                                                                                   <  2:346241/148‑1 (MQ=255)
                           gACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTc                                                                                                                   <  1:399712/150‑1 (MQ=255)
                           gACCTGCAACACCCGCACCAATCACCAGCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTc                                                                                                                   <  1:1143923/150‑1 (MQ=255)
                              cTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAg                                                                                                                <  2:546202/150‑1 (MQ=255)
                                         gCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCGAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTgcgc                                                                                                     <  4:256499/150‑1 (MQ=255)
                                         gCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTgcgc                                                                                                     >  2:987201/1‑150 (MQ=255)
                                         gCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTgcgc                                                                                                     <  2:271244/150‑1 (MQ=255)
                                          cACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGcgcg                                                                                                    >  3:176587/1‑150 (MQ=255)
                                           aCCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGt                                                                                                   <  2:1097752/150‑1 (MQ=255)
                                           aCCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGt                                                                                                   <  2:506198/150‑1 (MQ=255)
                                            ccAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGt                                                                                                   <  1:718210/150‑1 (MQ=255)
                                               aTCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGAt                                                                                                                                           >  3:637553/1‑106 (MQ=255)
                                               aTCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGAt                                                                                                                                           <  4:637553/106‑1 (MQ=255)
                                                tCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGAt                                                                                              >  3:652884/1‑150 (MQ=255)
                                                  aCCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATAc                                                                                            >  3:454059/1‑150 (MQ=255)
                                                   ccATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATAcg                                                                                           >  4:599343/1‑150 (MQ=255)
                                                    cATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATAcgc                                                                                          <  2:248225/150‑1 (MQ=255)
                                                     aTCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGt                                                                                                                        >  4:340515/1‑119 (MQ=255)
                                                     aTCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGt                                                                                                                        <  3:340515/119‑1 (MQ=255)
                                                     aTCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATAcgcg                                                                                         >  4:487423/1‑150 (MQ=255)
                                                     aTCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATAcgcg                                                                                         <  3:9056/150‑1 (MQ=255)
                                                      tCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCgg                                                                                        >  4:654308/1‑150 (MQ=255)
                                                          ttttGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGtt                                                                                                                                        <  4:34927/98‑1 (MQ=255)
                                                          ttttGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGtt                                                                                                                                        >  3:34927/1‑98 (MQ=255)
                                                           tttGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACag                                                                                   <  1:192035/150‑1 (MQ=255)
                                                               ccGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTc                                                                               <  2:1126105/150‑1 (MQ=255)
                                                                       aCTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCa                                                                       <  4:728/150‑1 (MQ=255)
                                                                          ttCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCa                                                                    <  4:319373/150‑1 (MQ=255)
                                                                             ccGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACgg                                                                 >  4:386760/1‑150 (MQ=255)
                                                                                 ccGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCac                                                             <  2:103645/150‑1 (MQ=255)
                                                                                   gCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCacgt                                                           >  1:1113708/1‑150 (MQ=255)
                                                                                      gTAATTTGCCCGGTAAAGAAGCGCCCAAATGCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCacgttac                                                        <  2:288752/150‑1 (MQ=255)
                                                                                       tAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTc                                                                                                                   >  4:269442/1‑90 (MQ=255)
                                                                                       tAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTc                                                                                                                   <  3:269442/90‑1 (MQ=255)
                                                                                        aaTTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCacgttacgt                                                      >  4:153588/1‑150 (MQ=255)
                                                                                            tGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCg                                                  <  2:79132/150‑1 (MQ=255)
                                                                                               ccGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCaa                                               <  1:76614/150‑1 (MQ=255)
                                                                                                cGGTAACGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTcc                                                   >  4:558030/1‑146 (MQ=255)
                                                                                                cGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTcc                                                   <  3:558030/146‑1 (MQ=255)
                                                                                                       gttacgcCCAAATCCAAGTGCAGCTTCAACAATCGCGCGATAA‑CGGCGATGTCGCCCGGCGAGCTTCGTGCGGCCAGCGATTGTCCGCGTGAGATACGCGGCACAGAGTACATCGCCTTCACGGTCACGTTACGTTCCCCCAGATGTTGc                                       >  2:467953/5‑150 (MQ=255)
                                                                                                       gaagCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGt                                                                                                                        >  2:713643/1‑69 (MQ=255)
                                                                                                       gaagCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGt                                                                                                                        <  1:713643/69‑1 (MQ=255)
                                                                                                       gaagCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGc                                       <  4:309001/150‑1 (MQ=255)
                                                                                                          gcgcCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCAtt                                    <  4:82208/150‑1 (MQ=255)
                                                                                                           accccAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTa                                   <  4:410079/148‑1 (MQ=255)
                                                                                                             cccAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAAt                                 <  1:482928/150‑1 (MQ=255)
                                                                                                             cccAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAAt                                 >  4:731412/1‑150 (MQ=255)
                                                                                                              ccAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAAtt                                >  3:536783/1‑150 (MQ=255)
                                                                                                               cAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTc                               >  1:487246/1‑150 (MQ=255)
                                                                                                               cAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTc                               >  2:800693/1‑150 (MQ=255)
                                                                                                                   ttCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGa                           >  3:512260/1‑150 (MQ=255)
                                                                                                                     cATGTGCCGCTTCAACAATGGGGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGAtt                         <  4:25332/150‑1 (MQ=255)
                                                                                                                      aTGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTc                        >  1:831417/1‑150 (MQ=255)
                                                                                                                        gtgCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTg                      <  2:1100342/150‑1 (MQ=255)
                                                                                                                         tgCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTgc                     >  1:539372/1‑150 (MQ=255)
                                                                                                                             gCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAg                 >  4:2531/1‑150 (MQ=255)
                                                                                                                                 caacaaTGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCa             <  2:487246/150‑1 (MQ=255)
                                                                                                                                  aacaaTGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCAGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAg            >  3:623454/1‑150 (MQ=255)
                                                                                                                                     aaTGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTgc                                                                                                       <  2:513076/56‑1 (MQ=255)
                                                                                                                                     aaTGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTgc                                                                                                       >  1:513076/1‑56 (MQ=255)
                                                                                                                                        ggCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATaaaa      <  3:113123/150‑1 (MQ=255)
                                                                                                                                        ggCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATaaaa      <  4:454059/150‑1 (MQ=255)
                                                                                                                                        ggCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATaaaa      <  3:654308/150‑1 (MQ=255)
                                                                                                                                         gcgcGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAAc     >  4:600907/1‑150 (MQ=255)
                                                                                                                                            cgATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCa  <  4:176587/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GACTGTTTGCTGCGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCA  >  NC_000913_3_pae_pgi/1677345‑1677635

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