Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,443,196 A→T 15.7% Q474L (CAG→CTG)  feaB → phenylacetaldehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,443,1960AT15.7% 660.2 / 104.6 242Q474L (CAG→CTG) feaBphenylacetaldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (113/91);  new base T (23/15);  total (136/106)
Fisher's exact test for biased strand distribution p-value = 5.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 249 total aligned reads displayed.

AAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTC  >  CP009273/1443047‑1443345
                                                                                                                                                     |                                                                                                                                                     
aaaCGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCa                                                                                                                                                       <  1:1235852/150‑1 (MQ=255)
aaaCGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCa                                                                                                                                                       <  1:3482275/150‑1 (MQ=255)
  aCGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGt                                                                                                                                                     >  2:1497841/1‑150 (MQ=255)
    gACACGGAATATGGCTTAACTGCCTGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCa                                                                                                                                                   <  2:3422641/150‑1 (MQ=255)
    gACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCa                                                                                                                                                   <  1:4146465/150‑1 (MQ=255)
      cacGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCAGACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAgg                                                                                                                                                 <  2:801214/150‑1 (MQ=255)
        cGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGaa                                                                                                                                               >  1:197160/1‑150 (MQ=255)
          gAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACg                                                                                                                                             <  2:3310842/150‑1 (MQ=255)
              atGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGcc                                                                                                                                         >  2:2201811/1‑150 (MQ=255)
                   ttAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGAt                                                                                                                                    >  1:1728411/1‑150 (MQ=255)
                     aaCTGCCAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGAttt                                                                                                                                  >  1:991389/1‑150 (MQ=255)
                       ccgCCCGTGTCTGGACCCAACACCTCACCCACGCTCTCGAAGAGAGCGATCGCTTACCGGCAGGGACGGTGTGCGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTg                                                                                                                                <  2:2587215/148‑1 (MQ=255)
                        tGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTgg                                                                                                                               >  2:3412619/1‑150 (MQ=255)
                                tCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGAct                                                                                                                       >  1:1714645/1‑150 (MQ=255)
                                tCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGAct                                                                                                                       >  1:1716158/1‑150 (MQ=255)
                                 cTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGActg                                                                                                                      >  1:1218123/1‑150 (MQ=255)
                                 cTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGActg                                                                                                                      <  2:2291089/150‑1 (MQ=255)
                                  tGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGActgg                                                                                                                     <  1:81541/150‑1 (MQ=255)
                                  tGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGActgg                                                                                                                     <  1:583146/150‑1 (MQ=255)
                                  tGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGActgg                                                                                                                     >  1:3207664/1‑150 (MQ=255)
                                   ggACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGActggc                                                                                                                    <  2:1522954/150‑1 (MQ=255)
                                     aCGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGActggctg                                                                                                                  <  2:1447666/150‑1 (MQ=255)
                                        cAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGAc                                                                                                               <  2:974206/150‑1 (MQ=255)
                                          aaaTCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGGGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACgg                                                                                                             >  2:3246051/1‑150 (MQ=255)
                                          aaaTATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACgg                                                                                                             >  1:2207353/1‑150 (MQ=255)
                                           caTCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGt                                                                                                            <  2:3490462/149‑1 (MQ=255)
                                           aaTCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGt                                                                                                            <  1:2665046/150‑1 (MQ=255)
                                           aaTCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGt                                                                                                            <  2:3490424/150‑1 (MQ=255)
                                                       gCTCTGGACTATAGCGATCGCTTACAGTCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACt                                                                                                <  2:4410061/150‑1 (MQ=255)
                                                       gCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACt                                                                                                <  2:2170108/150‑1 (MQ=255)
                                                        ctctGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTa                                                                                               >  2:3481666/1‑150 (MQ=255)
                                                          ctGGAATATAGAGATCGCTTCCAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAg                                                                                             <  2:4014706/150‑1 (MQ=255)
                                                           tGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGt                                                                                            <  2:1890892/150‑1 (MQ=255)
                                                            ggAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTc                                                                                           <  2:2508960/150‑1 (MQ=255)
                                                            ggAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTc                                                                                           <  2:3584682/150‑1 (MQ=255)
                                                            ggAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTc                                                                                           <  2:3584155/150‑1 (MQ=255)
                                                            ggAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTc                                                                                           <  1:4149627/150‑1 (MQ=255)
                                                             gAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCg                                                                                          <  2:4158823/150‑1 (MQ=255)
                                                               atatAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGgt                                                                                        >  2:348194/1‑150 (MQ=255)
                                                                tataGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGgtg                                                                                       >  2:1291211/1‑150 (MQ=255)
                                                                 ataGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGgtgt                                                                                      >  2:645031/1‑150 (MQ=255)
                                                                   aGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGgtgtgt                                                                                    >  1:3813804/1‑150 (MQ=255)
                                                                   aGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGgtgtgt                                                                                    >  2:3226995/1‑150 (MQ=255)
                                                                   aGCGATCGCTGACAGGCAGGGACGGTGTGGTTAAACAGCCTTACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTTATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGgtgtgt                                                                                    <  2:4156796/150‑1 (MQ=255)
                                                                    gCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGgtgtgtg                                                                                   >  2:2635302/1‑150 (MQ=255)
                                                                     cGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGgtgtgtgt                                                                                  >  1:1864745/1‑150 (MQ=255)
                                                                           cTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTa                                                                            <  2:640150/150‑1 (MQ=255)
                                                                           cTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTa                                                                            >  1:2072747/1‑150 (MQ=255)
                                                                               caggcaggGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTCTTaa                                                                        >  1:1376115/1‑150 (MQ=255)
                                                                                    aggGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTgg                                                                   >  1:542517/1‑150 (MQ=255)
                                                                                      ggACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGtt                                                                 <  2:2666123/150‑1 (MQ=255)
                                                                                       gACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTc                                                                <  1:1442524/150‑1 (MQ=255)
                                                                                       gACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTc                                                                <  1:3708393/150‑1 (MQ=255)
                                                                                           gtgtGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTc                                                            <  1:3031058/150‑1 (MQ=255)
                                                                                            tgtgGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCa                                                           >  2:3160527/1‑150 (MQ=255)
                                                                                            tgtgGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCa                                                           >  2:3160662/1‑150 (MQ=255)
                                                                                             gtgGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCAt                                                          >  1:1240979/1‑150 (MQ=255)
                                                                                               gggTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATaa                                                        >  1:3467123/1‑150 (MQ=255)
                                                                                               gggTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATaa                                                        >  1:2601198/1‑150 (MQ=255)
                                                                                                 gTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGt                                                      <  1:1888626/150‑1 (MQ=255)
                                                                                                  tAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTa                                                     >  1:811379/1‑150 (MQ=255)
                                                                                                    aaCAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTaaa                                                   >  1:1337924/1‑150 (MQ=255)
                                                                                                       aGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCAGGTTCGCTCATAAGTaaaaaa                                                >  2:2914747/1‑150 (MQ=255)
                                                                                                         ccATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACg                                              <  1:629127/150‑1 (MQ=255)
                                                                                                           aTACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGc                                            >  1:483045/1‑150 (MQ=255)
                                                                                                           aTACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGc                                            >  1:3512942/1‑150 (MQ=255)
                                                                                                               cTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCt                                        >  2:3249088/1‑150 (MQ=255)
                                                                                                                 tAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTgg                                      >  2:4463712/1‑150 (MQ=255)
                                                                                                                  aaTTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGt                                     >  1:4360896/1‑150 (MQ=255)
                                                                                                                   aTTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTg                                    >  1:2820521/1‑150 (MQ=255)
                                                                                                                     tGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGcc                                  <  1:1610114/150‑1 (MQ=255)
                                                                                                                        cGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGtt                               >  1:156035/1‑150 (MQ=255)
                                                                                                                         gCTAGCTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGttt                              >  1:2087605/1‑150 (MQ=255)
                                                                                                                         gCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGAGTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGttt                              >  1:3942360/1‑150 (MQ=255)
                                                                                                                           tAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGttttt                            <  2:3610955/150‑1 (MQ=255)
                                                                                                                             aCTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGatttttt                          <  2:991651/150‑1 (MQ=255)
                                                                                                                              cTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTg                         >  1:3006355/1‑150 (MQ=255)
                                                                                                                              cTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTg                         >  1:3090842/1‑150 (MQ=255)
                                                                                                                              cTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTg                         >  1:926387/1‑150 (MQ=255)
                                                                                                                                tACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTc                       <  2:219158/150‑1 (MQ=255)
                                                                                                                                tACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTc                       >  1:1046063/1‑150 (MQ=255)
                                                                                                                                tACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTc                       <  1:4306659/150‑1 (MQ=255)
                                                                                                                                 aCCGTGTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCt                      <  2:2964493/150‑1 (MQ=255)
                                                                                                                                  ccGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTg                     <  1:2201179/150‑1 (MQ=255)
                                                                                                                                  ccGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTg                     <  2:197160/150‑1 (MQ=255)
                                                                                                                                   cGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGa                    <  1:3588645/150‑1 (MQ=255)
                                                                                                                                    gTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGaa                   >  1:3823770/1‑150 (MQ=255)
                                                                                                                                     tttGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTCTTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGaaa                  >  2:1090667/1‑150 (MQ=255)
                                                                                                                                     tttGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGaaa                  <  2:1338521/150‑1 (MQ=255)
                                                                                                                                      ttGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAAc                 >  2:900612/1‑150 (MQ=255)
                                                                                                                                        ggtggGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACaa               >  1:2355515/1‑150 (MQ=255)
                                                                                                                                           gggATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCa            >  1:4363391/1‑150 (MQ=255)
                                                                                                                                             gATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACt          <  1:165281/150‑1 (MQ=255)
                                                                                                                                               tGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTa        >  1:4484643/1‑150 (MQ=255)
                                                                                                                                                gAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTAt       >  1:2399326/1‑150 (MQ=255)
                                                                                                                                                  aGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTAtct     <  2:1397570/150‑1 (MQ=255)
                                                                                                                                                   gCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTAtctt    <  2:4083917/150‑1 (MQ=255)
                                                                                                                                                    cAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTAtcttt   >  2:1696289/1‑150 (MQ=255)
                                                                                                                                                    cAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGGTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTAtcttt   >  2:3740864/1‑150 (MQ=255)
                                                                                                                                                     aGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTAtctttc  <  2:1766232/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTC  >  CP009273/1443047‑1443345

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