Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi2,892,0910CA32.0% 199.2 / 118.1 122intergenic (‑193/‑123)cysJ/queDsulfite reductase, alpha subunit, flavoprotein/6‑pyruvoyl tetrahydrobiopterin synthase (PTPS)
Reads supporting (aligned to +/- strand):  ref base C (44/39);  new base A (21/18);  total (65/57)
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.73e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 123 total aligned reads displayed.

AGCGCCCATAAAACAGAAGAGAAGGTAAGGTTAACGGGGCAAACGGTGTGGATTAAAGACGGGATAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAA  >  NC_000913_3_pae_pgi/2891943‑2892237
                                                                                                                                                    |                                                                                                                                                  
aGCGCCCATAAAACAGAAGAGAAGGTAAGGTTAACGGGGCAAACGGTGTGGATTAAAGACGGGATAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACa                                                                                                                                                   >  3:406503/1‑150 (MQ=255)
  cgcCCATAAAACAGAAGAGAAGGTAAGGTTAACGGGGCAAACGGTGTGGATTAAAGACGGGATAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACaaa                                                                                                                                                 <  1:713917/150‑1 (MQ=255)
      cATAAAACAGAAGAGAAGGTAAGGTTAACGGGGCAAACGGTGTGGATTAAAGACGGGATAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAAAAAATGaa                                                                                                                                             >  1:824113/1‑150 (MQ=255)
       aTAAAACAGAAGAGAAGGTAAGGTTAACGGGGCAAACGGTGTGGATTAAAGACGGGATAGCGATAACTAATAACCAAATCTACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAAAGAAAAAAAAATGAAt                                                                                                                                            <  2:998439/150‑1 (MQ=255)
       aTAAAACAGAAGAGAAGGTAAGGTTAACGGGGCAAACGGTGTGGATTAAAGACGTATAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAAAAAATGAAtt                                                                                                                                           >  2:995110/1‑150 (MQ=255)
       aTAAAACAGAAGAGAAGGTAAGGTTAACGGGGCAAACGGTGTGGATTAAAGACGTATAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAAAAAATGAAtt                                                                                                                                           >  2:995180/1‑150 (MQ=255)
       aTAAAACAGAAGAGAAGGTAAGGTTAACGGGGCAAACGGTGTGGATTAAAGAAAGGATAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATTGATTAACCAACAAAAAAAAAAAAAAATGAAt                                                                                                                                            <  2:998762/150‑1 (MQ=255)
            aCAGAAGAGAAGGTAAGGTTAACGGGGCAAACGGTGTGGATTAAAGACGGGATAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGc                                                                                                                                       <  4:614150/150‑1 (MQ=255)
              agaagagaagGTAAGGTTAACGGGGCAAACGGTGTGGATTAAAGACGGGATAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCa                                                                                                                                     >  4:449210/1‑150 (MQ=255)
                 agagaagGTAAGGTTAACGGGGCAAACGGTGTGGATTAAAGACGGGATAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATc                                                                                                                                  <  2:390819/150‑1 (MQ=255)
                 agagaagGTAAAGTTAACGGGGCAAACGGTGTGGATTCAAGACGTGATAGCGTTATCTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAAAAAATGAATTTAGcaatac                                                                                                                                  >  4:612050/1‑145 (MQ=255)
                     aaggtaaggtTAACGGGGCAAACGGTGTGGATTAAAGACGGGATAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTa                                                                                                                              >  4:159515/1‑150 (MQ=255)
                     aaggtaaggtTAACGGGGCAAACGGTGTGGATTAAAGACGGGATAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAAAAAATGAATTTAGCCAATCATTa                                                                                                                              >  1:906872/1‑150 (MQ=255)
                             gtTAACGGGGCAAACGGTGTGGATTAAAGACGGGATAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAAt                                                                                                                      <  4:588706/150‑1 (MQ=255)
                             gtTAACGGGGCAAACGGTGTGGATTAAAGACGGGATAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAAAAAATGAATTTAGCCAATCATTAAGATAAAt                                                                                                                      >  3:274446/1‑150 (MQ=255)
                              ttAACGGGGCAAACGGTGTGGATTAAAGACGGGATAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATc                                                                                                                     >  1:338621/1‑150 (MQ=255)
                                 aCGGGGCAAACGGTGTGGATTAAAGACGGGATAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGc                                                                                                                  >  4:480576/1‑150 (MQ=255)
                                 aCGGGGCAAACGGTGTGGATTAAAGACGGGATAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGc                                                                                                                  >  1:123799/1‑150 (MQ=255)
                                   ggggCAAACGGTGTGGATTAAAGACGGGATAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGa                                                                                                                >  1:458568/1‑150 (MQ=255)
                                   ggggCAAACGGTGTGGATTAAAGACGGGATAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGa                                                                                                                >  2:131444/1‑150 (MQ=255)
                                        aaaCGGTGTGGATTAAAGACGGGATAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTg                                                                                                           <  3:159515/150‑1 (MQ=255)
                                         aaCGGTGTGGATTAAAGACGGGATAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTgc                                                                                                          >  1:1145179/1‑150 (MQ=255)
                                           cGGTGTGGATTAAAGACGGGATAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAAAAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTgcgc                                                                                                        >  1:359188/1‑150 (MQ=255)
                                               gtgGATTAAAGACGGGATAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACGACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGcaaca                                                                                                    >  2:746828/1‑150 (MQ=255)
                                                 ggATTAAAGACGGGATAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACaaa                                                                                                  <  1:37092/150‑1 (MQ=255)
                                                   aTTAAAGACGGGATAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTaa                                                                                                                             >  3:205612/1‑121 (MQ=255)
                                                   aTTAAAGACGGGATAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTaa                                                                                                                             <  4:205612/121‑1 (MQ=255)
                                                    ttAAAGACGGGATAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTc                                                                                               <  3:253387/150‑1 (MQ=255)
                                                          aCGGGATAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTa                                                                                         >  2:181899/1‑150 (MQ=255)
                                                            gggATAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGa                                                                                       <  3:296673/150‑1 (MQ=255)
                                                            gggATAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGa                                                                                       <  3:480576/150‑1 (MQ=255)
                                                             ggATAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAAAAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGAt                                                                                      >  4:444215/1‑150 (MQ=255)
                                                              gATAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTaa                                                                                                                             >  4:196020/1‑110 (MQ=255)
                                                              gATAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTaa                                                                                                                             <  3:196020/110‑1 (MQ=255)
                                                              gATAGCGATAAATAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAAAATAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATa                                                                                     <  2:1092452/150‑1 (MQ=255)
                                                               aTAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTgcg                                                                                                         <  2:163159/129‑1 (MQ=255)
                                                                 aGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATCCATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATg                                                                                  <  1:855413/150‑1 (MQ=255)
                                                                   cGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCg                                                                                >  1:408206/1‑150 (MQ=255)
                                                                   cGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAAAAAATGAATTTAGCCAATCATTAAGATAAATCAGCGAt                                                                                                               >  4:616485/1‑119 (MQ=255)
                                                                   cGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAAAAAATGAATTTAGCCAATCATTAAGATAAATCAGCGAt                                                                                                               <  3:616485/119‑1 (MQ=255)
                                                                   cGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAAAAAACGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCg                                                                                >  2:824909/1‑150 (MQ=255)
                                                                    gATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGa                                                                               >  2:852026/1‑150 (MQ=255)
                                                                     aTAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAAAAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTa                                                                                         >  4:596315/1‑139 (MQ=255)
                                                                       aaCTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAAAAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGaaaa                                                                            <  4:86606/150‑1 (MQ=255)
                                                                         cTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGaaaaaa                                                                          <  1:652453/150‑1 (MQ=255)
                                                                          taataaCCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAAc                                                                         >  1:491804/1‑150 (MQ=255)
                                                                          taataaCCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAAc                                                                         >  2:534239/1‑150 (MQ=255)
                                                                           aataaCCAAATCAACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACa                                                                        <  3:720577/150‑1 (MQ=255)
                                                                             taaCCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGc                                                                      <  2:1122260/150‑1 (MQ=255)
                                                                                ccAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCCAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCtt                                                                   >  1:671461/1‑150 (MQ=255)
                                                                                   aaTCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAAAAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTcc                                                                <  2:906872/150‑1 (MQ=255)
                                                                                    aTCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCg                                                               <  4:56214/150‑1 (MQ=255)
                                                                                     tCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCgg                                                              >  2:156368/1‑150 (MQ=255)
                                                                                       gACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGtttt                                                                                          >  2:477811/1‑120 (MQ=255)
                                                                                       gACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGtttt                                                                                          <  1:477811/120‑1 (MQ=255)
                                                                                             ccTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCTAATGATTAAGATATATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTAc                                                      <  2:744234/150‑1 (MQ=255)
                                                                                              cTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTAcg                                                     <  1:534239/150‑1 (MQ=255)
                                                                                               tGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAAAAAATGAATTTAGCCAATCATTAAGTTAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTAcgg                                                    <  1:995180/150‑1 (MQ=255)
                                                                                                    tAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGtt                                               <  2:338621/150‑1 (MQ=255)
                                                                                                    tAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAAAAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGtt                                               <  4:274446/150‑1 (MQ=255)
                                                                                                     aGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGttt                                              >  3:332167/1‑150 (MQ=255)
                                                                                                       cAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTa                                            >  1:840060/1‑150 (MQ=255)
                                                                                                       cAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTa                                            >  1:991725/1‑150 (MQ=255)
                                                                                                         aTAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGaaaa                                                                            >  3:510373/1‑116 (MQ=255)
                                                                                                         aTAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGaaaa                                                                            <  4:510373/116‑1 (MQ=255)
                                                                                                              cTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCt                                     <  3:670484/150‑1 (MQ=255)
                                                                                                              cTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCt                                     <  1:1044384/150‑1 (MQ=255)
                                                                                                                 tAACAACAAAATAGATTAACCAACCTAATGAAAAAAAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGtt                                               <  1:286664/137‑1 (MQ=255)
                                                                                                                 tAACAACAAAATAGATTAACCAACCTAATGAAAAAAAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGtt                                               >  2:286664/1‑137 (MQ=255)
                                                                                                                   acaacaAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTaga                                <  2:870154/150‑1 (MQ=255)
                                                                                                                     aacaaAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAgaga                              >  1:151121/1‑150 (MQ=255)
                                                                                                                     aacaaAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAgaga                              >  4:244543/1‑150 (MQ=255)
                                                                                                                      acaaAATAGATTAACCAACCTAATGAAAAAAAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGaa                             >  2:228047/1‑150 (MQ=255)
                                                                                                                         aaaTAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGttt                                                                                           <  1:952453/85‑1 (MQ=255)
                                                                                                                          aaTAGATTAACCAACCTAATGAAAAAAAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTa                         >  1:882625/1‑150 (MQ=255)
                                                                                                                           aTAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTat                        <  3:91672/150‑1 (MQ=255)
                                                                                                                           aTAGATTAACCAACCTAATGAAAAAAAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTat                        <  3:444215/150‑1 (MQ=255)
                                                                                                                            tAGATTAACCAACCTAATGAAAAAAAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTatg                       <  1:228047/150‑1 (MQ=255)
                                                                                                                              gATTAACCAACCTAATGAAAAAAAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTc                                                                 <  2:319974/106‑1 (MQ=255)
                                                                                                                               aTTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTatgatg                    >  3:163128/1‑150 (MQ=255)
                                                                                                                                 tAACCAACCTAATGAAAAAAAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTc                  <  2:359188/150‑1 (MQ=255)
                                                                                                                                  atccaaccTAATGAACAAAAAATGAATTTAGCCAATAATTAAGCTAAATCAGCGACTTAGCGCAACAAAGTCGCTTTTGATAATTCTAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTCGAGAAATTATGATGTcc                 <  3:612050/148‑1 (MQ=255)
                                                                                                                                     caaccTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTccacc              >  4:47204/1‑150 (MQ=255)
                                                                                                                                        ccTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACg           >  3:720983/1‑150 (MQ=255)
                                                                                                                                        ccTAATGAAAAAAAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACg           >  1:1058431/1‑150 (MQ=255)
                                                                                                                                         cTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAg                                                                                        >  2:300094/1‑72 (MQ=255)
                                                                                                                                         cTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAg                                                                                        <  1:300094/72‑1 (MQ=255)
                                                                                                                                         cTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGt          <  1:859181/150‑1 (MQ=255)
                                                                                                                                         cTAATGAAAAAAAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGt          >  2:1066296/1‑150 (MQ=255)
                                                                                                                                          tAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGtt         >  3:345008/1‑150 (MQ=255)
                                                                                                                                          tAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGCTACTCTACGGCGGTCTTCTCGTCATGTAGAGAAATTATGATGTCCACCACGtt         >  1:817896/1‑150 (MQ=255)
                                                                                                                                           aaTGAAAAAAAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCATGTTa        >  4:17109/1‑150 (MQ=255)
                                                                                                                                             tGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTAtt      >  2:690858/1‑150 (MQ=255)
                                                                                                                                             tGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTAtt      >  1:227047/1‑150 (MQ=255)
                                                                                                                                              gAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTAttt     <  4:406503/150‑1 (MQ=255)
                                                                                                                                              gAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTAttt     >  1:380080/1‑150 (MQ=255)
                                                                                                                                              gAAAAAAAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAg                                                                       <  4:150923/84‑1 (MQ=255)
                                                                                                                                               aaaaaCAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTGTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTa    >  2:721121/1‑150 (MQ=255)
                                                                                                                                               aaaaaCAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGtt                                               >  2:646560/1‑107 (MQ=255)
                                                                                                                                                 aaaCAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTaaa  <  1:758419/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AGCGCCCATAAAACAGAAGAGAAGGTAAGGTTAACGGGGCAAACGGTGTGGATTAAAGACGGGATAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAA  >  NC_000913_3_pae_pgi/2891943‑2892237

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: