Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_pgi 3,440,924 C→A G36C (GGC→TGC)  rpoA ← RNA polymerase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi3,440,9240CA95.7% 394.7 / 6.5 116G36C (GGC→TGC) rpoARNA polymerase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (2/3);  new base A (53/58);  total (55/61)
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.52e-01

Only 100 of 117 total aligned reads displayed.

GGTTGAGCAGGATTTCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGG  >  NC_000913_3_pae_pgi/3440777‑3441062
                                                                                                                                                   |                                                                                                                                          
ggTTGAGCAGGATTTCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCaaa                                                                                                                                          <  1:158786/150‑1 (MQ=255)
       cAGGATTTCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCCGAGGATGGCAAAAGCCACg                                                                                                                                   >  4:704977/1‑150 (MQ=255)
          caTTTCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGAAGTACACCATCAATCTCTACCTCGGTCACCGCGCCACACGCCATCGCTGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGctc                                                                                                                                >  2:166824/2‑150 (MQ=255)
             ttCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTaa                                                                                                                             >  3:455025/1‑150 (MQ=255)
               ccAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAg                                                                                                                           >  3:601502/1‑150 (MQ=255)
                 aGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGc                                                                                                                         <  4:386932/150‑1 (MQ=255)
                 aGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGc                                                                                                                         <  1:335931/150‑1 (MQ=255)
                      atCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAg                                                                                                                    <  3:171270/150‑1 (MQ=255)
                       tCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCTCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAgg                                                                                                                   <  1:893416/150‑1 (MQ=255)
                       tCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAgg                                                                                                                   <  1:858253/150‑1 (MQ=255)
                         ttCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGt                                                                                                                 >  1:1213455/1‑150 (MQ=255)
                          tCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGTGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTc                                                                                                                <  4:613460/150‑1 (MQ=255)
                          tCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTc                                                                                                                >  2:558575/1‑150 (MQ=255)
                            cTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCAc                                                                                                              <  2:490160/150‑1 (MQ=255)
                              gAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCt                                                                                                            >  4:204574/1‑150 (MQ=255)
                               aaCGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCtt                                                                                                           >  1:638726/1‑150 (MQ=255)
                               aaCGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCGCAAGGGTCACCtt                                                                                                           <  2:722825/150‑1 (MQ=255)
                                aCGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTg                                                                                                          <  1:79846/150‑1 (MQ=255)
                                aCGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTg                                                                                                          <  3:536183/150‑1 (MQ=255)
                                aCGCCTTCTTTGGTGCTGTACTCATGGAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTg                                                                                                          >  2:711226/1‑150 (MQ=255)
                                   ccTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGgcg                                                                                                       <  4:447381/150‑1 (MQ=255)
                                    cttcttTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGgcgt                                                                                                      <  1:912905/150‑1 (MQ=255)
                                    cttcttTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGgcgt                                                                                                      <  1:914242/150‑1 (MQ=255)
                                     ttcttTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGgcgtg                                                                                                     >  4:546294/1‑150 (MQ=255)
                                         ttGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTc                                                                                                 <  3:705277/150‑1 (MQ=255)
                                           ggTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGa                                                                                               >  4:538467/1‑150 (MQ=255)
                                            ttGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGaa                                                                                              <  2:301917/149‑1 (MQ=255)
                                              gCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACt                                                                                            >  3:127240/1‑150 (MQ=255)
                                                tGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCa                                                                                          >  2:323659/1‑150 (MQ=255)
                                                tGTACTCATGCAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCa                                                                                          <  1:1134211/150‑1 (MQ=255)
                                                 gTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTTTGTCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCAc                                                                                         >  1:267983/1‑150 (MQ=255)
                                                  tACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGGATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACt                                                                                        >  4:24227/1‑150 (MQ=255)
                                                       aTGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTc                                                                                   >  2:985330/1‑150 (MQ=255)
                                                          tagtaCACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGat                                                                                <  2:1131855/150‑1 (MQ=255)
                                                          tagtaCACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGat                                                                                >  2:131646/1‑150 (MQ=255)
                                                           agtaCACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGata                                                                               >  3:491254/1‑150 (MQ=255)
                                                            gtaCACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGatat                                                                              >  1:1033150/1‑150 (MQ=255)
                                                              acacCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCa                                                                            >  3:186275/1‑150 (MQ=255)
                                                                acCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGTGTGCGTCGAACTCACTTGCTCGATATCAAc                                                                          >  3:643097/1‑150 (MQ=255)
                                                                acCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAAc                                                                          <  2:1033150/150‑1 (MQ=255)
                                                                  cATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCa                                                                        >  2:8938/1‑150 (MQ=255)
                                                                   aTCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAg                                                                       >  1:920307/1‑150 (MQ=255)
                                                                    tCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAgg                                                                      >  4:211595/1‑150 (MQ=255)
                                                                     cAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGgc                                                                     >  4:697084/1‑150 (MQ=255)
                                                                       aTCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGgcgc                                                                   <  3:25785/150‑1 (MQ=255)
                                                                        tctcAACCTCGGTCACCGCGCAACGCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATACCAACCAGGcgcg                                                                  >  4:156500/1‑150 (MQ=255)
                                                                        tctcAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGcgcg                                                                  >  3:664310/1‑150 (MQ=255)
                                                                        tctcAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGcgcg                                                                  >  1:361873/1‑150 (MQ=255)
                                                                           cAACCTCGGTCACCGCGCAACTCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGtt                                                               >  2:695187/1‑150 (MQ=255)
                                                                           cAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGtt                                                               >  3:321651/1‑150 (MQ=255)
                                                                            aaCCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGttt                                                              <  1:1022510/150‑1 (MQ=255)
                                                                            aaCCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGttt                                                              <  2:566659/150‑1 (MQ=255)
                                                                            aaCCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGttt                                                              <  1:634931/150‑1 (MQ=255)
                                                                               cTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAGGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAg                                                           <  3:441988/150‑1 (MQ=255)
                                                                                tCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGa                                                          <  4:321651/150‑1 (MQ=255)
                                                                                 cGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGaa                                                         <  2:145985/150‑1 (MQ=255)
                                                                                  ggTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGaaa                                                        <  1:1074428/150‑1 (MQ=255)
                                                                                   gTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAAc                                                       <  4:186275/150‑1 (MQ=255)
                                                                                    tCACCGGGCAACCGGGCATCGATGAGAGCAGGGTACGGCGCAGTGCGTTACCCATAGTGTGGCAAACGCCACGCACTAAAGGCTCAAGGTGCACCTTGGCGTGCGTCGAACTCCCTTGTGCGATATCGACCGGGCGCGGTGTTattgaca                                                      <  2:343269/150‑7 (MQ=255)
                                                                                       ccGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTg                                                   >  2:215808/1‑150 (MQ=255)
                                                                                        cgcgCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGt                                                  <  1:136911/150‑1 (MQ=255)
                                                                                        cgcgCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGt                                                  <  1:558575/150‑1 (MQ=255)
                                                                                        cgcgCAACACGGCATCGATGAGAGCAGGATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGAGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGt                                                  <  2:507833/150‑1 (MQ=255)
                                                                                          cgcAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTca                                                <  3:546294/150‑1 (MQ=255)
                                                                                            cAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTcaca                                              <  2:677841/150‑1 (MQ=255)
                                                                                              aCCCGGCATCGATGAGAGCAGAATAAGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGa                                            <  3:535836/150‑1 (MQ=255)
                                                                                               cccGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGGCACAGaa                                           <  2:1213455/150‑1 (MQ=255)
                                                                                                 cGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAAcc                                         >  2:498357/1‑150 (MQ=255)
                                                                                                     aTCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGc                                     <  1:8938/150‑1 (MQ=255)
                                                                                                         aTGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACtt                                                                                       >  3:317036/1‑96 (MQ=255)
                                                                                                         aTGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACtt                                                                                       <  4:317036/96‑1 (MQ=255)
                                                                                                          tGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCAttct                                >  2:223940/1‑148 (MQ=255)
                                                                                                           gagaGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATtgtg                               >  4:283872/1‑150 (MQ=255)
                                                                                                           gagaGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATtgtg                               >  3:451801/1‑150 (MQ=255)
                                                                                                           gagaGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATACCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATtgtg                               <  3:156500/150‑1 (MQ=255)
                                                                                                            agagCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTgtgt                              <  2:36549/150‑1 (MQ=255)
                                                                                                                 aGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTca                                                <  4:472098/127‑1 (MQ=255)
                                                                                                                 aGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTca                                                >  3:472098/1‑127 (MQ=255)
                                                                                                                    aTACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCtt                      <  4:491254/150‑1 (MQ=255)
                                                                                                                     tACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCttt                     <  4:297520/150‑1 (MQ=255)
                                                                                                                     tACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCttt                     <  1:61118/150‑1 (MQ=255)
                                                                                                                      aCGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTg                    <  1:711226/150‑1 (MQ=255)
                                                                                                                      aCGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGCCACAGAACCCTGCATTGTGTCCTCTCTTTg                    <  4:127240/150‑1 (MQ=255)
                                                                                                                        ggCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGt                  >  3:512724/1‑150 (MQ=255)
                                                                                                                         gcgcAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTa                 >  4:733637/1‑150 (MQ=255)
                                                                                                                             aGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCttt                     <  3:725532/142‑1 (MQ=255)
                                                                                                                             aGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCttt                     >  4:725532/1‑142 (MQ=255)
                                                                                                                               tGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGgc                                                                                                        <  1:318410/57‑1 (MQ=255)
                                                                                                                               tGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGgc                                                                                                        >  2:318410/1‑57 (MQ=255)
                                                                                                                                 cGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCTATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCtt         <  3:204574/150‑1 (MQ=255)
                                                                                                                                     aCCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACt     >  2:299614/1‑150 (MQ=255)
                                                                                                                                      cccAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTgtgt                              >  2:101279/1‑124 (MQ=255)
                                                                                                                                      cccAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTgtgt                              <  1:101279/124‑1 (MQ=255)
                                                                                                                                       ccAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTg   <  1:985330/150‑1 (MQ=255)
                                                                                                                                       ccAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTg   >  3:67609/1‑150 (MQ=255)
                                                                                                                                       ccAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTg   >  1:501623/1‑150 (MQ=255)
                                                                                                                                        cAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTgg  <  3:733637/150‑1 (MQ=255)
                                                                                                                                        cAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTgg  <  4:451801/150‑1 (MQ=255)
                                                                                                                                        cAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTgg  >  2:194658/1‑150 (MQ=255)
                                                                                                                                        cAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTCTCCTCTCTTTGGTACTAAGCTTTACTTgg  >  4:717043/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                          
GGTTGAGCAGGATTTCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGG  >  NC_000913_3_pae_pgi/3440777‑3441062

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: