Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,606,516 G→T 91.7% D420Y (GAC→TAC)  hyfR → hydrogenase‑4 transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,606,5160GT91.7% 976.8 / 75.7 334D420Y (GAC→TAC) hyfRhydrogenase‑4 transcriptional activator
Reads supporting (aligned to +/- strand):  ref base G (15/13);  new base T (144/162);  total (159/175)
Fisher's exact test for biased strand distribution p-value = 5.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

Only 100 of 334 total aligned reads displayed.

GTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATGACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGAAATTGGCGATCTGCCGTTAGAACTTCAGCCTAAACTGCTGCGCGTATTGCAGGAGCGGGAGATTG  >  CP009273/2606369‑2606663
                                                                                                                                                   |                                                                                                                                                   
 tgtTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACa                                                                                                                                                  <  1:1520042/150‑1 (MQ=255)
 tgtTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACa                                                                                                                                                  >  2:2007092/1‑150 (MQ=255)
  gtTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACaa                                                                                                                                                 >  2:434564/1‑150 (MQ=255)
   ttGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACaaa                                                                                                                                                <  1:3666491/150‑1 (MQ=255)
     gATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAgg                                                                                                                                              <  1:2265402/150‑1 (MQ=255)
      aTTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAggg                                                                                                                                             >  1:549378/1‑150 (MQ=255)
         tgtgGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCg                                                                                                                                          >  1:2631592/1‑150 (MQ=255)
          gtggtgAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGt                                                                                                                                         >  2:3352444/1‑150 (MQ=255)
          gtggtgAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATGACAAAGGGGCGt                                                                                                                                         <  1:441716/150‑1 (MQ=255)
           tggtgAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGtt                                                                                                                                        >  1:4738829/1‑150 (MQ=255)
              tgAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTAc                                                                                                                                     <  2:1321955/150‑1 (MQ=255)
              tgAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTAc                                                                                                                                     >  2:2935436/1‑150 (MQ=255)
              tgAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTAc                                                                                                                                     >  1:876317/1‑150 (MQ=255)
               gAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACt                                                                                                                                    >  2:2967919/1‑150 (MQ=255)
                aaaCCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTg                                                                                                                                   <  1:4481374/150‑1 (MQ=255)
                  aCCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGt                                                                                                                                 >  1:1168417/1‑150 (MQ=255)
                   ccGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTg                                                                                                                                >  1:4904699/1‑150 (MQ=255)
                     ggAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCg                                                                                                                              >  1:3318545/1‑150 (MQ=255)
                       aaCTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGAt                                                                                                                            >  1:1604971/1‑150 (MQ=255)
                        aCTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGAtt                                                                                                                           >  1:3142259/1‑150 (MQ=255)
                         cTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTa                                                                                                                          <  1:2584927/150‑1 (MQ=255)
                         cTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTa                                                                                                                          <  1:1844983/150‑1 (MQ=255)
                            gCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATa                                                                                                                       <  2:4035648/150‑1 (MQ=255)
                              aaaGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATAcc                                                                                                                     <  1:3355385/150‑1 (MQ=255)
                              aaaGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATAcc                                                                                                                     >  2:5138461/1‑150 (MQ=255)
                               aaGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATAccc                                                                                                                    <  1:1285696/150‑1 (MQ=255)
                                  agGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATc                                                                                                                 <  2:3492536/150‑1 (MQ=255)
                                      gATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTgg                                                                                                             <  1:4223957/150‑1 (MQ=255)
                                       aTCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGt                                                                                                            <  2:1881987/150‑1 (MQ=255)
                                       aTCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGt                                                                                                            >  1:4252171/1‑150 (MQ=255)
                                        tCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTc                                                                                                           <  1:2960711/150‑1 (MQ=255)
                                          gCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGt                                                                                                         <  2:2152015/150‑1 (MQ=255)
                                            cAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGttt                                                                                                       >  2:279623/1‑150 (MQ=255)
                                             agagCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGtttt                                                                                                      >  2:3821649/1‑150 (MQ=255)
                                               agCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATGACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGa                                                                                                    >  2:4483209/1‑150 (MQ=255)
                                                gCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGaa                                                                                                   >  1:2140412/1‑150 (MQ=255)
                                                  gATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAAt                                                                                                 >  2:1130976/1‑150 (MQ=255)
                                                  gATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATGACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAAt                                                                                                 >  2:2079751/1‑150 (MQ=255)
                                                      cATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGcc                                                                                             >  1:4058234/1‑150 (MQ=255)
                                                               aGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATggcggc                                                                                    <  2:412837/150‑1 (MQ=255)
                                                               aGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATggcggc                                                                                    <  2:4945471/150‑1 (MQ=255)
                                                               aGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATggcggc                                                                                    <  2:2232592/150‑1 (MQ=255)
                                                                gCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCa                                                                                   >  2:1510578/1‑150 (MQ=255)
                                                                gCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATGACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCa                                                                                   <  1:4634128/150‑1 (MQ=255)
                                                                 cccGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCAc                                                                                  <  2:2212804/150‑1 (MQ=255)
                                                                  ccGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACg                                                                                 <  1:3409504/150‑1 (MQ=255)
                                                                  ccGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACg                                                                                 >  1:2892548/1‑150 (MQ=255)
                                                                  ccGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACg                                                                                 <  2:5057224/150‑1 (MQ=255)
                                                                   cgcgacgcgaCAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACgt                                                                                <  1:938787/150‑1 (MQ=255)
                                                                     cgacgcgacAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACgttt                                                                              <  1:1603574/150‑2 (MQ=255)
                                                                     cgacgcgacAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACgttg                                                                              <  1:5027104/150‑1 (MQ=255)
                                                                      gacgcgacAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACgttgt                                                                             >  2:4205604/1‑150 (MQ=255)
                                                                       acgcgacAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACgttgtt                                                                            <  2:3541531/150‑1 (MQ=255)
                                                                       acgcgacAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACgttgtt                                                                            <  2:3925953/150‑1 (MQ=255)
                                                                        cgcgacAAGCCGCTGGTCAAAATCGACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACgttgttc                                                                           <  2:1776201/150‑2 (MQ=255)
                                                                          cgacAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTtct                                                                         <  1:4171943/150‑1 (MQ=255)
                                                                          cgacAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATGACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTtct                                                                         >  2:1703216/1‑150 (MQ=255)
                                                                           gacAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTtctc                                                                        <  2:3308891/150‑1 (MQ=255)
                                                                                gCCGCTGGTCAAAATCAACTGCGCGGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGa                                                                   <  2:3293746/150‑1 (MQ=255)
                                                                                 ccGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGaa                                                                  >  1:72632/1‑150 (MQ=255)
                                                                                  cGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGaaa                                                                 >  2:133021/1‑150 (MQ=255)
                                                                                   gCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGAAAt                                                                <  1:2743744/150‑1 (MQ=255)
                                                                                    cTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGAAAtt                                                               <  1:639979/150‑1 (MQ=255)
                                                                                     tGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGAAATTg                                                              >  2:58333/1‑150 (MQ=255)
                                                                                         cAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGAAATTGGCGa                                                          <  1:1941835/150‑1 (MQ=255)
                                                                                           aaaTCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGAAATTGGCGATc                                                        <  1:3904944/150‑1 (MQ=255)
                                                                                              tCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGAAATTGGCGATCTGc                                                     >  1:3486031/1‑150 (MQ=255)
                                                                                               cAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGAAATTGGCGATCTGcc                                                    >  2:2954984/1‑150 (MQ=255)
                                                                                                   tGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGAAATTGGCGATCTGCCGTTa                                                <  2:67427/150‑1 (MQ=255)
                                                                                                   tGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGAAATTGGCGATCTGCCGTTa                                                <  1:3173150/150‑1 (MQ=255)
                                                                                                       cTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGAAATTGGCGATCTGCCGTTAGAAc                                            <  1:4311532/150‑1 (MQ=255)
                                                                                                          ccATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGAAATTGGCGATCTGCCGTTAGAACTTc                                         >  2:1717055/1‑150 (MQ=255)
                                                                                                          ccATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATGACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGAAATTGGCGATCTGCCGTTAGAACTTc                                         >  1:3207789/1‑150 (MQ=255)
                                                                                                          ccATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATGACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGAAATTGGCGATCTGCCGTTAGAACTTc                                         >  2:3714129/1‑150 (MQ=255)
                                                                                                             tCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGAAATTGGCGATCTGCCGTTAGAACTTCAGc                                      <  1:307187/150‑1 (MQ=255)
                                                                                                               cccGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGAAATTGGCGATCTGCCGTTAGAACTTCAGCCt                                    <  2:619337/150‑1 (MQ=255)
                                                                                                                ccGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGAAATTGGCGATCTGCCGTTAGAACTTCAGCCTa                                   <  1:991072/150‑1 (MQ=255)
                                                                                                                  gCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGAAATTGGCGATCTGCCGTTAGAACTTCAGCCTaaa                                 <  2:3288509/150‑1 (MQ=255)
                                                                                                                    cAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGAAATTGGCGATCTGCCGTTAGAACTTCAGCCTAAAct                               >  1:3550874/1‑150 (MQ=255)
                                                                                                                    cAGTCTTCTGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGAAATTGGCGATCTGCCGTTAGAACTTCAGCCTAAAct                               <  2:1063043/150‑1 (MQ=255)
                                                                                                                       tcttctGGAAAGTGAGTTATTCGGTCATGACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGAAATTGGCGATCTGCCGTTAGAACTTCAGCCTAAActgct                            <  1:1689628/150‑1 (MQ=255)
                                                                                                                        cttctGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGAAATTGGCGATCTGCCGTTAGAACTTCAGCCTAAActgctg                           <  1:4005450/150‑1 (MQ=255)
                                                                                                                         ttctGGAAAGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGAAATTGGCGATCTGCCGTTAGAACTTCAGCCTAAACtgctgc                          >  2:2724563/1‑150 (MQ=255)
                                                                                                                                aaGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGAAATTGGCGATCTGCCGTTAGAACTTCAGCCTAAACTGCTGCGCGTAtt                   >  1:153161/1‑150 (MQ=255)
                                                                                                                                aaGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGAAATTGGCGATCTGCCGTTAGAACTTCAGCCTAAACTGCTGCGCGTAtt                   >  2:4218719/1‑150 (MQ=255)
                                                                                                                                aaGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGAAATTGGCGATCTGCCGTTAGAACTTCAGCCTAAACTGCTGCGCGTAtt                   <  2:1956183/150‑1 (MQ=255)
                                                                                                                                 aGTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGAAATTGGCGATCTGCCGTTAGAACTTCAGCCTAAACTGCTGCGCGTATTg                  >  2:3126820/1‑150 (MQ=255)
                                                                                                                                  gTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGAAATTGGCGATCTGCCGTTAGAACTTCAGCCTAAACTGCTGCGCGTATTGc                 >  2:161919/1‑150 (MQ=255)
                                                                                                                                  gTGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGAAATTGGCGATCTGCCGTTAGAACTTCAGCCTAAACTGCTGCGCGTATTGc                 >  2:2591653/1‑150 (MQ=255)
                                                                                                                                   tGAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGAAATTGGCGATCTGCCGTTAGAACTTCAGCCTAAACTGCTGCGCGTATTGCa                <  1:142451/150‑1 (MQ=255)
                                                                                                                                    gAGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGAAATTGGCGATCTGCCGTTAGAACTTCAGCCTAAACTGCTGCGCGTATTGCAg               <  1:4701468/150‑1 (MQ=255)
                                                                                                                                     aGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGAAATTGGCGATCTGCCGTTAGAACTTCAGCCTAAACTGCTGCGCGTATTGCAgg              >  2:4680690/1‑150 (MQ=255)
                                                                                                                                     aGTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGACGGCACGTTGTTTCTCGATGAAATTGGCGATCTGCCGTTAGAACTTCAGCCTAAACTGCTGCGCGTATTGCAgg              >  1:3682435/1‑150 (MQ=255)
                                                                                                                                      gTTATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGAAATTGGCGATCTGCCGTTAGAACTTCAGCCTAAACTGCTGCGCGTATTGCAGGa             <  1:3529059/150‑1 (MQ=255)
                                                                                                                                       ttATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGAAATTGGCGATCTGCCGTTAGAACTTCAGCCTAAACTGCTGCGCGTATTGCAGGAg            >  1:1006997/1‑150 (MQ=255)
                                                                                                                                        tATTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGAAATTGGCGATCTGCCGTTAGAACTTCAGCCTAAACTGCTGCGCGTATTGCAGGAGc           <  2:2817002/150‑1 (MQ=255)
                                                                                                                                         aTTCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGAAATTGGCGATCTGCCGTTAGAACTTCAGCCTAAACTGCTGCGCGTATTGCAGGAGCg          >  1:5077023/1‑150 (MQ=255)
                                                                                                                                          ttCGGTCATTACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGAAATTGGCGATCTGCCGTTAGAACTTCAGCCTAAACTGCTGCGCGTATTGCAGGAGCgg         <  2:4576030/150‑1 (MQ=255)
                                                                                                                                                 aTGACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGAAATTGGCGATCTGCCGTTAGAACTTCAGCCTAAACTGCTGCGCGTATTGCAGGAGCGGGAGATTg  >  2:3483037/1‑150 (MQ=255)
                                                                                                                                                 aTGACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGAAATTGGCGATCTGACGTTAGAACTTCATCCTAAACTGCTGCGCGTATTGCAGGAGCGGGAGATTg  >  2:3196253/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GTGTTGATTTGTGGTGAAACCGGAACTGGCAAAGAGGTGATCGCCAGAGCGATCCATCAACTTAGCCCGCGACGCGACAAGCCGCTGGTCAAAATCAACTGCGCTGCCATCCCCGCCAGTCTTCTGGAAAGTGAGTTATTCGGTCATGACAAAGGGGCGTTTACTGGTGCGATTAATACCCATCGTGGTCGTTTTGAAATTGCCGATGGCGGCACGTTGTTTCTCGATGAAATTGGCGATCTGCCGTTAGAACTTCAGCCTAAACTGCTGCGCGTATTGCAGGAGCGGGAGATTG  >  CP009273/2606369‑2606663

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: